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Abstract 

Moving from Beisert-Staudacher equations, the complete set of Asymptotic Bethe Ansatz 
equations and 5-matrix for the excitations over the GKP vacuum is found. The resulting 
model on this new vacuum is an integrable spin chain of length i? = 2 In s (s = spin) with 
particle rapidities as inhomogeneities, two (purely transmitting) defects and SU{A) (residual 
R-)symmetry. The non-trivial dynamics of AA = 4 SYM appears in elaborated dressing factors 
of the 2D two-particle scattering factors, all depending on the ’fundamental’ one between two 
scalar excitations. From scattering factors we determine bound states. In particular, we 
study the strong coupling limit, in the non-perturbative, perturbative and giant hole regimes. 
Eventually, from these scattering data we construct the 4D pentagon transition amplitudes 
(perturbative regime). In this manner, we detail the multi-particle contributions (flux tube) 
to the MHV gluon scattering amplitudes/Wilson loops (OPE or BSV series) and re-sum them 
to the Thermodynamic Bubble Ansatz. 
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1 Introduction 


The study of the energy of the excitations on a suitably chosen vacuum state is a problem which 
is common to very many physical theories. It often happens that most intriguing excitations arise 
over a vacuum state which is an intricate superposition of ’basic’ states, i.e. a sort of Fermi sea 
of interacting ’pseudoparticles’. In general, this vacuum may be dubbed antiferromagnetic as the 
prototypical example in the realm of integrable models is the antiferromagnetic vacuum state of 
the Heisenberg spin chain. In their turn, important excitations on it are called spinons or solitons, 
whilst magnons are the (pseudo)particles forming the see on the ferromagnetic vacuum. In an easy 
Bethe Ansatz perspective [1], spinons may appear as holes in a distribution of a large number of 
real Bethe roots. As a consequence, these holes are constrained by quantisation conditions for their 
rapidities, which may anew be seen as Bethe (-Yang) equations for these new ’fundamental’ parti¬ 
cles. Of course, we expect this phenomenon to be of non-perturbative nature, so that integrability 
is the right realm to exploit it. 

A similar, but obviously much richer situation, arises in the framework of Beisert-Staudacher 
Asymptotic Bethe Ansatz (ABA) equations determining, via a specihc root configuration, the 
anomalous dimension (energy) of the single trace helds in planar A/" = 4 SYM [2]. In this context 
one can choose as ’antiferromagnetic’ vacuum the configuration which contains a large number, 
s, of type-4 roots and which describes, up to wrapping corrections [3], high spin (= s) twist two 
operators, namely, sketchily, 

0 = TrZDlZ + ... , (1.1) 

where Z is one of the three (complex) scalars of the theory. In fact, this is likely the ’simplest’ 
example of Wilson twist operator. It belongs to the paradigmatic sl(2) sector of scalar operators, 
which are made up of only one (out of three) complex scalar Z and the (light-cone) covariant 
derivative Z1+, so enjoying the sketchy form 

Tr(D;Z^)+ ...., (1.2) 

where dots stand for permutations. Built up in this selected way they result to be perturbatively 
closed under renormalisation, so forming a sector. These composite single trace operator have of 
course Lorentz spin s and twist (or length in the ferromagnetic/half-BPS vacuum perspective) L, 
with minimum value L = 2 for which (a descendant^ of) the Gubser-Klebanov-Polyakov (GKP) 
’vacuum’ solution is realised [4]. Also, the AdS/GFT correspondence [5] relates an operator (1.2) 
to a spinning folded closed strings on AdSs x spacetime, with angular momenta s/y/X and Lj\f\ 
on each space respectively, the’t Hooft coupling in the multi-color A'), —>■ oo (planar) regime 

= ^cQym 'i 9^ = ; (1-3) 

being connected to the string tension T = ^ [4, 6]. On the other hand we may think of the 
operators (1.2) as obtained from the GKP vacuum (1.1) by adding scalar excitations on top of 

^In this case, the spin may be shifted by a finite amount which does not affect our analysis and results at high 
spin. 
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it. Of course, when L > 2 we can realise states with different energies, at fixed L, and typically 
the minimal energy has been more extensively studied for ’large size’ s —?■ +oo. In particular, 
the minimal anomalous dimension of (1.2) has been proven to enjoy at one loop the same leading 
behaviour ~ In s at high spin (and fixed L) [7], as in the opposite string regime (strong coupling) [4]. 
Later on, the coefficient of this term^ was obtained at all loops from the solution of a linear integral 
equation directly derived from the Beisert-Staudacher equations via the root density approach [10]. 
In very brief summary, as computed in [11], the high spin (asymptotic) expansion (at fixed g and 
L) enjoys the peculiar form 

CX) 

' 1 { 9 ,s,L) = f (g) Ins + f,i{g,L) + (Ins)"" + O ((lns)/s) , (1.4) 

n=l 

in inverse integer powers of the size^ i? ~ Ins, except the sub-leading (lns)° contribution fsi{g,L) 
(defect contribution). The latter, which reduces to the so-called virtual scaling function for L = 2, 
has been captured in [12] by a Non-Linear Integral Equation (NLIE) and in [13] by a linear integral 
equation (by means of which explicit strong coupling expansions can be performed [14], along the 
lines of those for the cusp [15]). Up to this order, we can be sure that this expansion enjoys the 
same form at all perturbative orders in QCD, or its Mellin transform, i.e. the evolution kernels [16]. 
Moreover, in the supersymmetric case similar linear integral equations hold for all the coefficients 
in (1.4) [11] and also for the next order O ((lns)/s) [17], and all these, - importantly the first two 
f{g) and fsi{g,L), - are now believed to be exactly given by the ABA without wrapping^, also 
thanks to these recent studies. 

The latter were focussed on the same scalar, Z, added to (1.1), but we can generalise to the 
other helds: indeed, elementary one-particle excitations may correspond to inserting one of the 
other helds i.e., besides the other two scalars, a gauge held (gluon) or a Fermi held (gaugino)^. In 
other words, they are the lowest twist (=three) operators/states with the form 

Oi-parucie = Tr ZDl-^'ipD^Z + ... , (1.5) 

where (p = Z,W,X, the scalars, or ip = F+±, F+_l, the two components of the gauge held, or 

ip = 4/_|_, 'L_|_, the 4-1-4 (anti-)fermions, respectively. Besides the energy, one can determine also 

the momentum of an operator through the Beisert-Staudacher ABA equations. Along this line, 
the one-particle dispersion relations of the excitations (1.5) have been receiving much attention in 
the diherent coupling regimes (c/. for instance [19] and references therein); but recently they have 
been summarised, corrected and put forward in an illuminating work by Basso [20] (also reference 
therein). 

^This is the so-called universal scaling function, f{g), which does not depend on L and equals twice the cusp 
anomalous dimension (renormalisation divergence [8]) of a light-like Wilson cusp, as in QCD [9]. 

^In fact, it is consistent with the length of the long classical string R ^ ln(s/-\/A) [4, 6]. 

^For instance in [18] wrapping corrections to ABA start to contribute at order e~^ = = 1/s^, inducing 

to think of a factor 2 in the size of the folded string i? = 2 In s -|- .... 

^Notice that in the half-BPS vacuum description this state would belong to a longer spin chain of length L = 3. 
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On the same footing, we start wondering in [21] about the scattering 5'-matrix which may be 
attached to the two-particle states (of, at least, twist-4) 

02-parUcles = Tv Z LfiD^ ip 2 D^ Z + . . . , (1-6) 

where </?i and (p 2 niay be any general elementary local field as (p in (1.5), whereas in [21] we confined 
our attention to the peculiar (c/. below) case pi = (p 2 = Z. In fact, as argued above, we expect 
the Beisert-Staudacher quantisation conditions to give correct results at leading Ins and next to 
leading order (lns)°. And then, regarding i? ~ Ins as the size of the system, these orders are 
exactly the ones we need to write down 2D (many-particle) scattering amplitudes, i.e. (on-shell) 
quantisation conditions, for rapidities of excitations on the GKP vacuum. Generalising to all the 
other scalars, [22] have deduced the entire AO(6) scattering, while we have computed in [23] all the 
gf-depending scalar factors of the different scattering channels, neglecting the SU (4) representation 
structure. 

Moving from this lack, we shall make here our analysis deeper, by computing explicitly the 
matrix structures of the different A17(4) representations carried by the ’elementary’ particles and 
by their bound states. We will not only consider the two-body scattering, but also in general the 
multi-particle 2D scattering amplitudes. As a byproduct we will see a well know characterisation of 
integrable theories, namely the elasticity and factorisation, i.e. the determination of many-particle 
scattering by the two-particle one. Besides the traditional name of Bethe-Yang equations, we can 
call these quantisation conditions Asymptotic Bethe Ansatz equations as well, but now the term 
’Asymptotic’ refers to the new length ~ Ins, which measures the validity of the equations (and 
to the ’new’ vacuum). More precisely, from the BMN (ferromagnetic) vacuum [24] (no roots) we 
will switch on, in the Beisert-Staudacher equations, the configurations corresponding to the GKP 
(antiferromagnetic) vacuum and to all possible ’elementary’ excitations over the GKP vacuum; 
to accomplish this, we will be using the idea of converting many (Bethe) algebraic equations 
describing an excited state into few non-linear integral equations (NLIEs) [25, 26, 27, 12, 13]. In 
this way, we will obtain the quantisation conditions of all the ’elementary’ excitations over the 
GKP vacuum and show that the structure of these equations coincides with Bethe equations of 
a inhomogeneous spin chain of length i? = 2 In s with two identical (purely transmitting) defects 
and a Af/(4) symmetry in different representations (where the particle rapidities represent the 
inhomogeneities). Of course, the scalar pre-factors in front of the above At/(4) matrix structure 
are dependent on g and characteristic of the theory (and GKP vacuum). Nevertheless, we can 
express all in terms of the scalar-scalar one [23]. Moreover, we will discuss in many details the 
consequences of switching to a different vacuum which basically means that any elementary particle 
interacts with the sea of covariant derivatives namely the type-4 roots. For instance, the poles 
of the new 2D scattering factors of these particle imply the entrance of bound states thereof into 
the spectrum and then the existence of new scattering amplitudes for the latter particles. As 
anticipated, not only the 2D scattering amplitudes, but also many physical quantities assume 
novel expressions, as for instance the energy, momentum [20] and all the other conserved charges 
carried by a single elementary or composite excitation (c/. below). Furthermore, the scattering 
of any particle onto two defects arises, as anticipated in [22, 23], though they were absent in the 
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ferromagnetic setup, and they are likely to be associated to the two external holes (or tips of the 
GKP string). Together with the change of length, these are somehow unprecedented features in 
the theory of (quantum) integrable systems though their common origin can be traced back to 
the sl{2) spin chain (describing the one-scalar sector (1.2) at one-loop): nevertheless, we are used 
to insert the defect ab initio in the theory on the ferromagnetic vacuum and then hnding the 
anti-ferromagnetic dynamics with defect (possibly characterised by a different scattering factor) 
and the same length. 

As a consequence of this new ABA, also the exact Thermodynamic Bethe Ansatz (TBA) [28, 29] 
for the spectrum of anomalous dimensions as derived from its mirror version should ostensibly look 
very differently® than the usual one on the BMN vacuum [30], although they should give the same 
spectrum after all. Even more interestingly, a recent series of papers by Basso, Sever and Vieira 
(BSV) extended to all terms the operator product expansion of [41] and thus proposed a non- 
perturbative approach to 4D gluon scattering amplitudes/null polygonal Wilson loops (which are 
allegedly the same [31, 32, 33]) in W = 4 SYM, which relies on these 2D scattering factors as 
input data or building blocks [34, 35, 36, 37, 38]. In this perspective, light-like polygonal Wilson 
loops (WLs) can thought of as an inhnite sum over more fundamental polygons, namely square 
and pentagonal WLs, whose knowledge relies on the GKP scattering factors. By virtue of the 
AdS/GFT strong/weak duality, this superposition of pentagons and squares should lead, at large 
coupling to the classical string regime, namely the minimisation of the supersymmetric string 
action [31]. In general, this is a complicated problem of minimal area (string action) subtending 
a polygon living on the boundary of AdS^, and results in a set of non-linear coupled integral 
equations [39, 40, 41, 42]. For some still hidden reason, their form resembles that of a relativistic 
Thermodynamic Bethe (or Bubble, in this case!) Ansatz system whose free energy yields the area 
[29, 28]^. Instead, we wish in this paper to construct this TBA set-up by summing the inhnite 
BSV series and performing a saddle point evaluation. For this aim, we will perform a propaedeutic 
analysis of all the different strong coupling regimes. 

The article is organised according to the following plan. In section 2 we derive the ABA 
equations, hrst at one-loop as exemplifying case so to highlight all the relevant features, then for any 
value of the coupling. In section 3 the conserved charges of the excitations (on the GKP vacuum) 
are computed. In section 4 the strong coupling limit of the scattering factors is considered, in the 
different dynamical regimes, i. e. non-perturbative, perturbative and giant hole regimes. Section 5 
contains equivalent forms for the momentum associated to any elementary particle excitation, in 
particular that elaborated in [20]. Section 6 is a study of the strong coupling behaviour of (the 
scattering factor for) the spin chain defects. In section 7 the properties of the different kinds of 
particles under the SU(4) symmetry are taken into exam, so that in section 8 we are able to describe 
the structure of the overall S-matrix. In section 9 the so-called string hypothesis is used on the 
GKP ABA, in order to survey the on-shell states and, in particular, the bound states of elementary 
particles; for later purposes, an accurate study is devoted to the behaviour of bound states of 

®The attentive reader may guess many aspects of it from the form of string/stack solutions as reported in section 

9. 

^In the particular case of the hexagon WL [39] the system does coincide with the usual TBA one [43] 
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gluons at any coupling and fermion-antifermion=meson and its bound states which, instead, do 
appear only at the leading order of perturbative strong coupling regime, i.e. in the classical string 
theory. Section 10 computes all the string perturbative expressions of the pentagonal amplitudes 
(contributing to the BSV series). Finally, these inhnite contributions are summed up exactly in 
section 11 as for the hexagonal Wilson loop; the result is remarkably coinciding with the Yang- 
Yang functional and Thermodynamic Bubble Ansatz (TBA) equations for minimal area solution 
[39, 41, 40]. After the conclusions (section 12), several appendices follow. In appendices A and B 
the dehnition of the functions employed throughout the text, as well as some useful formulae and 
integrals are listed. Appendix C is a synopsis of scattering factors, displayed at arbitrary coupling, 
one-loop, strong coupling (mirror transformed too), while appendix D gives some details about 
their derivation. Finally, in appendix E all the ABA equations are listed. 


2 General equations 

2.1 Excitations 


The main aim of this section is to write Bethe equations describing ’elementary’ excitations over the 
long GKP string, in the more general case when H scalars {uh,h = Np large fermions, 

Np large antifermions, Nf small fermions, Nj small antifermions, Ng gauge helds and Ng 
gauge helds F+_l are present. 

In the notation of Beisert-Staudacher equations [2] the GKP vacuum is described by a large 
even number s of type-4 roots hlling the interval [—b, b] of the real axis, together with two external 
holes [44, 7]. In the large s limit b is approximated with s/2 and the positions of the two external 
holes with ±s/\/2: corrections to those estimates give rise to 0(l/s^) terms in the hnal Bethe 
equations, that we will neglect. 

It is a general fact that, in order to deal with a large number of Bethe roots, it is convenient 
to use their counting function Z^{y), which satishes a nonlinear integral equation [25, 26]. We 
found then natural to apply that strategy to the study of GKP vacuum and its excitations. In 
this approach scalar excitations, which are represented by holes in the distribution of type-4 roots 
in [—6, 5], are classihed by quantisation conditions for Z/^{v). The same function Z^{v) governs the 
interaction between roots with different flavour and scalars. 

Goming in specihc to the classihcation of the various excitations [20], we have already said 
that scalars are represented by holes in the distribution of type-4 roots. Large (small) fermions 
are described by Us-type (ui-type) roots and large (small) antifermions are described by Us-type 
(uy-type) roots. Rapidity of large fermions is the function xp{u^) = x^u^), where {g is related to 
the’t Hooft coupling A by A = Svr^^f^) 


x{u) 



u 


>2g^ 


( 2 . 1 ) 


with the arithmetic dehnition of the square root. Therefore, rapidity of large fermions satishes the 
inequality \xp\ > g/V2. On the other hand rapidity of small fermions is the function Xf{ui) = 





2 

2 ^^, with definition (2.1) for x{u) and, consequently, it is constrained by the inequality \xf \ < 
gj^- Changing uj, —)■ and Ui —?■ allows to describe large and small antifermions, respectively. 

Gauge fields F+x with rapidity correspond to stacks, 

U 2 ,j=u^j, U3j = u^j±i/2, j = l,...,Ng, (2.2) 

with real centres while gauge fields F+± with rapidity are described by stacks, 

uej=u^j, U5j =u^j±i/2, j = l,...,Ng, (2.3) 

with real centres 

We consider also the presence of isotopic roots, which do not carry momentum and energy, 
but take into account internal degrees of freedom, i.e. the residual SU (4) symmetry of the GKP 
vacuum. In specific, we have Ka roots Uaj of type U 2 , 

^a,j 1 j 1; •••; F^a ) (2"^) 


Kc roots Ucj of type Ug 

^c,j ) J Fq , ( 2 - 5 ) 

and Kf, stacks, 

U 4 .J ^b,j i ~ , ^b,j ) j 1) ■■■■> F}j , (2.6) 

with centers Ubj. 

We are now going to present our derivation of the full set of Bethe-Yang equations for excitations 
on the GKP vacuum. For excitations with rapidity Um belonging to the representation p of the 
symmetry group SU (4) of the GKP vacuum such equations will appear in the following form 


n 


'^m 4 “ lOiq • Wp 




'^q,k 


lag ■ Wp 


K, 


q 


F' 


u. 




nn 


^q^k P 4~ 


q ^q' 


laq ■ aq Uq^k 


^q'p 


laq ■ aq! 


1 = 


_ ^iRP {Urn^~\'‘^iD {Urn^ 


nn 

q k=l 


^q.k lOtq * W q 


Ur 


^q.k 


Ky.q * Wp 

^ ^ m'^m 


Um'^ 1 


. (2.7) 

( 2 . 8 ) 


where {aq} are the set of simple roots of S'f/(4), Uq^k are the isotopic roots associated and Wp the 
highest weight of the representation p. The structure of these equations agrees with the general 
pattern shown in [45]. While the first equation (2.7) comes from the symmetry properties of the 
vacuum, the second one (2.8) is a quantisation condition for the rapidity Um of an excitation moving 
in a one dimensional chain. Within this interpretation, R is given the meaning of the physical 
length of the chain, P{um) that of the momentum of an excitation with rapidity Um- The extra 
term 2D{um) in the exponent is interpreted as the effect of two purely transmitting [i.e. without 
reflection [46]) defects related to the tips of the GKP string. The rational factor in the right hand 
side of (2.8) takes into account the internal degrees of freedom: solving (2.7) one obtains Uq^k 
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in terms of Um'- plugging this result in the rational term in ( 2 . 8 ) one obtains, together with the 
products over the various S{um, Um'), the phase change due to the scattering between an excitation 
with rapidity Um and the other excitations with rapidity Um'- 

We will start with the one loop case, where all factors entering equations (2.8) are written in 
an explicit form, specihcally in terms of products of Euler Gamma functions. The general all loops 
case will appear as a technical complication of the one loop, since the building blocks of the various 
equations ( 2 . 8 ) will be obtained after solving a linear integral equation. 


2.2 Equations at one loop 


Scalars 


In order to show how our strategy works, we hrst concentrate on the one loop case. We 
start from the fourth of the Beisert-Staudacher equations, in the presence of a large number s of 
real type-4 roots, together with the general pattern of excitations and isotopic roots described in 
previous section. We remark that in the one loop case only large fermions and large antifermions 
are present: for uniformity of notations in this subsection we will denote large fermions rapidity 
xp with Up and large antifermions rapidity Xp with up. 

We introduce the counting function 


Z4{V) 


i K 4 

• r 1 9 ~ 1 t — t/ -r 

= iLm4 -h * > In r- 

7 - 1 - 7 ) — 


i — V + U4J 
^ + v ^ i + v - U 4 J 




t=l 

Upj 


Kb 

+ 2 i ^ In 
i=i 


+ V - Ub,j 
- V + UbJ 


Np 2 I 

W In 

^ f - U + Upj 


^ t - U -f Upj 


Na 


+ i In r— 


i + V — Ui 


Na 


I — V + U^j 


-I- i In 7 
i=i 


i + v 




l — V + Uj 


+ 


(2,9) 


where the sum up to i ^4 = s -|- 2Kb is a sum over s real type-4 roots and over the 2Kb complex 
type-4 roots contained in the stack (2.6). In terms of Z 4 the fourth of Beisert-Staudacher equations 
reads as 






( 2 . 10 ) 


In addition, in the large s limit the behaviour of Z 4 {v) is dominated by the second term in the 
right hand side, which implies that for v real Z'^{v) < 0. With this information we can prove that 
the length L is not independent of the total number of excitations. Indeed, it is widely known 
[44, 7] that condition (2.10) is satisfied on the real axis not only by the type-4 real roots, but also 
by -|- 2 real numbers, called holes. H holes are inside the interval [—h,h] {h = s/2 + 0(l/s^)) 
and, consequently, are often called ’internal’ or ’small’; the remaining two holes have positions 
±s /\/2 -|- 0(l/s^) and then are dubbed ’external’ or ’large’. Since Z'^{v) < 0 for n real, the 
difference between the extremal values on the real axis Z 4 {+oo) — Z 4 (—oo) has to count the total 
number of real roots and holes, i.e. 


Z 4 (-|-oo) — Z 4 (— 00 ) = —27r(s -|- if -|- 2 ) , ( 2 - 11 ) 

On the other hand, the definition (2.9) implies the asymptotic behaviours 
Z 4 (±cx>) = +7r{L + K 4 - 2Kb -Np-Np-Ng-Ng) = +tt{L + s - Np - Np-Ng-Ng) . (2.12) 
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Comparison between (2.11) and (2.12) gives the desired connexion 

L = H -\- 2 Np + Np + Ng + Ng , (2.13) 

between the length L and the total number of excitations. 

Relation (2.13) once plugged in the exponent of the right hand side of (2.10), provides a 
simplihcation of quantisation condition for type-4 roots and real holes. Restricting to holes, whose 
position we call Uh, h = 1 ,..., H, we get the compact formula (remember that s is always even) 

e"*^4K) ^ (-1)^-1. (2.14) 


After hxing these preliminary aspects, we come back to equation (2.9). In order to get manageable 
expressions, we convert the sum over real type-4 roots into an integral by means of the master 
equation [26]: 


real '^f(u 4 ,j) = -f 


i=i 


(^/2 


■/ 






where L^iy) = Imln[l -|- (—l)^e®^4h specialise formula (2.15) to our case and get 


(2.15) 


Ka 


Kb 


< ^ 7 J- 7 > — 7 / .. • < ^ 


l \ V — U4J 


{^ + v- Ubj){y + v- Ubj) {i + v - ■^){i + V + ^) 


V2^ 


H 


- ‘E'‘'7 


i - V + Uh 


f + oo 


+ 


dw 


h=l 


I + V — Uh 


TT 1 -f (u — w)‘ 


:lZi(w) - 2Li{w)] 


(2.16) 


Eventually, putting together (2.9) and (2.16), we find that 

/ + CXD 

dw G{v — w) L^iw ), 

•CXD 

where F[v),G{v) satisfy the linear integral equations 

IG, ■ o • 

I + f - Uhj f -V + UbJ 


(2.17) 


F{v) = + 




^-V + Ubjf+V-Ub,j ^ 


No 




i + V — m 


j=l 2 
9 Kp 


i + v - 

+ i 2^ In T — g + 


i - V + m 
j=i j 


j=l 2 


^ + U - Up 


Ni 


E l 9 4” n '^F,i 

In^-^ 

. , l-V+Upj 


+ 


H 


V + .n„ ElzilEAA + 

i — 


^+oo 


dw 


G{u — v) = 


h=l 

1 1 


t — V + Uh 


b+cx) 


TT 1 + {u — vy 


+ 


^-^ + 7=2^-^-V^2 
dw 1 


TT 1 + (u — w)' 


TT 1 + (m — w] 


;G{W — v) 


(2.18) 

:F(y>) . 

(2.19) 
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Equation (2.18) is solved exactly; however, we remember that the pattern of excitations discussed 
before holds only in the large spin limit. To be consistent with that, we have to use the large s 
asymptotic behaviour 


i In 




T[l-iv + i^ 


i In 




—4nln —= + 0(l/s^), 
v 2 


( 2 . 20 ) 


and write the hnal result for F{v) as 

r(i + ,v) 


F{v) = —iL\n 


r (i - iv) 


- Av\n—= + 1 } In —-: 

^/2 ^ TTI -9 


H 


r(l + in — iufi) 


h=l 


r(l - iv + iuh) 


+ 


^ { + v-Uhj T{l + iv-M) ^ T{l + iv-iu^A 

ln|-^ + i> ln- 7 --^-r^ + i> In^-K+ 2.21 

^ ^ T{l-iv + iu^.) ^ r(l-in + mj) 

^ ^ r(| + ^f 1 „ rg + n. - mp,) ^ 

^ rll - 971 -J- P/r, .I ^ rtl-97- ' ^ \ / J 


N-, 


i=i 


rO 


tv 


+ iuFj) r(2 - in + iupj) 


On the other hand, the solution of (2.19) for G reads as 


G(n — w) = 


27t' 


(1 + in — iw) + '0(1 — iv + iw)]. 


( 2 . 22 ) 


Then, we notice that in the large s limit Z 4 {v) ~ F{v) ~ —4nln(s/\/2) + 0(s°). This means that 
the leading behaviour of the nonlinear term in (2.17) is the same as that of the analogous term 
for the GKP vacuum: therefore, we can use results of [12] and approximate the non linear term in 
(2.17) as 


<•+00 


dw G{v — w) Li{w) = —2n In 2 + 0(l/s^). 


(2.23) 


Plugging (2.21) and (2.23) into (2.17), we eventually get 


Za{v) 


P (I + in) 

= —iL In —-h — 4n In s 


H 


Kb 

*^ln 

1=1 

Nf 

i 

1=1 


P (1 - iv) 

I + n - Ub,j 
I - n + Ubj 




h=l 


P(1 — in + iUh) 


Ng 

1=1 


P(1 + in — iuj) 
P(1 — in + in^) 


N-g 


+ i ^ In 




tUp. 


Ng 


P(| — in + iupj) 


1=1 

E , P(l + in-injfO 
In^l-^ 

j=i r(2-*n + *nFj) 


P(1 + in — iUj) 
P(1 — in + in^) 


+ (2.24) 


+ 0(l/s' 


Now, it is clear that imposing quantisation condition (2.14) on (2.24) provides a constraint between 
the rapidity Uh of a scalar and the rapidities of all the other excitations. As in all integrable models, 
this constraint has the general form (2.8): therefore we could be tempted to use (2.14) to dehne 
momenta and scattering factors of excitations, as well as the effective length of the chain. Such 
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procedure, however, will provide scattering factors i In S' which diverge as u* In u* when the rapidity 
M* of a generic excitation becomes very large. Fortunately, it happens that this problem can be 
avoided if we make use of the zero momentum condition, which is a selection rule to extract physical 
states out of the Beisert-Staudacher equations. To be specific, all physical states have to satisfy 
the condition = 1, where 


Since 


Ki 

P = i In 
i=i 

s + 2Ki, is even, we can also write 


'^4,j + I 

~ 2 


(2.25) 


Ki 


p = 


i In ^- — + 27rZ . 

^ 1 - 7/ . .. 


j = l 


Ua 


(2.26) 


This expression is regular for U 4 j = 0 and, therefore, it is more convenient for our calculations: 


P = 


K 4 i . H 

9 + U 4 j 

I > in 4 - - = I 


i Vln4- 

^ f - “4J 


XT?+ 7 


Ku 


i + 


i=i 
vrA'b — i 


f + CO 


dv 


In 


h=l 2 

h+V d 


+ Uh 


+ 


i=i 


^6, 


' -00 27r f - n dv 

Na 


[Ziiv) - 2X4(1;)] = vriFb - 1 In ^ ^ 

1 r (2 - lUh) 


h=l 


— t 


^ r (I - inf) ^ 


i=i 

Nf 


i=i 


T (I - inf) 


^ T (1 - inf) ^ 


in: 


0 


i=i 


i=i 


T (1 - inf) 


0{l/s^) . 


(2.27) 


As a technical remark, we notice that nonlinear terms give no contributions at the orders In s and 
(In s)°. 

Putting together (2.24) and (2.27) we obtain the equality 


Z 4 (n)-P = 


Kt, 

ij^ln 

i=i 

H 

+ i yf In 

h=l 

Ng 

+ i yf In 

i=i 

Np 

+ i In 
i=i 

T 

— 2i In — 


n - Ub,j + I 


r + 


V Ubj 2 

T (I - in) T (i + iuh) r(l + i 


in 


iwft) 


T (I + in) T (i - iuh) r(l - in + in^) 

T (I — in) T (I + inf) r(l + in — inf) 

T (I + in) T (I - inf) r(l - in + inf) 


+ 


+ (^ ^ ^) + 


(2.28) 


T (I - in) r(l + inf )r (I + in - inf) 

r(i + in)r(l-inf)r(i-in + inf) ^ ^ 


I-in) 


— 4nlns + 0(l/s^) 
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Therefore, we have gained the possibility to write the condition e d^4(«/i) p] _ (—1)^ ^ 
convenient alternative to (2.14): 


as a 


Vr(i + m,j)/ Y [^ + iuh)T [l-iuh')T{l-iuh + iuh') ^JiUh-Ub,j - ^ 

h'^h 

r (i + i{uh - u^j)) r (i - iuh) r (I + mj) r (i + i{uh - wf)) r (| - iuh) r (| + wf) 

r (1 - t{uh - nf)) r (1 + r (I - zu^) f}^ r (i - - ^f)) r Q + r (| - tu^) 

-j^ r(| + i{uh - UF,j))T{l + iuF,j)T{^ - iuh) r(| + i{uh - upj))T{l + iupj)T{l - iuh) 

r(i - i{uh - MFj))r(l - iuF,j)T{\ + iUh) r(i - i{uh - Ufj))T{ 1 - iuFj)T{l + iuh) 

(2.29) 


We take (2.29) as Bethe-Yang equations for scalars. In the spirit of (2.8) we make the following 
identihcations: 

• Length of the chain i? = 2 In s 

• Momentum of a scalar PQ^\uh) = 2uh 

The terms in (2.29) depending on two rapidities have the natural interpretation of scattering 
factors between scalars and other excitations. Using notations given in Appendix C, we write 


1 




T(| - iuhY ^ Uh - Ubj + f 
T{l + iUh)\ fJ^Uh-Ubj - ^ 


H 

JJ S^o'’\uh,Uh') ■ 

h' = l 
h'^h 


No 


No 


i=i 


(sg) 

0 


(uh, uf 


Ifs, 

i=i 


(«ff) 

0 


Np p 

{Uh, U^j) n Ufj) S^o^\uh, Upj 


(2.30) 


i=i 


i=i 


where denotes the scattering factors between a scalar and a generic excitation. We remark 
that i In behaves like In u* when the rapidity of an excitation becomes large. Eventually, the 
last term 


r(i - iuh) 
r(| + iuh) 


(2.31) 


has the form of the phase delay due to two purely transmitting defects. 


Finally, in view of generalisations to all loops we hnd convenient to identify the various pieces 
entering the function Z 4 ^{v) 


H 


Kb 


Z^iv) = 0 o(n,s/\/ 2 )+ 0 o(n,-s/\/ 2 )+ ^ 0 o(n,M/,)+i^ln 


h=l 

Np 


I + n - Ub,j 


^=1 2 - ^ + ^b,j 


+ 


Y1 “i) + “i) + Y1 ^0 “Bi) + Y1 ^0 '^F,j) - 2n In 2 , (2.32) 

j=i j=i j=i j=i 
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as solutions of integral equations 


0 o(n, u) = (j)o{v -u) + $o(n) - 

r (1 + zn — iu) r(l/2 — iv) 
r (1 — zn + iu) r(l/2 + iv) ’ 


^+oo 


dwtpoiv — w)Q'q{w, u) = 


= z In 


Fq{v,u) =xo(w-M|l) + $o(t') 


<• + 00 


dwipo{v — w)Fq (w, u) = 


= z In 


= Xo{v - u\2) + <|)o(n) 


r (1/2 + zn — iu) r(l/2 — iv) 
r (1/2 — zn + iu) r(l/2 + iv) ’ 


f + CO 


dw^po{v — w)Fq(w, u) = 


= z In 


r (1 + zn — iu) r(l/2 — iv) 
r (1 — zn + iu) r(l/2 + iv) ’ 


(2.33) 


(2.34) 


(2.35) 


where $ 0 , 4>o, Xo are defined in Appendix A. In (2.32) the large s limit has to be taken in the first 
two terms in the right hand side. This limit gives 


0 q(z;, s/\/2) + 0 o(z;, —s/\fi) —)■ —4z; In — 2 z In 

v 2 


r (i + iv) 
r (I - iv) 


(2.36) 


Fermions 

The equations for (large) fermions with rapidity XF,k = UF,k come from the (inverse of the) 
third of the Beisert-Staudacher equations. We have 


1 = 


Ka 

n 


'^F.k '^a,j “1“ 


1 Q I 

Upk -u^i + 


K 4 


n 


- Ml - ^ 


j=l “Bfc - ^a,j - 2 -j 2 

Ka I 7 Q I i 

UF,k — UaJ + 2 TT '^F,k ~ '^j F 2 




+ UAj — UF,k 


n 


n 


= 1 t + ~ '“4,j 

Kb TT ^ ~ UF,k + UbJ 


UF,k — Ua,j — 2 '^F,k ~ '^j ~ 2 




H 

n 


+ UF,k — Uh 


i + UF,k - UbJ f + M/i - UF,k 


exp 


27r 


V — UF,k d 


F- In — 2Li{v)) 


(2.37) 


I - n + UF,k dv 

since 7^4 = s + 2Kb is even. Then, we evaluate the integral (L 4 contributes with subleading 0(l/s^) 
terms) 


b + CX) 


dv 


In 


+ V - UF,k d 
— V + UF,k dv 


{Z 4 {v) — 2L4 {v)) = Lin 


r(l -iuF,k) 
r(l + iuF,k) 


+ 4:iuF,k In s — 


H 


j ^(f iUF^k) j r ( 

^ r(| - iuh + iuF,k) ^ 


+ tUj — tUF, 


iUF,k) 


Nb 


In 


T (I + zzzj — iuf, 


\k) 


j=i 

Ni, 


T [I - iu^j + iuF,k) ^ T - iu^j + iuF,k) 


h=l 
Np 

Vln 

^ l\l - lUF,j +iUF,k) ^ l\l - lupj + lUF,k) ^ i-UF,k + Ub,j 


r(l + iUF,j — iUF,k) 1 r(l + — ^'^F,k) 

r(l - iUF,j + iUF,k) 


Vln-)- 

^ rn - i 




Z T UF^k Mfoj' 


(2.38) 
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Now, in order to reproduce scattering factors already appearing in equations for scalars (2.29, 
2.30), we use the zero momentum condition: we multiply (2.37) with 1 = e®'^, with P given by 
(2.27). Using notations defined in Appendix C we write the final Bethe equations for fermionic 
excitations as 


1 = 


,iRP^^\uF,k) r(l - iUF,k) ] ^ 7^ UF,k - Ua,j +i/2 ^ ^(^Fs) 


n 


_r(l + iUF,k)\ UF,k - Ua,j - i/2 




Na 


Na 


Nf 




JJ Si,^^\uF,k, U^j) n So^’^\'^F,k, U^) n So^^\'^F,k, Ufj) JJ Si^^^\uF,k, Upj) , (2.39) 

j=l j=l j=l j=l 


where we introduced the length i? = 21 ns of the chain and the momentum PQ^\uF,k) = ‘2uF,k of 
a fermionic excitation. As for scalars, the term 

r(i - iuF,k) 

_r(i + iuF,k) 

stands for phase delay due to purely transmitting defects. 

Equations for large antifermions come from the (inverse of the) fifth of the Beisert-Staudacher 
equations and are obtained in a completely similar way as in the fermions case. The final result is: 


(2.40) 


1 ^ 


T{l-iuF,k) 




2 Kc 

n 


RF,k ~ RcJ + \ 


U. 


Na 


Na 


cj 2 


Na 


n ^o\^F,k, U^j) n ^0^\^F,k, Uj) n *^0 ^Fj) JJ Sl^^\uF,k, UF,j) , (2.41) 

j=l j=l j=l j=l 


where we introduced the length i? = 21 ns of the chain, the momentum Pq^\up i.) = 2up j^ of a 
antifermionic excitation and the ’defect’ term 


-iuF^k) y 

r(i + iup^k) _ 


(2.42) 


Gluons 

In the presence of a large number s of real type-4 roots, a gluon with rapidity u^. is described 
by a stack composed of a single type-2 root U 2 ^k = ul ^ two-string formed by two type-3 roots 
U 3 ,k = Uk i Rapidity is then constrained by the equation obtained by multiplying together 
the (inverse of the) second of the Beisert-Staudacher equations with U 2 ^k = uf. with the (inverse of 
the) third for u^^k = he. 


1 = 


Ng 

n 

j=i 

j¥=k 


ul-Uj + I 

9 9 ^ 

Ul — Uj — % 


+ 2 


1=1 


Ut 


UbJ 


Nf 

n 

1=1 


9 I i ^4: 

"^k ~ '^F,j I 


u 


Rf-, 


1=1 


uj - Ui,j - i 
Uk - UA,j + i ' 
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We concentrate on the last term in (2.43), which we rewrite as (remember that K 4 , is even) 


Ki 

n 


ul 


q 

J = 1 


M4j 


+ i 


Ki 

n 

i=i 


i-ul + U4,j 


Kb rSi 


n 

i=i 


H 

n 


i + ul-Uh 


h=l 


Ul + Uh 

The integral term equals 


exp 


^ + K- “4j 

dv ^ i — ul + v d 

In 


f-ul + Ubj ^-ul + Ub,j 
'^b,j 2 ~ 




2% i + ul — V dv 


{Zi{v) - 2Li{v)) 


dv ^ l-ul + V d or rmr (2 


— — In „ 

J 2ti i + Ul — V dv 


{Z 4 {v) — 2L4 {v)) = Aiullns — Lin 


" + O 


+ 


r (2 — iuh + iuD 
&''r(2 + m. 


No 


r (2 - iul + iul) 


tut 


^T.^'^Y{2 + iu^. 

j=i j 


tut 


Ks 

E 

i=i 


In 


r (I - O 

r (2 - + iul) 

r (2 + iui - iul) 


^(2 + 

r (I + iuF,j - iul) 


Nb 


^ il - ^^F,j + K) 


V In I ~ 


Putting the last two formulae together, we have 


Ki g 

Ul - U4,j - t AiuUns 

- P K 

ul - U4,j + i 
9 


n 

i=i 

Ns 

n 


n L 


r( 


i + O 


Kb 

n 


+ ■Qr(i 


iUh + iul) 


r (2 - iul + 


r (2 


n 


iul + 


J i=i 2 

'■ r(| 


+ ul 


Ub,j r(l + iuh - iul) 


iuF,j + iul) TT ^ (i 


n 


(2.44) 


+ 


(2.45) 


=1 r(2 + iul ~ 7=1 r(2 + iul ~ f=i ^ (I + 7=i ^ (I + ^'^fj - K) 


n 


iupj + iul) 


We now plug such expression in (2.43) and multiply the resulting expression by 1 = e*'^, where 
P is given by (2.27). We observe the exact cancelation of the term depending on type-5 isotopic 
roots and get the hnal set of equations, written in terms of scattering factors listed in Appendix 
C: 




r (i - O 

r (I + 


H 


Ylsl^"\ul,Uh) 


h=l 

Ng Ng Np 

n «?) n “f) If 

j=i j=i j=i 

j¥=k 


Ng 


U 


ky 


UF,j)Ylsif^\uluFj] 


(2.46) 


i=i 


where we introduced the length i? = 21 ns of the chain and the momentum Po^\ul) = 2ul of a 
gluon F+±. In this case, the effect of the two transmitting defects on gluons is 


r (I - ^ul) 
r (I + ^ul) 


(2.47) 
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In analogous fashion, we obtain the equation for the gluon held F+±: 




r (I - i4} 

r (I + m|) 


n 2 


H 




h=l 


Na 


iVs 


n “f) n “f) 

i=i 


j¥=k 


lf"\4 

t=i i=i 


n UF,j) n '^F,j) , (2.48) 


where again i? = 2 In s is the length of the chain and Pq yul.) = 2 m^ is the momentum of the gluon 
excitation P+±- 

Isotopic roots 

We remember (see (2.4, 2.5, 2.6) the dehnition of the three sets of isotopic roots, which do not 
carry momentum and energy, but take into account the sm(4) symmetry of the GKP vacuum. 

We have Ka roots Ua,j of type U 2 , Kc roots Ucj of type uq and Kf, stacks, = Ubj ± |, u^j = 
Ubj = Ubj with centers Ubj- 

The equations for the isotopic roots Ua and Uc come directly from the second and the sixth 
of the Beisert-Staudacher equations: we observe the cancelation of the contributions coming from 
gauge held stacks: 


1 = 


1 = 


Ka • Np i Kf) 2 

'^a,k ^ 1 r ^F,j 2 T T ^a,/c ^b,j 2 


n 


n 


n 


'^a.,k UPJ + 2 j=i '^a,k '^h,j + 2 

Kc . i i 

^c,k ^c,j ^ 1 r ^c^k ^_F,j 2 T r 2 


n 


n 


n 


(2.49) 


(2.50) 


'^c,k Uc,j i Uc,k Upj + 2 Uc,k Ub,j + 2 

For what concerns the isotopic roots Ub, we consider the product of the third Beisert-Staudacher 
equation for = Ub,k with the hfth for u^^k = Ub,k, the fourth for u^^k = Ub,k + V 2 and the fourth 
for u^^k = Ub,k — f/2. We arrive at the following equation. 


1 = 


» L i Kc ^ ^ i 

'^b,k ^ \ T r 2 T r 2 


n 


n 


^b,k ^ ^ "^b^k ^CL^j “ 1 “ 2 j=i "^b^k “ 1 “ 2 


K 4 

n 


Ki 


'bf^b,k ^4,j H” 2 T r ^b,k ^4,_7 2 


Si 


n 


II 7 11 i*^7lll 

; = 1 ^t',k 2 j=l ^b,k T ~2 j=l \^b,k ^b,j 

Kg Q Qy Kq Q ‘Xj N TT . K p 

Ub,k S 2 TT ^b,k S 2 TT Rb,k — ^F,j "F 1 TT Ub,k — ^F,j T I 


n 


^bj ^ 


(2,61) 


n 


3i 


n 


3i 


n 


n 


= 1 Ub,k -Uj- ^ 1 Ub,k - - Y 7=1 '^b,k - UF,j - I Ub,k - Upj 


We have 


K4 

n 


^ 6 ,/c ^4,j ~l~ 2 


K4 

n 


f + Ub,k ~ R-AJ 


Kb 

n 


'^b^k ^bj ^ 


. H 


n 


Rb,k T kLh 


Rb,k Raj 2 2—1 2 Rb,k T Raj Rb,k Rbj * 2 P Rb,k Rh 


j= 

j¥=k 


exp 


'■+°° dv ^ ^ + Uhk - V 


— In 
271 1 


Rb,k + R 


(Zjv) - 2Ljv)) 


[l+0(l/s")] 


(2.52) 
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Plugging (2.9) into the integral in the last term of (2.52), we hnd that 


exp 


dv ^ + Ub,k-v ' 

— hi 4 - - 2L^{v)) 


Ka 3i 


f —oo 2‘7r 2 ^6,/c "h 

3i 


* + Ubh 


L Kt, 

n 




^ ^b,k J 1 *h ^b,j 

’ ' .7 = 1 ^ 


,9 ^3 3i 


g Np . 


Ne 


n 


~2 + '^b.k ~ 'W4J -p-r ~2 ~ '^b.k + -p-r y “ '^b,k + -p-r i — Ufj^k + TT ^ “ '^b,k + 'Up 


2 ^b,k 2 '^b^k 2 "^b^k ^ ~l~ ^b,k ^F,j ^ “1“ ^F,j 


n 


3i 


n 


3i 


n 


n 


Putting all together, we eventually get the following equation, for the third isotopic root Ub- 


Ka 


i=n 


"^b^k ^CL,j 


Kc 

n 


^b,k ^c,j 


Ku 


n 


H 


^b,k '^b,j ^ 1 r ^b,k 2 


n 


'^b,k UaJ + 2 Ub,k Uc,j + 2 '^b,k UbJ * Ub,k Uh -\- 2 


(2.53) 




which does not depend on the roots associated to gluons. 


2.3 The general (all loops) case 


We now generalise all the results discussed in the one loop case to the most general all loops case. 
For the sake of clarity, the complete set of equations is summarised in Appendix E. 

As we did in the one loop case, we start from scalar excitations. 

Scalars 


Let us introduce the counting function for the type-4 roots 

Ki 


Z 4 {v) = iLln 

Kt, 


X n 


x^(v 


+ i y~~^in 

j¥=k 


X {v) — x'l ^ 


‘'4,i 


2x+(v)x,^^^ 2 


a;+(n) - x.. 1 - — 


a (n,M4,i) + (2.54) 


2x {v)x\ 


4,j 


E , X'^{v)—Xbi X'^{v)—XF'i X'^{v)—Xfa 

In ^ + * V In ^ i V In ^ + 

nri ■ — 'T" I I f ^ T • — X~ [V] ^ ^ 


i=i 

Nf 


Xb,j - X-[V) 


1 - 




i=i 

% 1 


xf,. 


X+('u) 


i=i 


^F,j 


x~{v 


I _ 1 _ 

j=l x~{v) j=l x~(v) 

Na 


+ 


X~^ (v) — X^,'^ X^, —X^{v) ^ X~^ (v) — X^,'^ X^: —x~^{v) 

+ f > In- FT— -^—I-f > In- 

^ x~ {v) — x^,'^ x~ {v) — x^, ^ 


i=i 


i=i 


x~{v') — x^~^ x~(v) — x^~ 


where a^(v,u) is the so-called dressing factor [47, 10], x^(v) = x(v ± i/2) and we are using the 
notations 




X 


± 

4J 


= X 



^b,j ^i^b,j) • 


(2.55) 
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The property = (—1)^“^ follows from the dehnition (2.54) and from the relation (2.13) 

between L and the number of the various excitations: the condition identihes 

the H internal holes, i.e. the scalar excitations. 

It is convenient to write (2.54) in terms of functions <I>, y, introduced in the Appendix A 


Z^{v) 


Ki Kf, 

L$(n) — 0 (n, Uj) + 2i In 

i=i j=i 


a;+(n) - Xbj 
Xb,j -x-{v) 


+ 


Ng N-g Np Np 

+ + ^x(n,Mf|l) + ^Xf(w,mfj) + ^Xf(w,Mfj) - (2.56) 

i=i i=i i=i i=i 

Nf Nf 

- '^XH{v,Ufj) -'^XHiv, U fj ) , 

i=i i=i 

where s of the type-4 roots involved in the sum are real, while 2Kb are part of the stack dehning 
the isotopic root Ub- We concentrate on the real type-4 roots and write the sum over them as an 
integral by means of (2.15), getting 


Kb + 

L^{v) -F 2i V] In - — -^ - y'[0(n, Ubj + i/2) -F 0(n, Ubj - i/2)] -F 

^ XbJ - X- [v) ^ 

j=i 'j j=i 

Ng Ng Np Np 

+ + + '^XFiv,UF,j)+ '^XFiv,upj) - (2.57) 

j=i j=i j=i j=i 

Nf Nf H 

- Xh{v, Ufj) - xh{v, Ufj) + Y '^h) + 
j=l j=l h=l 

r+oo 

+ (j){v,s/V2) + (j){v,-s/\/2) - / dw(p{v,w)[Z 4 {w) - 2 L 4 {w)] + 0{l/s‘^), 


Kb 


ZJv) = 


where (p is dehned in (A.4). Then, we can write 


<• + 00 


Z4{v) = F{v) + 2 / dwG{v,w)L4{w), 

J —oo 

where F{v) satishes the linear integral equation 


(2.58) 


F{v) = L^{v) + 2iY^^- — -^ ~ '^bj + i/2) + 4>{v, Ubj - i/2)] + 

rp, . - |7?1 


. Xb,j - X (v) 

J=1 J=1 


+ + Yx{x,U^j\l) + YxF{v,UF,j) + 5]^XF(t^,MFj) - 

j=l j=l j=l j=l 

Nf Nf- H 

- Y Xniv^Ufj) - Y Xh{v, Ufj) + Y '^h) + 

j=l j=l h=l 

r+oo 

+ (/(n, s/\/2) 4 -(/(f,-s/\/2) - / dw(p{v,w)F{w) + 0{l/s^), 


(2.59) 
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and 


<• + 00 


G{v,w) = Lp{v,w)— / dw'Lp{v,w')G{w',w) 


(2.60) 


We now work out the solution to (2.59). The part depending on the isotopic roots is written in an 
explicit form. For the remaining parts we remember that L = H + 2 + Ng + Ng + Np + Np: this 
allows to put L<F(n) together with the other functions in the right hand side of (2.59). Eventually, 
the solution to (2.59) is written in terms of solutions of linear integral equations. In specific, we 
have 


H 


Kb 


F{v) = Q'{v,s/V2) + Q'{v,-s/V2) + '^Q'{v,Uh) 


i/2 + V — Up. 


where 


Na 


h=l 

Np 


^ i/2-v + Ubj 


Nil 


+ 


F^{v, u^) + Y “f) + ^ ^F,j) + Y 


i=i 

Nf 


J =1 

Nr 


i=i 


j=i 


Y “Ft) + Y “hi) + V, 

i=i i=i 


Q'{v,u) = <p{v,u) + $(n) 


™+CXD 


dw(p{v, w)Q'{w, u ), 


b+cx) 


F^ {v,u) = xf{v,u) + ^{v) — / dwip{v,w)F^ {w,u), 


b+co 


F^{y,u) = —xh{v,u) — / dw(p{v,w)F\w,u), 


F^{v,u) = xiv,u\l) + ^{v) — I dw(p{v,w)F^{w,u). 


0 + 00 


G/ 


(2.61) 


(2.62) 

(2.63) 

(2.64) 

(2.65) 


We hrst analyse the s-depending terms. A tedious calculation® shows that in the large s limit 

0'(n, s/V2) + 0'(n, -s/ \/2) = In [-Av + Zbes{v)] - 2P{v) + 0(1/ s ^), (2.66) 

v 2 

where Zbes{v) = -Zbes{-v) and -FZbes{v) = gbes{v), ctbes{v) being the famous BES density 


[10] and P{v) is the solution of the integral equation® 

/■+°° dw 


P{v) = —<F(n) 


[</?(n, w) — (p{v, —w)]P{w). 


(2.67) 


Then, we pass to study the nonlinear term NL{v) = 2 dwG{v,w)L 4 {w). The same term was 
computed in [21], where only real type-4 roots were present. Here we can use the same results, since 

^Formula (2.66) clarifies the origin of the length 21ns and the two ’defects’: they are both due to the interaction 
with the two heavy large holes. 

®We write the kernel of equation (2.67) in an explicitly antisymmetric form in order to avoid one loop divergencies. 
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in the large s limit the leading behaviour of Z 4 {v) does not depend on the presence of excitations. 
We have 

In 9 

( 2 . 68 ) 


NL{v) = -2v In 2 + —Zbes{v) + 0(1/ s^). 


Putting everything together, we arrive at 


H 


Kb 


ZJv) = 


In s [-Av + Zbes{v)] - 2P{v) + ^ Uh)+i^\n- 


i/2 + V — Ubj 


N, 


N-o 


h=l 

Nf 


^ i/2 -V + Ubj 


+ 




^ F°(v ,«») + F°(k, «f) + Ufj) + 5^ 


upj 


(2.69) 


i=i 

Nf 


i=i 

Nf 


i=i 


i=i 


+ “hi) + , 

i=i i=i 

Now, as in the one loop case, before imposing the quantisation condition for holes, we introduce 
the momentum of the chain 


Ki 


P = i In 


X 


4,j 


,=1 ^ 4 ,, 


K 4 , H Kf, 

*^(“4j) = ^{Uh) + f ^ In 

j=l h=l j=l 


rv. + + 

■^b,j 


X 


b,j 


dv , , d , , , , 

:$(n) —[Z4(n) -2L4 (w)] , 


2e 


dv 


(2.70) 


where x^^ = x{ubj ± i). Terms containing L4 give no contributions at the orders Ins and (lns)°. 
Terms containing Ub produce only a term nKb- The dependence on excitations is worked out after 
inserting for Z 4 {v) expression (2.69). However, for our convenience we prefer to work directly on 
the expression Z^iy) — P: after some calculation (see Appendix D for details) we arrive at the 
expression 


Kb _ . i H 

Z^iv) - P = \n.s[-Av + Zbes{v)]-‘^P{v ) ^ Q('^)+ 


j=i “ 2 h=i 

Nf Nf Nf 

i In S^'^^\v,UF,j) + i In up j) + i In + (2.71) 

j=i j=i j=i 


Nf 


i=i 


Nf, 

V,Ufj) + E i\nS^ 

i=i 


N-a 




^Hn S^' 

i=i 




where we introduced the scalar-scalar phase 


0 (u,m) = e'{v,u)+P{u) =t\n[-S^^^'>{v,u)] (2.72) 

and the scattering factors S^^*\v,u*) between scalars and other excitations, which are listed in 
Appendix C. Imposing the quantisation condition e-d^4(«h)-P] _ ^Pg gnaj Bethe 
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equation for scalars 


1 -ff 

f n 

2 h'=l 

h'^h 

Nf ^f 

■ n n «Fj) «/,.) n®''”*"'- “/j) • < 2.73) 

j=i j=i j=i j=i 

where we introduced the length of the chain i? = 2 Ins, the momentum of a scalar excitation with 
rapidity u 

P^^\u) = 2u-^Zbes{u) (2.74) 

and the effect of the two purely transmitting defects 

2D^^\u) = 2p{u). (2.75) 


i=i i=i 


Kb 


1 = 


_ {ub)+2iDM (ub) JJ" ^^4 


Ufi ^b,j 


Important properties of the scalar-scalar phase (2.72), which can be proven using equations (2.62, 
2.67) are 

Q{u,v) = —Q{—u,—v), Q{u,v) =—Q{v,u). (2.76) 

Eventually, we remember an efficient way proposed in [23] to compute the scalar-scalar phase. We 
found that 

Q{u,v) = M{u,v) — M{v,u), (2.77) 

where M{u,v) = + Z{u;v) and Zd)(-u), Z{u;v) are univocally defined by the conditions 

—Z^^\u)=a^^\u), -^Z(m; u) = (t(m; u), Z^^\u) =—Z^^\—u), Z{u]v) =—Z{—u]v) ^ 
du du 

(2.78) 

with the functions a‘d'>{u)^a{u]v) solutions of equations (B.22, B.23), respectively. 

This procedure provides an alternative (with respect to solving equation (2.67)) way to deter¬ 
mine the function P(m), once Zbes{u), Z^^\u) and Z{u;v) are known. Indeed using (2.66) we 
have 


2p{v) 


lim 

5^ + CXD 


lim 

5^ + CXD 


In ^ [-4u -h Zbes{v)] - 0(u, s/V2) - Q(v, -s/V2) = 
v2 

In ^ [-4u + Zbes{v)] - 2Z«(u) - Z{v, s/V2) - Z{v, -s/V2) 
v2 


(2.79) 


A hnal alternative to compute P{u) is to look at equation (2.69) when no excitations nor isotopic 
roots are present. Then we see that —2P{u) represent the contribution (^(Ins'’) to the twist two 
counting function of the pure sl{2) sector. This function has been analysed in [13, 14] (in notations 
of the second of [14] it is connected to the function 5'®^*'’“). 

Fermions 
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The equations for large fermions come from the (inverse) of the third of the Beisert-Staudacher 
equations. We have 


1 = 


-r-x UF,k - Ugj + i/2 J-J UF,k -Uj+i/2 XF,k ^ 

7=1 “Bfc - “aj - V2 UF,k - - V2 ^ XF^k - 

TT ~ ^gj + ^/2 t4 ^F,/c U^+l/l y-j- X^ - Xp^k yr Xp^k ~ 

UF,k - Ua,j - */2 UF,k -U]-%I2 XF^k “ X~^J x+ - XF,k 


i=i 


exp 


i=i 

XF{v,UF,k)-^(zA{v) - 2 La{v)^ 


' —OO 


27r^ 


(2.80) 


since = s + 2Kb is even, where x^ = x (uh ± f)- As in the one loop case, we multiply such 
expression by 1 = e*^. Then we use expression (2.69) for Za{v) and remember that the term 
containing L^iv) gives subleading 0{l/s^) contributions. We get 


1 = ^F,k Ua,+l/2 " . 

UF,k - Ua,j - 112 


(2.81) 


J=1 


Np Np _ Nf Njr 

■ JJ S^^^\uF,k, UF,j) S^^^\uF,k, Upj) JJ S^^^\uF,k,Ufj) JJ S^^^\uF,k, Uf,j) ' 
i=i i=i i=i i=i 


i=i 

Ng Ng 

i=i i=i 




where i? = 2 In s is the length of the chain and 
P^^\u) = 

2D^^\u) = 


^ [xf{v,u) + xf{-v,u)] 


' —CX) 

/•H-oo 


(XEEsix) 


^ [XFiv, u) + Xp(-u, m)] -pPix) 


(2.82) 

(2.83) 


are the momentum of a fermion and the effect on it of the two defects. 

The equations for large antifermions come from the (inverse of the) hfth of the Beisert- 
Staudacher equations. Their derivation is analogous to the fermionic case: 


1 = “F,t Ucj + a/2 TT g(F.)(^ , 

fJi Upx - Uc,j - 112. 

Np ^ Np ^ ^ Nf ^ Nf 

■ JJ {Up^k, UF,j) {up^k, Ufj) n (^F,k, Uf,j) n ^'^F,k, Ufj] 

j=l j=l j=l j=l 


(2.84) 


j=i i=i 


u 


‘j! ■ 
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Equations for small fermions are obtained starting from the (inverse of the) hrst of the Beisert- 
Staudacher equations. We have 


1 



^f,k “1“ f/2 

Uf,k - U2,j -il2 


Ki I- 


n 

i=i 


x+ 

4,j 


1 - 




TT Uf^k — Ugj + i/2 -r-j- Uf^k 

Uf,k - Ua,j - I /2 Uf^k - - i/2 


n 


i=i 


1 - 


^f,k 

b,3 


^f,k 

x7~ 

0,3 


H \ 

n g;. exp 



dv , , d 




(2.85) 


In contrast with the fermionic case, we do not multiply this equality by 1 = e*'^. We use expression 
(2.69) for Z^^iv). Working out the various terms we get 


1 


where 


,iiJP(/)K.fc)+2iZ)(/)(n^,,) - '^a,j + t/2 Uh) ' 


( 2 . 86 ) 


Nf 


Np 


Nf 


Nf 


i=i 


i=i 


n uF,j) n upj) n Uf,j) n s^^^~\uf,k, ufj) 

i=i i=i 

Ng Ng 

n (“hfc’ “f) n Uj) > 


t=i 


j=i 


P^f\u) 

2D^f\u) 



[xh{v,u) + xh{-v,u)] 
[XHiv,u) + Xh{-v,u)] 


1 

d_ 

dv 


O'BEsiv) 


P{v). 


(2.87) 

( 2 . 88 ) 


In a completely analogous way we work on the (inverse of the) seventh of the Beisert-Staudacher 
equations, which gives the quantisation condition for small antifermions: 


^ ^^RP(fHnfP+2iD(fHnF,) fj ^cj + ^/2 ^ ) . 


1 = e 


'll f U lie j ^ /2 
7 = 1 pk c,j / 


(2.89) 


Nf 


Nf 


Nf 


n ufj) jj upj) n Ufj) jj ufj) 

j=l j=l j=l j=l 


Na 


Na 


n s<*' («/.», u?) n («/>.«?) 


i=i 


i=i 


Gluons 

As in the one loop case, we multiply (the inverse of) the second of the Beisert-Staudacher 
equations for U 2 ,k = u^. with (the inverse of) the third for i/2 and (the inverse of) the 
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third for — i/2. We then get the following equations for the center of the gluonic string 


q 


N, 


1 = 


n 


g q 

ui 


“ “i + ^ TT - ^tj) fr “fc “ + ^2 '^f,3 + V2 


=1 “fc 


n 


j=i 


K - “i - * fJi fJi 4 - ^FJ -i/2 ul - Ufj -i/2 


n 


n 


j¥=k 

TT 4 ~ '^b,j + i/2^ 
ul - Ubj - i/2 ■ 

Making explicit the type-4 roots, we arrive at 


(2.90) 


N, 


9 ^.9 


4~4 ^4~ + i/2 Y\ 4~ '^f,j + V2 4~ '^b,j + i/2 


N 


1 = 


n 

i=i 


H 


4-4 


n 


n 


^ n,9 


Ki 


- * fJl 4 - UF,j - i/2 ul - Ufj - i/2 A A ul - Ubj - i/2 


n 


3=^ " 3^ 

^^ 94 - H—hW^H—h 


(xi - x-jixr - ft, - 4jmj - <~) 


(^9+ _ ™-l-W™+ _ ™9-\ 111 ) / O- _- \(rr9+ 

b^iK-^k -^h/K-^h -^k ) j=i \-^k -^bj Jy-^k 


X 


bj 


exp 


‘L |(xf’.«l|l)4(z4W-2L4W) 


(2.91) 


Following what we did for (large) fermions, we multiply such expression by 1 = e*^. Then we use 
expression (2.69) for Z 4 {v). We observe the exact cancelation of terms involving the isotopic root 
Ub and eventually for the field we obtain the equations 


Na 


H 


I ^ giflp(9)K)+2iD(9)(n9) jj (2.92) 




i=i 


h=l 


Nf Np _ Nf Nj 

n Ufj) J]^ S^>^^\ul, Ufj) JJ Ufj) JJ S^sD^ul, Uf-j) , 

j=l j=l j=l j=l 


where 


P^^\u) = - 

2D^s\u) = - 


/■+°° dv 

I ^ [x(^^,m|1)+ x(-p,m|1)] 

/■+°° dv d ~ 

[x{v,u\i) + x{-v, m|i)] ■ 


(FBEsiv) 
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dv 


(2.93) 

(2.94) 


The procedure for the field F+x is completely analogous, hence we give only the final equations 

Ng Ng H 


i=i 


h=i 

Nf 


j=l,jj^k 

Np F[p Nf J V f 

n ^Fj) n“p,j) n ■«'’'’(“?> “«) n ^^>(4, ■ ■ 

j=i j=i j=i j=i 
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Isotopic roots 

The equations for the isotopic roots Ua and Uc come directly from the second and the sixth 
of the Beisert-Staudacher equations and their derivation is completely analogous to the one loop 
case: the only difference is that in the general all loops case also small fermions are present. 


1 


1 


Ka . ■ Np ■ /n ■ /o • /o 

1 T ^a,k ^a.,j T ^ I T '^a,k ^F,j I j2, -i r Ua^k % j2 -i r Ua^k ^b,j 112 

^ j—\ ^a,,k ^F,j T i/2 Ua^k '^f,j T i/2 '^a,k ^b,j T i/2 

iCc . ■ • /o ■ ic Kb . 

TT Uc,k — UcJ + I J-T Uc,k ~ Upj — ij 2 j-j Mc,fc ” Ufj — ij 2 j-r Uc^k — UbJ — ij 2 

i j—i ^c,k ^F,j T i/2 ^c,fc ^f,j i/2 '^c,k ^b,j T i/2 


(2.96) 

(2.97) 


Then, we consider the product of the third equation for = Ub,k with the hfth for 

the fourth for = Ub,k + i/2 and the fourth for u^^k = Ub,k — i/2. We arrive at the following 

equation 


Ka 


• /r\ • /r\ -^4 • /r\ 9 Q ' /O 

^b,k ^/2 1 r tlb,k ^/2 1 r Ub^k T ^/2 -i r "blb^k % f 2 


Ka 


No 


1 = n 


n 


^b,k ^/2 ^6,/c H“ ^/2 z/2 “i“ '2'/2 


n 


n 


n 


Ub,k - - i/2 / x” 


^<"4 1 - 


\^b,h 


n 




+ +-- — 


2xrra; 


6,fc 4j 2 




++ 

b,fc 


^4,i 1 


^ {Ub,k T i/2, 'U4j)(J {Ub^k i/2, U4.j') 


‘^^b,k ^4j 


n 


.,++ 


^b,k ^b,j 


Nf H—h H—h 

■^b,k -^F,] -p-r J^b,k ■^F,j 


1 - ^ 1 - ^ 


n 


n 


=1 \%,fc ^b,j j Xpj Xpj 


iiTv^n 


^b,k 


Kg ++ a+ ++ a- Kg ++ 


n 


■'b,k 


■'b,k 


g-\- - 

rtp /y'-' /y rtp' 

Jjh h -F ^ ^ Jj 


„9- 


7 = 1 ■^b,k 


n 


-'b^k 


g-\- -\—|- 

^b,k 


J = 

g- 


1 _ TL 
^77 i=i 


1 - 


^b,k 


j = l ^b,k 


g~\~ — 

rtp /y'-' /y rtp' 

i, O/^- ^ O/ 


where 

We have 


L = H + 2 + Np + Np + N, + Ns. 


(2.98) 


(2.99) 


Ka 

n 


^b,k ^4,_7 “i“ ^/2 


n 


. n 


^ ^b,k '^/2 


'^b,k ^4,j z/2 y^b^k "^bj ^ h—1 '^/2 

j¥=k 

dv Uh^k — V + i/2 d 


( 1 + 0 ( 1 /*^)) 


exp 


In 

27r Ub,k — V — i/2dv 


(Ziiv) - 2L4 {v)) 


( 2 . 100 ) 


where for Z 4 (n) it is convenient to use form (2.54). It is remarkable that, plugging (2.54) into the 
integral in the last term of (2.100), we hnd that 


exp 


/■+°° dv Ub,k - n + i/2 d ' 

/ — - 1 — l~ZAv) 

l_^ 2 f Ub^k — v — i/2dv 


( 2 . 101 ) 
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produces massive cancelations in (2.98). On the other hand, the nonlinear term containing L'^{v) 
gives a negligible 0(l/s^) contribution. Eventually, for the third isotopic root Uf, we obtain the 
same equation as in the one loop: 

^ 1 r ^o,,j 1!2, 1 r U\) ]^ ^c,j %!2 1 r ^b,j “I” ^ 1 T ^b,k 1 j2 102) 

• 'bl‘b,k ^a,j “b */2 'lt‘b,k ^c,j "h */2 ^b,j * '^b,k '^h "b i/2 

3 Conserved observables 

Momentum was already obtained in previous sections: therefore, we concentrate on higher charges 
Qr and in particular on anomalous dimensions 7 = Q 2 - Let us introduce the function 


qr{u) = 


ig 


x+{u) 


r—1 


X~[U 


r—1 


r > 2 


whose Fourier transform reads 


V — 1 

2 / \/2 \ -\M Jr-i{y/2gk) 


qr{k) = 2TTig | — 1 e 


k 


(3.1) 


(3.2) 


The r-th charge of an excited state over the GKP vacuum enjoys the expression 

K4 Ki H 

Qr ^ ^ qr{Ui,j) ^ ^ \qr i,^b,j “b i/2'^ T Qri.'^bJ ^/2)] ^ ^ 

j=l j=l h=l 

/‘+°° dv d ( \ 

~ / ~ 2L4(ti)] + O 


(3,3) 


^°This cancelation was already noticed and proven by Basso in Appendix C.2 of [20]. 













where for Z 4 (v) we use expression (2.69). Doing this, we observe the exact cancelation of the 
dependence on the isotopic root Ub and we are left with the formula 


Qr 


^gr(Mft)+lns j ^qr{v)-^[4v - Zbes{v)] + j ^qr{v)^2p{v) - 


h=l 


, , 1^9 Ng Nf 

^^r{v)—^^Q'{v,Uh) + '^Fg{v,U^j) + '^Fg{v,U^j) + '^ Ff{v , Upj) + 
h=l j=l j=l j=l 


Ng 


Nf 


Nf 


Ff{v,uj j) + 2 / dwG{v, w)L 4 {w) — 2L4{v) 

j=i j=i j=i d 

, f dv , . d , „ , ,, f dv , . d , 


H 


No 


No 


Z. <5t’ (“O + E ‘3‘"’ ('“?) + E ('“f)+E M+ E cf ’('“F,y) + 

h=l j=l j=l j=l j=l 

N, Nf 

E ‘3t’(«/,,) + E QV\N.i) + 0(l/i>"). (3.4) 


i=i 


i=i 


The first two terms in the right hand side of (3.4) are contributions from the GKP background. 
The remaining terms in (3.4) are the contributions that any single particle brings to the overall 
value of the r-th charge. 

• For scalars we have 

Qr\'P = -qAu) - j ^^r{v)^Q'{v,u). (3.5) 

Restricting to r even, we use relation (2.21) of [23] to write 

/ +00 fjU 

—qr{k)[a^^\k) + a{k]u)], (3.6) 

where the functions a'd'>{k), a{k;u) satisfy equations (2.19), (2.20) of [23], respectively. It is 
convenient to introduce the functions, dehned for k > 0, 


S^^\k) 

S{k] u) 


sinh I 
Tik 


a^^\k) + 


sinh I 
irk 


d{k] u) — 


sinh I ' 
27re-*^ 

1 — e~^ 



(cos ku — 1) 


and to expand them in Neumann series 


^W(fc) = S{t,u) = 


P=1 


P=1 


(3.7) 

(3.8) 


(3.9) 
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We use defining equations (3.10) of [48] for and (3.15) of [21] for Sp{u) to simplify (3.6) as 
follows 


r—1 


r — 1 \ zg 


(3.10) 


We could not find a formula analogous to (3.10) in the case r odd. 

When r = 2, this simple expression can be connected with the first of (4.6) of [20]. We indeed 
remember formula (4.36) of [49] and that 


+°° / 1 \n^,2n 


(3,11) 


where for s[^'^ we use formula (4.35) of [49] to connect with the solution of the BES equation. 
Operating in this way we get, after some algebra, 


e<‘>(ti) = 7'''(u) = 


r’+co 


dt 


e 2 — cos tu a , , cos tu — e '2 « , /- , 

^- ^Y_{V2gt) + ^— ^Y+i'^gt) 


62 — e 2 


62—6 2 


(3.12) 


i.e. the hrst of (4.6) of [20]. Functions 7 ± are defined in (B.19, B.20). 

It can be of interest to express in the 0(6) limit [50]. We introduce 

, , 287r4 1 _Z9. r ^ /l 

^{ 9 ) = —1 + 0 - 

rd)" [ \g 


(3.13) 


In the 0(6) limit \/2gS^'* = m{g) — 1 and \/2gS[{u) = m{g) (cosh — l). Substituting in (3.10) 
we get 

(3.14) 


TT 


7 ^®^(m) = m{g) cosh -^u — 1, 
and for the complete anomalous dimension in presence only of scalar excitations 

H 


liis/( 5 f) + fsi{g) + ^ {rn{g) cosh^Uh - l) 


h=i 


For gluons we have 


q ( 5 )( m ) = - I ^q^(v)^F^{v,u). 


Using equations (2.62, 2.63), we arrive at the formula 




q7(m) = 

When r = 2 

we have 


Qi^\u) 

= 7*'®^(m) = 

dk 

y_oo 47r2 g 

2% 

CX) / 

6-^ V 1 

( ^ 

ik 

n=l \ 

\^\/2ix{u + 


^ [x{v,u\l) + <f>(n)] ^Qr\'v) • 


Y+{V2gk) - sgn{k)Y-{V2gk)) 
9 


27r 

ik 

271 




(3.15) 


(3,16) 


(3.17) 


(3.18) 


/n- , i\ I I 7n(V^a>c) - — J(,(V^glc)e 1 
V2zx(u - i) I / IK 
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By means of the relation (3.40) in [20], and making use of the identity 
Jo k " \ oPiix(u T i) 


(3.19) 


the expression above becomes 




fc 1-e-^ ^ ^ 


dk 'y'^{y/2gk) 


k 


[cos kue 2 — 1]. (3.20) 


Large fermions and antifermions with rapidity u carry an amount of r-charge equal to 

Qf\u) = - I ^^r{v)-^FF{v,u). (3.21) 


Using equations (2.62, 2.63), we arrive at the formula 

dv 


Qi^\u) = - 

For small fermions we have 


d 


2^ [Xf{v,u) + <F(u)] ^Qr\v) 


Ql^\u)= I ^XH{v,u)^Q[f\v) 


(3.22) 


(3.23) 


/ p\ 

When r = 2, the very same reasonings outlined above apply to large fermions, so that (Q^ yu) = 
dk H ^ Y+{V2gk) [-cos ku - 1]. (3.24) 




— 1 


Analogously for r = 2 and small fermions 

Q2\u) = = 


= _1 f ^ 'y^(y/2gk) cos ku . 


(3.25) 


3.1 One loop 

All the previous expression are explicitly computed at one loop, upon introducing the following 
notation for the derivatives of the digamma function 

( 3 . 26 ) 

= ij{z) . ( 3 . 27 ) 

• Scalars 

(?'■'(«) = y<’'-">(l/2 - tu) + (-l)'V>'’'-"*(l/2 + tu) - 7’’-"'(l)(l + (-1)')] ; (3.28) 

• Fermions and antifermions 

Qf’(«) = -7)^^ - «■) + (-l)'V-<'-"'(l + «■) - + (-1)')] ; (3.29) 

• Gluons 

Q7(«) = P'’-">(3/2 - m) + (-l)’'^(’-^)(3/2 + in) - lf<’-=>(l)(l + (-I)')] ■ (3.30) 
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4 Strong coupling regimes of 2D scattering factors 


In this section we want to give a detailed analysis of the different strong coupling limits of the 2D 
scattering factors S**'{u,v) of sub-section 2.3. In fact there are different ways of performing the 
g —?■ -|-oo limit as these give rise to different results or regimes, so paralleling what happens to the 
energy/momentum dispersion relations [20]. 

First, we shall discuss the regime, relevant only for scalars (as the other excitations decouples to¬ 
wards very high energy), where we keep their rapidities fixed, namely the so-called non-perturbative 
regime. In this case integrations inside the expressions for the various scattering factors receive 
the leading contributions from the region where the integration variables are fixed (while sending 
g -|-oo). This regime is dominated by scalars which are the only ones to have a non-trivial 
(hnite) 5'-factor, whilst the other S'-factors involving other excitations reduce to one. Here we 
find out the (usual) 0(6) non-linear sigma model scattering theory as low energy string theory 
[50, 48, 51] Alternatively, we can hrst rescale the external rapidities u = \/2gu, v = y/2gv and 
then send g —?■ -|-oo. If the rescaled variables, u and v, have modulus smaller than one we are 
(with the exception of scalars, see discussion below) in the perturbative string regime (where the 
irrelevant and relevant perturbations of the 0(6) non-linear sigma prevail on it putting at zero its 
mass); while if their modulus is greater than one we are in the so-called giant hole (semiclassical 
soliton) regime. In both cases, in order to have the maximum contribution to the integrals, after 
rescaling external rapidities, we have to perform the same rescaling of the integration variables 
Ui = y/2gui and eventually take the limit g —)■ -|-cxo. 


4.1 Scalars 


Scalars in the non-perturbative regime 

We report the strong coupling limit of the scalar-scalar scattering factor in the non-perturbative 
regime, i.e. g —)■ -|-oo, with u,v hxed (details on the calculation can be found in [23], see also [21] 
and [22]): 


g —)■ -|-oo 


Q{u,v) Qnp{u-v) 


r(l-^V)r(i + ^V) ( n{u-v) 

r(i + z^)r(i-z^) V 2 


(4.1) 


which depends only on the difference of the rapidities and coincide with the pre-factor of the 
A-matrix, as derived in [53], of the 0(6) non-linear sigma model upon the identification (of the 
hyperbolic rapidities) 9 = nu/2 and 6' = ■nv/2. This definitely supports the proposal of the latter 
model by [50] as that describing the string at low energy (see also subsequent studies [48, 51]). 


Scalars in the perturbative regime 

next approximation it would be perturbed by irrelevant fields as suggested by the expansion of the energy 
in inverse powers of the size (i?) [21] (c/. also the dispersion relation in [20] and the effective field theory of [52]). 
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As in [20] the perturbative regime for scalars is recovered by introducing a new rapidity z as 


M = — In — 

TT m{g) 


Z 2 7Tl ( 0 ) 

, for M > 0 (right mover); m = — In-, for u < 0 (left mover) , (4.2) 

71 Z 


m{g) being the non-perturbative mass (3.13). The rapidity is kept hxed in the region m{g) < z < 
1 as g ^ + 00 . If m{g) < z < 1 formula (4.1) is valid, therefore the function 0 in the perturbative 
regime is obtained by plugging (4.2) into (4.1). 

Scalars in the scaling regimes 

We rescale the rapidities u = y/2gu, v = y/2gv and then send g ^ +oo, with u, v fixed and 
|m| > 1, |h| > 1. This is the so-called giant hole regime. Details on the calculation can be found in 
[23]. We give the final result for the double derivative, which will be useful for next computations 


-^-^e{V2gu, V2gv) = V2g 
du dv 


/s+iDA -I- fH+ii 

W-l/ V»;+l/ W+l/ Vii-l/ 

u — V 


1/4 


+ 0(9°). 


(4,3) 


This result (4.3) agrees with corresponding formula coming from using the scattering phase (2.34) 
of [54]. 

Another possibility is to define rescaled rapidities (with a bar) u = \/2gu, v = y/2gv and 
then send g —?■ -|-cxd, with u, v hxed and Ihl < 1, |h| < 1. Although for the other particles this 
second possibility gives rise to the perturbative string regime (giving for understood an obvious 
modihcation u —)■ x{u) for the (small) fermion, cf. below), it does not in the case of scalars as given 
in [20], because of their non-perturbative, dynamically generated mass. Yet, we need to consider 
the scalar 0 in this regime at least to access the other S'-matrix elements (depending on it). In 
fact, we may write the limiting value 


dudv 


0(m, v) = 27t—6{u — u) -|- 
du 


1 1 

y/2g dv"^ V — u 


+ 0{l/g^). 


(4.4) 


Importantly, formula (4.4) is valid also in the domains Ihl < 1, |h| > 1 and lul > 1, |h| < 1. We 
will make frequent use of (4.4) in this Section. 

From (4.4) we can infer 

Q(u,v) = —7rsgn(u — v) -^ ^ ^ -|- 0(1/g'^) ^ (4.5) 

V^.gv - u 


4.2 Gluons 


Gluons in the pertnrbative regime 

We want to study the gluon-gluon scattering factor (C.13) in the limit g -|-oo, with u = u\/2g, 
V = vV2g, u, V hxed and < 1, < 1 (perturbative regime). 

We have 

i In u)) = Xi + J 2 + X 3 , (4.6) 
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where 


Xi = x{u,v\l,l) =—2aictaii\/2g(v — u) = 

., 3 -, 


= — 7 rsgn(u — m) + 


OiMg^) . 


g{v - u) 

Passing to study X 2 and X 3 , we first remark that in the perturbative regime 


(4.7) 


x^{u) 


g 1 

— H- , 

\/2 4Vl - 


{u =F — u^) + 0{1/g). 


(4.8) 


Since we have to work out M|l) + <h(ta), in addition to (4.8) we need also to know the behaviour 
of x^{w) when w = w\/2g and g —?■ +cx). When |tD| < 1 we can use (4.8). On the other hand, for 
Itcl > 1 we have 


x'^{w) = y/2gx{w) ± 



+ 0{l/g), 


x(w) 



(4.9) 


Using results (4.8, 4.9), we arrive at the relations, valid for w = w\/2g and g +cxo: 

x(m;,m| 1) + $(«;) = ^ -— + 0(1/^^), when |w|>l, (4.10) 

V2.g{u - w) 

x(t(;, m| 1) + $(tc) = 0{l/g), when |t(;| < 1. (4-11) 

Therefore, we have 

/•+°° g 

X 2 = -y —[^(yj^u\l) + ^w)]—[x{w,v\l) + ^w)] = 0{l/g^). (4.12) 

For what concerns the last term X 3 in the rhs of (C.13), we hnd convenient to perform the change 
of variables w = '/2gw, z = V2gz: 


dw 


X 3 = 


dz 


7^/ 7^ [x(w^,m| 1) + $(w)] 


27r 


27r 


■^^Q{V2gw, V2gz) 
dw dz 


[x(z,n|l) + <l>(^)]. (4.13) 


Now, from formulae (4.3, 4.4), we deduce that the leading behaviour of the double derivative of the 
scalar-scalar phase is realised in the giant hole regime |ti}| > 1, \z\ > 1. Therefore, we can write 


^ f dw 

X3= / — 


dz 


d d 


J|tD|>i 27r J\ 2 \>i 2 vr ^/2gw — y/2gu y/2gz — \/2gv dw dz 
Plugging (4.3) into (4.14) and performing the integrations we arrive at 


Q{V2gw, V2gz). 


(4.14) 


X 3 = 


2y/2g{u — v) 


2 - 


1 + u 
1 — u 


1/4 


1 — V 
1 + v 


1/4 


1 — u 
I + u 


1/4 


1 + 

l-v) 


1 / 4 ' 


(4.15) 
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Now, summing up (4.7, 4.12, 4.15) we obtain the final result for the gluon-gluon scattering phase 
at the order 0{l/g)-. 


^(99) = exp 


i 

V^g{u 


v) 



n + u 

\l-u 


V4 

[l + vj 


1 n-uY^ + 
2Vi + h; u-h; 


The expression above agrees with the correspondent result of Basso, Sever, Vieira [34]. 


i/4y 

(4.16) 


Gluons in the giant hole regime 

We now want to compute the gluon-gluon scattering factor (C.13) in the limit g —)■ -|-oo, with 
u = u\/2g, V = v\/2g, u, v fixed and > 1 , > 1 (giant hole regime). 

As a preliminary calculation we consider the quantity x(u,m| 1) -|- <h(u). For its expression we 
refer to (A.7). Computing the scaling limit u = u\/2g, v = v\/2g, g —)■ -|-oo, h, v fixed and > 1, 
> 1 , we find that 


x(u, m| 1) + *h(u) = 7 rsgn(u — u) — vrsgnu -|- 0(1/g) . 


(4.17) 


Since x(u,m|1) -|- <F(u) is at most 0{g^), properties (4.3, 4.4) imply that the following part of the 
integral X 3 , 


7> 

-^3 


dw 


dz 


■^-^Q{V2gw, V2gz) 
dw dz 


[x(z,u|l) + 4>(z)], (4.18) 


gives actually the dominant contribution (proportional to g) to i In (u, v)). The integrations 

in (4.18) are easily performed and the final result^^ is 


iln (—u)) = iln (^—S^^^\u,v)) = Q{u,v) + 0{g^) = iln (^—S^^^\u,v)) + 0{g^). (4.19) 


Gluons in the non-scaling regime 

In this regime we send g -|-cxo keeping the excitations rapidities fixed. For what concerns 
gluons, if we send g —?■ -|-oo, with gluons and scalar rapidities, u, v respectively, fixed, we get that 


X{v,u\l) + $(u) = 0(1//). 


(4.20) 


In order to get (4.20), relation 




q —)■ -l-cxo , u fixed =7 x^iu) = ±—^ 

/2 


u± 




1 )' 

4^29 


+ 0(1/9= 


(4.21) 


is useful. Result (4.20) means that in this regime the gluon-gluon scattering phase S^^^\u,v) 
reduces to 


4^In order to get (4.19) we use the properties 0(u, ±\/2(/) = Q{±^g, v) = 0 which are proven using expressions 
given in [54]. 
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4.3 Fermions 


Fermions in the perturbative regime 

We want to find the strong coupling limit of the fermion-fermion scattering factor in the per¬ 
turbative regime. As we will show in a moment, this regime fits in the small fermion case. We 
start from 


5 '(//)(.y^ n) = exp 1 
where 


dw . . d , , . 

-Xh{w,u)—xh{w,v) -1 


2ti' 


dw' 


1 _ 3:/(u) 

( \ X+{w] /X U 

XH{w,u) = -t\n - Xf{u) = - 


(4.22) 

(4.23) 


X (w) 






In the perturbative regime the fermion rapidity scales as Xf{u) = \/2gxf{u), u = \/2gu, with 


Xf{u) 




m| > 1 , 


k/(w)l < 2 • 


(4.24) 


It is then clear that we are in the small fermion case. 

For what concerns scalar rapidity, we make the rescaling w = \/2gw and we develop at strong 
coupling. We have to distinguish two cases. 

• If |t(}| > 1, then 


x^{w) = \/2gx(w) ± 


1 + \/l — 


1 - ^ 


+ 0{l/g), x{w) = 


w 


1+f 1 




In this case we have 


^ _ Xf{u) 


( \ -1 x+{w) ^ ^fW) 1 

XH{w,u) = -i\n- -^ = 


X {w) 


V2g ^Ji_4_Xf{u)-x{w) 


If |ta| < 1, then 


In this second case 




x^{w) = —j=[w ± i\/l — w'^] -f 0 (( 7 °), 


/ N •. 1 - 2xf{w — i\/l - w^) , , 

Xh{w, u) = -t In --+ 0{l/g). 

1 — 2xf[w + 1 — w ) 

In general we write —i In = Ji -|- J 2 , where 

dw d \ T> T< 

^ = A +h > 


dw' 


(4.25) 


+ 0(1//). (4.26) 


(4.27) 


(4.28) 


(4.29) 
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where 


and 


where 


ir = 


= 


'|ui|>l 


|ui|<l 


^Xh{w,u)-^Xh{w,v) ~ 0{l/g ^), 
^Xh{w,u)^Xh{w,v), 


dw' 


h = - I ^^_XH{w,u)-^e{w,z)xH{z,v) = 


271 271' 


dwdz 


I2 = - 


jrest 

^2 


'hl,hl>i 


r]'7lj rj r rj^ 

-_Xh{w, u)——Q{w, z)xh{z, v) , 


271 27i' 


dwdz 


' |ti;|or|2;|<l 


u)J—eiw, z)xh{z, v) 


277 277 ' 


dwdz 


(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 


In order to evaluate I 2 we need first to compute (the second derivative of) Q{w, z). This happens 
to depend on the domain of w, z. When |t(;| and 1^1 are both greater than one, we can use formula 
(2.33) of the letter [23]: in particular, in this domain z) = 0{g). In the remaining 

domains (i.e |tc| and \z\ are not both greater than one), we have formula (4.4): 


d^ 

dwdz 


Q(w, z) = 277——d(w — 
dw 


2 ) + 


1 1 

y/2g dz'^ z — w 


+ 0(1V). 


(4.35) 


Using this formula we can estimate We have 


jrest 

-'2 


' |ui|<l 


dw , . d , . 

-XHyw,u)—XHKw,v) - 


V2g 


277 '^^"^ ’ dw' 

dw dz , , „ 1 

u)-^P- -r 

277 277 dw‘^ z — w 


VlTI<i 


Xh{z,v) + 0{l/g‘^), 


(4.36) 


where we used the fact that xh{w,u) is 0{l/g) when |r(;| > 1. The first term in (4.36) cancels Ij. 
The second term equals 


2y/2 


dw dz d? 1 

tt^XhKw, u)-—P- -^ 

27r 27r dw^ z — w 


{Xf{u)-Xf{v)) 


Xh{.z,v) = 
l+Axf{u)xf{y) Xf{u)xf{v) 


g “ " l — Axf{u)xf{v){l — Axf{uY){l — Axf{vY) 

In order to get this result, we made use of the approximation 


d . , w — 2xf(u) 1 

-t^Xhkw.u) = , -Z- -^- 

dw \/l — w'^{w — u) 2y/2g{w — uY 




(4.37) 


(4.38) 
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and of the integral (B.31). For what concerns we have 

If dw f dz 1 

V2g 27 r y|,-|>i 27 r^ A _ iXf{u) -x{w) 

V 


I 2 = 


( 4 ( z+l \ 4 \ ( H±if 4 / z -1 \ ; 

ViD+l/ \z-lj ViD-l/ \z+lj 


Xf{v) 


W — Z 


^yrrT%(h)-x(f) ■ 


(4.39) 


Now, we use the identity 


4 »-l t 4 ( Z + l \ 4 _|_ t + 4 ( z-i ) 4 

iD+lJ \z-lJ Vii-l/ Vs+l/ 


1 + \ h- ~ \ h- ~ ^ 


W — Z 


and arrive at 


= 


x{w)-x{z) 

dz 1 ^/(h) 


(4.40) 


1 /■ dw j 

V2g 27r w^l - 


1 + \ h- ~ \ h- ~ ^ 


Xf{v) 


x{w)-x(i] 


(4.41) 


We now use the symmetry properties of the integrand under the exchange w with 2 and factorise 
the integral as 


l2 = 


2V2g 


Xfiu)xf{v)[xf{v) - Xf(u)]2{u,vf, 


where 


'2{u,v) = 


dw 


^ \h- ~ ^ 


'h>i 27 r {xf{u) - x{w)){xf{v) - x{w)) 

We change variable of integration from tv to x{tv) = y. We get 


(4.42) 


(4.43) 


!I(M,n) = 


dy 


\>i/2 27 r 2 / Y 1 - 4^ ^/(“) - y - y ’ 

which can be exactly computed by means of (B.31): 


(4.44) 


^(m,^) = 


^/2 


Xf{u) - Xf{v) 


\/l — 4:Xf{uy -^/l — 4:Xf{vy 


(4.45) 


Therefore, we obtain 


l2 = 


1 Xf{u)Xf{v) 

V 2 gXf{v) -Xf(u) 


Vl - 4x/(n)2 y/1 - 4a;/(M)2 


(4.46) 
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Adding (4.37, 4.46) we arrive at the final formula 




exp 


2iV2 


{xf{u) - Xf{v)) 


l + 4 :Xf{u)xf{v) 


X 


fiu)xfii 


i Xf{u)xf{v) 

y/ 2 gXf{v) -Xf{u) 


4xf{u)xf{v) (1 - 4x/(m) 2)(1 - 4x/(n)2) 

2 

1 


— 4x/(n)2 — 4:Xf{uy 


+ o(W 


(4.47) 


Fermions in the giant hole regime 

In the giant hole regime which fits into the large fermion case the fermion rapidity scales as 
xf{u) = \/2gx{u), u = \/2gu, with 


x(u) = 


u 

2 



m| > 1 , 


x{u)\ > 


1 

2 ’ 


(4.48) 


where x{u) was already defined in the second of (4.9). Referring then to formula (C.ll) for large 
fermions, we first show that ii w = w\/2g, u = u\/2g and g +cxo, with |t()| > 1, then 


Xf{.w, u) + ^{w) = nsgia{w — u) — 7rsgii{w) + 0(1/5^). 


(4.49) 


In order to prove (4.49), it is convenient to start from (A.9) and then use (4.9). Therefore, the 
situation is completely analogous to the gluon case: relation (4.49) implies that in (C.ll) the 
dominant contribution comes from integrations in the second term in the region |r7| > 1, |x| > 1, 
where the scalar-scalar factor 0 is proportional to g. The final result is 

= Q{u,v) + 0{g^) = iln(—n)) + 0{g^) = + 0{g^). (4.50) 


Fermions in the non-scaling regime 

We send g —?■ -|-oo by keeping fixed all the rapidities. For fermions, rapidities are the variables 
x: therefore if we keep x fixed, we are necessarily in the small fermion case, i.e. \xf \ < g/V2. We 
can then show that 


lim xh{u,v) = Xf{v) 

g—^+oo 


2^2 

9 


1 


0{i/g^) 


(4.51) 


This means that in this regime ihiS^^^\u^v) = 0(1/5'^). 


4.4 Mixed factors 

Scalar-gluon 

• Perturbative regime 

We start from the exact expression (C.18) 

*ln[R^®^)(M,v)] = x{u,v\l) + $( m ) - j ^^{u,w)[x{w,v\l) + $(«;)], (4.52) 
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where both the scalar and the gluon are in the perturbative regime. This means that the u rapidity 
is parametrised as (4.2) and the v rapidity is scaled as n = \/2gv. In these hypothesis we have 
that 

X{u,v\l) + $(«) = 

In addition, in first approximation, we can integrate in the region |t(;| < ;^ Inm and use for 0 the 
expression (4.1) in which rapidities are parametrised as (4.2). The hnal result is 


\/2gv - sgn(M) 


(4.53) 


*ln[^(*»)(n,n)] 


1 1 
27iy/2gv - sgn{u) 


^ ^np 



(4.54) 


where |m| = ^ In m{g) < z < 1 being the scalar rapidity in the perturbative regime. 

• Giant hole regime 

In the giant hole regime we use formula (4.17) for the limiting expression of x(n, ■u|l)+<l)(n) when 
both |n| and |h| are greater than one. Then the leading (i.e. 0(g)) contribution to i\n S^^^\u,v) 
comes from integration in the second term of (C.18) in the region |t(;| > 1. This integration is 
easily done and the result is 


Gn^(*^^)(M,n) = *ln^("5)(n,n) = 0(n, n) + 0(^°). 


(4.55) 


• Non-scaling regime 

In order to compute the scattering phase S^^^\u,v) in the non-scaling regime, we have to plug 
the expressions (4.4) and (4.20) into (C.18). Since x(M,n|l) -|- $(«) is of order 0{l/g‘^), we claim 
that 

ilnS'^^^\u,v) = [x{u,v\l) + ^{u)]- j ^ ^Q{u,w) [x{w,v\l) + ^{w)] = O . (4.56) 

Gluons-fermions 

• Perturbative regime 

We study the scattering factor between gluons with rapidity u and (small) fermions with ra¬ 
pidity Xf{v) in perturbative regime of the strong coupling limit, i.e. u = \/2gu, with |h| < 1 and 
Xf{v) = \/2gxf{v), V = V^gv, with |h| > 1, \xf{v)\ < 1/2. 

We start from 

*ln {-S^30^u,v)) = h + l2 + h, (4.57) 

where 

h = 2 arctan 2 (m — v ), (4.58) 

/•+°° rlw d 

h = J —[x{w,u\1) + ^{w)]—xh{w,v), (4.59) 

^3 = - / u\1) + z)xh{z, v) . (4.60) 
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In the perturbative regime 


1 


Ji = 7rsgn(M — v) 


\/ 2 g{u — v) 


+ 0 {l/g^), 


h = f ^[x{w,u\l) + <l>{w)]^XH{w,v) + 0{l/g^), 


' |u)|<i 


27r 


dw' 


(4.61) 

(4.62) 


h = 


dwdz 1 + %(«) 


'|tu|,| 2 |>i 27r 27r M w 
f dw 


w — z 




d 


|ui|<l 


27r 


[x(w,m| 1 ) + $(«;)]—^//(w,^) + 0 ( 1 /^ ) 
dw 


(4.63) 


We evaluate the sum I 2 + I 3 by hrst performing integration in w with the help of (B.33), getting 


h + h — 


Xf{v) 


dz 


( l+n "\ 4 ( -z-l f 4 _ ( l-n 'j 4 ( Z+1 \ 

Vl-il/ V 5 + 1 / Vl+S/ \z-lJ 


u — z 


2^ J\2\>1 27r ^ Xf{v) - x{z) 

Then, we integrate in z, using (B.34). We obtain 

l l-2xf{v) , l l+2xf(v) n-jz\i /7T 

1 Y l+2xf(v) \1—S/ Y l—2xf{v) Vl+M/ ^ 

h + h = -=-=- 

Ag V — u 

Summing up I 1 + I 2 + I 3 we get the hnal result 

P) \ [TE^Ei. ( i+hi 4 _i_ / i+2xj^ / inHi 4 

i Y l+2xf Vl-S/ Y l-2xf Vi+m/ 


+ 0 {l/g^ 


+ 0 ( 1 //). 

(4.64) 

(4.65) 


S^3f)^u,v) = 


exp 


4^ 


u — V 


+ 0 ( 1 // 


(4.66) 


For what concerns ilnS'+0(M, w) = I 2 + I 3 , we have 
S'*'^+(m, v) = exp 


l+2xf 




^9 


u — V 


+ 0 ( 1 // 


(4.67) 


• Giant hole regime 

Since both x(v,m|1) + $(f) and xf{v,u) + $(f) have the same limiting nonzero expression 
(4.17) when |7| > 1, Ihl > 1, the leading expressions for i ln(—i ln(S'^®0(^y^ 
coincide with the one for i\n S^^^\u,v). Therefore, 

t \n{-S^^^\u, v)) = i ln(^(50(^^ y)) + o{g°) = Q{u, v) + 0{g°). (4.68) 


• Non-scaling regime 
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As written before, in this regime fermions are necessarily small. Then, since x(tc,M|l) + <h(tc) 
is Oil/g"^) and xh{w,v) is 0{l/g), integrals hyh are both 0{l/g^). For what concerns Ji, since 
fermionic rapidities Uf^k are bounded by the inequality u‘jj^ > 2g^, we can safely approximate 
_g*/i = 1 -|- 0{l/g). Therefore, in the non-scaling regime S^^^\u,v) = S^^^\u,v) = 1 -|- 0(1/5^) 
and v) = n) = 1 -|- 0(1/ g^). 

Scalars-fermions 


• Perturbative regime 

We start from the exact expression (C.16), 


- t In v) = xh{u, v)- j ^)Xh{.w, v) , 


(4.69) 


and make the parametrisations (4.2) for u and v = \/2gv. At leading order we have 


Xh{u,v) 
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1 

h — 1 -1- — 1 


(4.70) 


Then, as for the scalar-gluon case, in hrst approximation, we can integrate in the region |tc| < 
^\nm{g) and use for 0 the expression (4.1) in which rapidities are parametrised as (4.2). The 
hnal result is 


— i In v) 


i In m) = 2 4 



1 

h — 1 -|- \/h^ — 1 



(4.71) 


with ImI = - In —fr- 

I I 7r m(g) 

• Giant hole regime 

Since Xf^v, u) + <h(n) and x(n, m| 1) -|- <h(n) go to the same limit, i.e. 7rsgn(n — u) — 7rsgn(n), by 
comparing (C.16) and (C.18) we get the equality GnS'^^ ^Ku,v) = —i In S^^^\v,u) + 0{g^), which, 
together with (4.55), gives 

iln^(*^)(M,n) = 0(M,n) + 0(^°). (4.72) 


• Non-scaling regime 

We perform the non-perturbative limit of the scalar-fermion scattering phase 

ilnS^"^\u,v) =-xh{u,v) + j ^ Xh{w,v) , (4.73) 

by taking g —> oo while keeping the scalar rapidities hnite, whereas the modulus of the fermionic 
rapidities Xf must be \xf \ < g/\/2. Under these assumptions, we can make use of the approxima¬ 
tions (4.51) and (4.1) for xh{u,v) and Q{u,v)-, eventually, we hnd: 

Gn^("^)(M,n) = O . (4.74) 
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4.5 Remark on the non-scaling regime 

We showed that in the non-scaling regime all the factors S**'{u,v) go as 1 -|- 0 ( 1 / 5 ^^), with the 
exception of the scalar-scalar one which goes as 


otel/ ^ r ( 1 + ;!!=!!) r (i - ;!!=!!) 

5' '(u, v) — - ^^ -^exp 


, vr 


[1 + 0 ( 1 /^)] 


(4.75) 


In addition to that, we recall that the fermionic rapidities Uf^k satisfy the inequalities u^fk> 2g‘^. 
Therefore in this regime all the rational factors involving fermionic rapidities (which appear in the 
quantisation conditions for fermions and in the equations for isotopic roots Ua and Uc) go to one. 
In addition all the exponentials of momenta and defect (= e*(P + D)) go to 1, with the exception 
of those for scalars. Summarising, in the limit g —)■ -|-oo with rapidities hxed (and hnite), the 
non-trivial equations are 


1 _ gii?pWK)+iDWK) TT ^ - ^b,j + 2 -|-r ^ 

fJi un - u,, - ^ 5 A 

h'^h 

^ _ "TT ^ "TT ~ 

j^k ^ "^CL^k ^b,j i 

H K 

n Ub,k — Uh + i/2 _ y-T Ub^k — Uag 

gij.i. — 111. — i lO A A 


K. 


^/2 1 I" Ub^k ^c,j 


i/2 -1 r Ub^k ^b,j T i 


'^b,k Uh i/2 


J =1 


UbM 


n 


'^/2 y^b^k ^cj i/2 Ul)_^k "^bj i 


n 




I • ■ /Q 

1 1 r ^ 1 r ^c,/c ^b,j 1 1 2, 

XX^ , — 1 ! . — ^xx 


No 


'.,k Uc,j i Uq ^ Ub,j T i/2 


YTn-^j + ^ 

J^k ^ 1 

^u\-vP,+i 

i=n-i—i—' 

fit - “3 -' 


(4.76) 


with given by (4.75). Since equations for gluons have no solutions for hnite rapidities (i.e. 
Ng = Ng = {]), equations (4.76) show that in the non-perturbative regime the only active excita¬ 
tions are the six scalars. The other excitations are obliged to assume inhnite rapidities and thus 
decouple to very high energy from the scalars. The latter satisfy the above ABA (4.76) which is 
the same we can derive from the 0(6) non-linear sigma model S'-matrix of [53]. Therefore, also 
the exact TBA would be that of the 0(6) model (if we can neglect the exchange of the g —)■ -|-cxo 
limit with the thermodynamics). 
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5 Particle momentum in different forms 


Momentum was already thoroughly discussed by Basso in [20]. The aim of this section is to show 
that the expressions for momenta of the various excitations we found (in our notations) in previous 
sections agree with corresponding formulae of [ 20 ]. 

• Scalars 

We found (2.74) that the momentum of a scalar excitation is 

P^^\u) = 2u —-Zbes{u) ■ (5.1) 

Now, using the mapping (B.19), valid for k > 0, 


- — - ; -, 

TT k 

between our quantities and quantities used in [20], we immediately write (5.1) in the form reported 
in [20] (second of the (4.6)). 

• Gluons 

For a gluon with rapidity u we found for the momentum the expression (2.93), which we can 
write in Fourier space as 




^+oo 


dk 


-27r 


sin ku 
k ' 


-FI 




Tl=l 


g \ 271 _ 1 M T , /K 

~7n —FT ' Jni^gk) + 

iy/2x~(u)J k 


I . /K , . ^ ' 2 'Jn(\/2^/c)l(-47r(5(/c) + ^(JBE5(fc)) = 

\iy2x+[u)/ K J 2 


x~ x^ 


+00 


"+00 


dk 




\iy/2x+{u)J k 

n odd 


—00 

n 


- 


^ sin /cm 

—27r—;—e ' ' + i 


+CXD 


>E 


72 = 1 

n odd 


\ ' 2 'J„( V2^/c) + 

iyj2x~{u)j k 


2 e _iq /— 1 

e 2 J^{^/2gk) asEsik) 


Now, we use the equalities 


r‘+oo 

—e~~ Jn{'/2gk)dBEs{k) = AE[V2g6n,i - it], n odd,, 


and 


E 

72 = 1 

n odd 

to eventually obtain 


\/2x J \i\/2x^ 


. dk . ^ _k ^ 

In = / —Sill /cm e 2 2gk ), 

Jo k 


.,0 


(5.3) 


(5.4) 


= 2u 


/•+°° dk 

/ sm /cm e 2 

-ft{\/2gk) -^1{y/2gk) 

lo k 

1 — e~^ — 1 


(5.5) 
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which agrees with the second of (4.9) of [20]. 

• Large fermions 

The momentum associated to a large fermion with rapidity u enjoys the expression (2.82). In 
Fourier space it reads 




= 2u — 


sin ku _1M 
k ® 

/•+“ dk , 

' — sm ku 


E 


9 


i\pl. 


x[u\ 


n _1M 

e 2 


ik 


Jn{V2gk) 


n=l 

-iX{y/2gk) +7®(a/25(/c) 




'0 


1 

2 


— sin ku'^t{\/2gk) , 

1 rC 


k — 1 

which recalls the second of (4.10) of [20]. In order to obtain the equation in the last line, we made 
use of the relation 

2 n-l 

7;„-i . (5.6) 

\ A / '/g'y'i 91 1 / 

n=l 


[ ^ sin ku 7 ® {V2gk) = i ^ - 

Jo 2fc ^AV2tx{u) 


which holds for u"^ > 2g‘^. 

• Small fermions 

The reasonings for the momentum of a small fermion with rapidity u mimic very closely the 
large fermion case. We start from our expression (2.87) and in Fourier space we eventually get the 
result (m^ > 2 ( 7 ^) 

1 dk 

P^^\u) = - / — sin{ku)'^t{'\/2gk), (5.7) 

2 Jo ^ 

therefore matching the second of (4.12) of [20]. 


6 Strong coupling analysis of the defect term 

We now perform a quantitative analysis of the strong coupling limit of the defect which appears 
in the Bethe equations on the GKP vacuum. 


6.1 Scalars 


It is convenient to concentrate on the function Z 4 (2.69) in absence of excitations, which equals 




( 6 . 1 ) 


The study of this function, which relies also on previous results, provides information on both the 
momentum and the defect of the scalar. 


Non-perturbative regime 


In this regime we send g —?■ + 00 , keeping the rapidity u fixed. We can use results from [21] 
where the non perturbative regime for the pure s/(2) sector is studied. We found that the function 
( 6 . 1 ) has the form 


Z 4 {u)\]yE = -‘^m{g) In 


2V2s 

9 


TT 


sinh —u + 0{m{g )^), 


( 6 . 2 ) 
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where m{g) is given by (3.13). Therefore, the contribution of the two defect is proportional to the 
momentum 

D^^\u) = m{g) sinh. + 0{m{g)^) , (6.3) 

so allowing us to fully re-absorb them into a re-definition of the size Ri^g) as in [55, 21] 

Perturbative regime 

We obtain interesting formulae in this regime by plugging (4.2) in (6.3), namely: 

D(*)(M) = ^ln^ + ... , (6.4) 

^ 9 

for right movers (m > 0), while 

Z1(^)(m) =-^ln^ + ... , (6.5) 

^ 9 

for left movers (m < 0), where evidently the dots ... mean sub-leading corrections. 

Scaling regimes 

We introduce the density (Ti{u)\NE = ■§^Zi{u)\]siE and rescale the rapidity u = y/2gu. By using 
techniques developed in [15], we eventually find 


/ r -M dh 

o'4:[v2gu)\NE = / ^ cos hu 


+ 


Jo V2g 
271 

cosh 7r\/25f-u 


_ h h 

p 2y/2g _ p 2\/2g __ 

j-T4h) +—-r7-Ch) 






41ns -|- 


+ 2'ip + iuy/2g^ + 2'ijj ^ - iuV2g^ + 


( 6 . 6 ) 




where the functions r± satisfy the relations, valid when |m| < 1: 
dh 

Jo V2g 


■ sin hu[T_(h) + r_|_(h)] = —2i[-^(l — iuV2g) — -^(1 + iuV2g)] 


(6.7) 


b+co 


dh 


cosh.M[r_(h.) — r+(h)] = 41ns — 2['0(1 — iu\/2g) + -0(1 -(- iu\/2g)]. (6.8) 


Jo V2g 

Going to the strong coupling limit g —)■ -|-cxo, with u fixed, we find 


a4{V2gu)\NE = 


b + CO 


dh 

VTg 


cos hu 


_ h h 

p 2\/2g _ p 2\/2g _ 

T^{h) + ^^T.{h) 


cosh 

V2g 


cosh 

V2g 


41n- + 0(^°), (6.9) 
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where 


0 + 00 


dh 


sin hM[r_(h) + r+(/i)] = 0 (( 7 °) , |m| < 1 , 


Jo V^g 

r+°° dh - - - s 

/ —^ cos h'u[r_(/i) — r+(h)] = 4 In—h 0((?°), !«! < 1. 

Jo v'^g g 

Solntions to these eqnations go differently according to the valne of Ihl. If !«! < 1 we have 


( 6 . 10 ) 

( 6 . 11 ) 




-S(u) + 0(1/+), 


( 6 . 12 ) 


which means that P^^\u) is exponentially small and that = —^^6{u) + 0 ( 1 / 5 '^). 

On the other hand, if Ihl > 1, we have 

(^4{u)\ne = - 21 n -^ p (^)( M ) + 0{g^) , 
g du 

with = 0{g) or, alternatively, the proportionality to the momentnm 

D^^\u) = -P^^\u) hig + 0{g), 


(6.13) 


(6.14) 


which would allow us to re-absorb, at this order (only), fully the defect into a simple redefinition 
of the size R. Using now (5.27) of [20] we can then express in terms of the rapidity 

x{u) = \u + \u\J 1 “ ^ (defined in the second of (4.9)) 


as 


D^^\u) = —V2g\ng 


2 x{u) -i /1 — „ - ' — arctan [ 2 x{u) \ I — 


Ax‘^{u) 


4:x‘^{u) 


+ 0{g). (6.15) 


6.2 Gluons 

Perturbative and giant hole regimes 

We start from the formula 

/ I 1 

^[x{v,u\l) + x(-n,M|l)] -a4(n)|jvE, (6.16) 

Scaling u = \/2gu and v = \/2gv, we have that x(n,M|l) + $(n) is 0{l/g) if |h| < 1 and (at 
most) 0{g^) if Ihl > 1. Referring to (6.12, 6.13), we remark that the integration receives leading 
contribution from the region |h| > 1. We conclude that 


|h| < 1 ^ 2 \n sP^<^\u) + 2D^<^\u) = 2\n-P^<^\u) + 0{g°), 

g 

|«| > 1 => 21nsP“(!i)+2C“(M) = 21tijP“(M) + 0(+. 


(6.17) 

(6.18) 
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We can now refer to formulae of [20] for the momentum of the gauge held and arrive at the hnal 
expressions 


|h| < 1 ^ = - In gP^>^\u) + 0{g°) = 


\ng 

'V2 


1 + u 
1 — u 


1 — u 
1 + u 


(6.19) 


|h| > 1 ^ =-\ngP^3\u)+ 0{g) = 

-V2g\ng 


- 4i^ - 2*(«) Jl - 


+ 0{g), (6.20) 


x(u) being dehned in the second of (4.9). 


6.3 Fermions 

Perturbative regime 

In this regime the scaled rapidity needs to satisfy \xf{u)\ < 1/2, so that it belongs to a small 
fermion. Therefore, we start from the formula 

/ rji) 1 H 

7^[Xh{v,u) + xh{-v,u)\ 2 ^^ 4 ( w )| vi <;- ( 6 . 21 ) 

Scaling v = \/2gv, we have that xh{v,u) is 0{l/g) if |h| > 1 and 0{g^) if |h| < 1. Referring to 
(6.12, 6.13), we remark that the integration receives leading contribution from the region |n| > 1. 
We conclude that 

2\nsP^^\u) + 2D^f\u) = 21n^ P^^\u) + 0{g^) (6.22) 

and 

D^P[u) = -\ngP^^\u) + 0{g^) = -2\ng^=^i^== + 0{g^). (6.23) 

v/l-4xj(h)2 

Giant hole regime 

On the contrary, here the scaled rapidity is |a:i;’('u)| > 1/2, so to characterise a large fermion. 
Therefore, we start from the formula 


2\ns P^^\u)+ 2D^^\u) = f ^[xf{v,u) + XF{-v,u)]^a4{v)\]yE ■ 


(6.24) 


Scaling v = y/2gv, we now have that xh{v,u) is 0{g^) for all |h|. Referring to (6.12, 6.13), we 
remark that the integration receives leading contribution from the region |h| > 1. We conclude 
that 

2\nsP^^\u) + 2D^^\u) = 21n- P^^\u) + 0{g) (6.25) 

9 

and, consequently, that 

D^^\u) =-In gP^^\u) = 


= —\/2g\n. 


9 


2x{u)a 1 


4x^(m) 


arctan 2x{u)a 1 


4:x'^{u) 


+ 0{g). (6.26) 
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7 The SU(4) symmetry 


The particles we are addressing to (scalars, gluons, fermions and anti-fermions) belong to a specific 
multiplet under the SU{A) symmetry ( 6 , 1, 4 and 4, respectively). In fact, the scattering matrix 
possess this symmetry. Starting from the scattering matrices derived in the previous section, the 
Bethe equations may be assembled for every sort of excitation; anyway, they are actually able to 
catch only a single state in each multiplet, precisely the one corresponding to the highest weight 
state of the representation. In this section the focus moves to a few sectors of the complete theory, 
which include just one type (or two at most) of excitations along with the set of isotopic roots, 
aiming at elucidating the behaviour of the different kinds of particle under SU{4:). 

Following [45], a set of Bethe equations can be formulated for any spin chain associated to 
a simple Lie algebra. Therefore, given the set of simple roots of a simple Lie algebra{aq}, and 
chosen a representation p by hxing its highest weight Wp (or equivalently a tern of positive integer 
Dynkin labels), the relative Bethe equations arise, with a further generalization stemming from 
the introduction of a set of inhomogeneities (labelled by their rapidities Um) along the spin chain; 


1 r T * 'iCp 

lOCq * Wp 


K, 

n 


"^q-ij lO^q * O^q 






'^q,k t(y.q * CXq' 

^q,k ^q* ij lOtq * Otqf 


(7.1) 


Turning to the sm(4) algebra, we perform a choice of three simple roots along with the three 


simple roots (p^, resulting from the dehning condition —4 


2 Q!j • ip}^ 


= 5, 


[04 


kj 


CKi = 


02 = 


Q!3 = - 


i-f.o 


2’2,/3’ Vey 


¥^ 1 = 0 


2 ’ 2x/3’ 2^6 


¥^2 = - 


1 \ 




= 0 , 0 , 


2 ’ 2\/3’V6/ 

3 




(7.2) 


To sum up, the Bethe equations in (7.1) specialize to the sm(4) algebra: 


1 r ^a,fc '^m T 'ioi • Wp 

"^a^k ioti * Wp 

TT — Um + i02 ' Wp 
Ub,k — Um — i02 ■ Wp 

1 r ^c,A: “ 1 “ ^^3 * 

'^c,k * Wp 

m ’ 


I ■ ■ /o 

1 r "^a^k 1 1 r y^a^k ^b,j %j2, 

1 r yb^k ^b,j “i“ ^ T r yb,k 2 T r yb,k %/2 

■ jj yb,k ^bj ^ yb,k "^CLij ^/2 yb,k ^c,?' ^/2 

j^k ’ j=l ’ ' j=l ’ ' 


Kc , . Kt 

TT yc,k ~ ycj + ^ TT 

“^r yck — yci — i 

j^k j=l 


yc,k ^b,j '2'/2 
yc,k ^b,j “1“ %l2 


(7.3) 


• Scalar sector: 
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When considering a system composed only of H scalar excitations with rapidities together 

with Ka roots Ki, roots Uf, and K^. roots Mc, the equations for the isotopic roots (2.49, 2.50, 
2.53) take the form : 


1 


H 


n 

h=l 


/ Ub,k - Mfe + I \ 

\'^b,k '^h 2 J 


1 


Ka .Kb 

n Ua,k — Ua,j + I -TT Ua^k 

11 1 — 'll ■ — ■? A A 


Ub.j 


Kb 


^a^k 


“J 


* Ua,k 


- Ub,j 

Ua,j ■ 


U. 


TT ~ '^b,j + 1 TT '^b,k 

“4 Ub,k — Ub,j — * 7“ Ub,k 

Kc \ ‘ ^ 

n '^c^k ^c,j ~r ^ 1 r 2 

_ , — 7 XX 


2 

2 

2 TT ~ ~ 2 


n 




u 




+ - 
^ 2 




^c.k 




* j=i W,fc + 2 


(7.4) 


A comparison with (7.3) promptly reveals that the equations (7.4) coincide with those for a spin 
chain associated to the antisymmetric (6) representation of sm(4), whose highest weight is, ac¬ 
cording to our convention, w% = 02 - The H hole rapidities (m/i, h = 1,..., 77) can be read as 
inhomogeneities along the spin chain, and their dynamics are regulated by the equations (2.73), 
suitably adapted to the case at hand. 


• (Large) Fermionic sector 

Let us stick now to a system composed of Np large fermions j = l,...,iV^, together with 
Ka roots Ua, Kb roots Ub and Kc roots Uc- While the fermions satisfy Bethe equations (2.81), the 
auxiliary roots obey relations: 


Nf 

n 


'^a,k ^F,j T 2 


\'^a,k ^F,j 2 


Ub,j 


1 = 


1 = 


Ka ■ Kb 

TT Ug^k — UgJ I T-r Ug^k 

Ua,k — Ua,i — * Ug^k — Ub,j + | 

Kb ■ Ka 

JJ- Ub^k — Ubj -|- I J-J- Ub.k 


(7.5) 


u. 


a,3 


1 Ka 

2 TT “^>7 
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M, 




^b,k ^b,j ^ - ^ '^b,k ^a,j T 2 7=1 ^b,k '^c,j T 2 


Ka 

n 


^c.k 


U 


c,J 


^c.k 


U. 




3 

■ Kb 

n 


^c,/c "^bj 2 


^ j=l ^c,k ^b,j 2 


A look at (7.3) suggests the equations (7.5) should be associated to a spin chain related to the 
fundamental representation (4) of sm( 4) (with highest weight W 4 = 0i), where the large fermions 
behave as inhomogeneities, with rapidities upj, j = 1,..., Np. 

Otherwise, when only large antifermions (in number of Np) appear in the vacuum, again 
accompanied by Kg isotopic roots Ug, Kb roots Ub and Kg roots Ug, the system is described by the 
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set of Bethe equations (2.84) together with the isotopic roots equations: 




n 


Uc,k ~ Up j + I 


\'^c,k ^F,j 2 


Fa • F}) i 

^ _ TT ~ Uaj + I TT Ua^k ~ 

jjLj^ — Ua,j — i Ua,k — Ub,j + | 

a % % 

^ Ub,k ~ ^ YJ ~ ~ 2 TT ~ ~ 2 

j—l ^b,k "^bj i '^b^k ~l~ 2 j—i 2 

Fc • Ff, i 

TT '*^c,fc ~ 1 TT '^c,k ~ UbJ — 2 

Uc,k ^c,j * Uc^k ^b,j 4“ 2 




(7.6) 


The (7.6) are in fact the equations for 4 spin chain (highest weight wi = (p^), as may be read from 
(7.3). 

Some interest should be paid to a system including both Np (large) fermions and Np (large) 
antifermions; in this case, the isotopic roots satisfy the relations: 
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i=i 


^ "^a^k ^FJ H“ 2 

\'^a,k ^FJ 2 
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iV_p 
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/ Uc,k ~ '^F,j + t 
\Uc,k ^F,j 2 


Ka , ■ Ki i 

TT Ua^k — Ua,j + * TT '^a,k “ “ 2 
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(7.7) 
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Ub,k — Ubj + I TT Ub^k — Uaj 
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1 7, 

2 T r 
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^c.k 


U. 


cj 


• A'b i 

+ I TT Uc^k Ub^j 2 




On the basis of sn(4) simple roots and fundamental weights (7.2), we can claim that the equations 
(7.7) are associated to a spin chain, related to the representation of sm(4) whose Dynkin labels are 
(1,0,1); in other terms, we found as its highest weight tcis = pi + and that leads to the 15. 
The reason lies in the way how fermions (in the 4) and antifermions (in the 4) scatter, since the 
process can be decomposed into two channels, according to the rule 


4(8)4 = 10 15 ; (7.8) 

the singlet 1 channel is not explicitly appearing in (7.7), but it can be revealed upon imposing 
some constraints on the isotopic roots (see next section). 

We eventually remark that analogous considerations and expressions hold if we replace large 
fermions/antifermions with small fermions/antifermions. 

• Gauge field sector 

When only Ng gluons (with rapidities up are excited over the vacuum, the isotopic roots decouple 
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from them, since 
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1 


1 


Ka 

n 

ya,k 

Ua,j 

+ i 

Kb 

n 

11 

ya,k 

- UaJ 

— i 

11 

i=i 

Kb 



+ i 

Ka 

n 

'^b,k 

- Ub,j 

n 

11 

i=i 

yb,k 

- UbJ 

— i 

11 

i=i 




+ i 

Kb 

n 

yc,k ~ 

- 'blcj 

n 

11 

Uf~ - 

- Uc i 

— i 

11 

j¥=k 




i=i 


Ua,k UbJ 2 
'Ua,k ^b,j “1“ 2 


^b,k 

^b,k ^CL,j 


i K, 



2 j=i 


'^c,k "^bj 2 
'^c,k '^bj “ 1 “ 2 


'^b,k ^cj 2 

^6,A; ^c,_j 2 


(7.9) 


Therefore, gluon excitations behave like singlets (1) under SU{A). The very same reasoning applies 
to barred-gluons. 


8 Eigenvalues 


While commenting on the equations (7.7), we hinted the role the S'7/(4) symmetry takes in the 
scattering between fermions and antifermions. Now we are going to examine in some more detail 
several scattering processes involving different kinds of particles. In general, given two types of 
particles a and /3, transforming under the representations of sm(4) pa and p/?, which act respectively 
on the spaces Vk and Vg, their scattering decomposes according to the Clebsch-Gordan rule 

Pa®Pii = ^^pn- ( 8 . 1 ) 

n 

Recalling [45], the scattering matrix (dehned on Va<^Vp) between excitations a and /9 enjoys 
the spectral decomposition 

= ( 8 . 2 ) 

n 

where are the eigenvalues of the matrix relatives to the (normalized) eigenvectors Pq, 

which act as projectors onto the space Vh, i.e. PniVa (8) Vg) = Vq. In this section we list the 
eigenvalues corresponding to the scattering between excitations on the top of the GKP string. 

Scalar-scalar 

The scalar-scalar scattering was completely clarihed in [22], here we list eigenvalues and corre¬ 
sponding isotopic roots; since scalars belong to the 6 , the decomposition follows: 

6G)6 = 1©15©20 . (8.3) 


The singlet 1 channel involves two type-6 isotopic roots, which shall be related to the hole rapidities 

Uh, Uh' according to = \ {^h + Uh' - and Ubp = I (^h + Uh' + 

together with a and c roots Ua = Uc = These constraints on the isotopic roots lead us to 

the eigenvalue 







yh yw y}i yfi' 


(8.4) 
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where the scalar factor S^^^\uh, Uh') can be read from (C.IO). The adjoint channel hi = 15 requires 
one 6 -type root, satisfying ui, = and no a nor c roots Ka = Kc = 0. Eventually, the resulting 

eigenvalue follows 

Si^^\uh,Uh') = — — S^''"\uh,Uh') . ( 8 . 5 ) 

Uh - Uh' - I 

Finally, the 12 = 20 channel request no isotopic roots {Ka = Kh = = 0), so that the eigenvalue 

simply coincides with (C.IO) 

S2o\uh,Uh') = S^'‘''\uh,Uh') . ( 8 . 6 ) 

Fermion-fermion 

In the fermion-fermion scattering, we have two eigenvalues corresponding to the decomposition 
4(8)4 = 10©6. The hrst one, for 12 = 10, corresponds to no isotopic roots and therefore it holds 

S[o^\uF^l,Up^2) = S^^^\uf^i,Uf^2) , 


where the scalar factor corresponds to (C.ll). The second one, for the 12 = 6 channel, is obtained 
from the solution with Ka = 1, K^ = 2, Kc = 0, such that Ua = while Ub,i = uf,i and 

Ub ,2 = Up, 2 ] consequently, we find that: 


qiFF) 

'^6 


{uf, 1 , Up,2) 


Uf,i-Uf^ 2 + 1 c(FF)r ^ 

-: S'- ’{uf,i,Uf^2) ■ 

Up,I — Mf ,2 — * 


Fermion-antifermion 

As previously mentioned, the fermion-antifermion scattering is associated to the decomposition 
4 (g) 4 = 15 © 1. Turning to the 12 = 15 channel, no isotopic roots are involved, therefore the 
eigenvalue equals the scalar factor: 




(8,7) 


The singlet channel instead is obtained from the solution with Ka = Kb = Kc = 1, where the 
isotopic roots satisfy the constraints = jUf,i + \uf i, Ub,i = ^up^i + ^up i and Mc,i = + 

As a consequence, we obtain as the eigenvalue for the 12 = 1 channel: 


S{FF) 


{uf,i,Uf^i) 


Up^i — Up I + 2i 

iip^Y — Up — 2t 




( 8 . 8 ) 


The same result, i.e. 



^/ 1 T 

“/,i - “/,i - 2 * 




(8.9) 


holds for small fermions. We highlighted (8.9), since this factor is responsible of the appearance 
of a mass two particle at 5 ^ = +cx). This issue will be discussed in next section. 
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9 Classification of possible bound states 

9.1 String solutions at large size 

In the large size limit, R —?■ +cx3, solutions to Bethe Ansatz equations show many (numerical 
and analytic) evidences that they organise into strings or stacks (generalised strings with different 
isospin or nested degrees of freedom). Their derivation follows as customary [56]. Let a complex 
rapidity u\ exists, whose imaginary part be different from zero, but with sign so that the factor 
^iRP{ui) gQgg 2 ero (infinity) in the large R limit; then another rapidity u*- must exist with the 
same real part, but imaginary part lowered if Imw^ > 0 or, otherwise, raised if \m.u\ < 0, by 
an appropriate quantity, in order to drive rational factors in S**'{ul,u*') to infinity (zero), thus 
balancing the ABA equations. The process can continue by involving a finite number of extra 
rapidities displaced at regular distances until a string of m roots disposed around a real ’center’ 
is closed. Since the ’wave function’ of a string of m roots is by construction rapidly decreasing at 
±cx3, we naturally associate this configuration with a bound state of m ’elementary’ excitations. 

In this section we discuss some possible bound states, with the important caveat that the list 
below is not meant to provide a complete classification of the particles living in the theory. This 
is indeed an interesting problem in itself and will be possibly dealt with in a future publication. 

We also remark that, strictly speaking, the complexes of solutions we provide below are meant 
to be valid for finite values of the coupling constant, i.e. g ^ 0 and g < oo. At g = +cxo the situation 
ought to be different, as it can be inferred from considerations on the classical (quadratic) string 
theory action. Indeed, its small fluctuations in the bosonic sector consists of two mass = \/2 (real) 
bosons and one mass 2 (real) boson, besides the five massless bosons (of the 0(6) non-linear sigma 
model) [6, 50]. Seemingly, this mass 2 boson degree of freedom is missing in the gauge theory 
but, in the following, we find evidence that, with this mass, there is indeed a composite state made 
up of a small fermion and a small antifermion in their singlet channel. By means of the string 
mechanism discussed before, this bound state cannot exist at finite g, since divergences (zeroes) of 
the phase factor Q^RAf)(uf) small fermions for complex rapidities go together with divergences 
(zeroes) of the S matrices, instead of compensating each other: thus this is a resonance. Indeed, 
after specialising the isotopic roots to the values defining the singlet channel of formula (8.9) and 
after switching off for simplicity’s sake all the excitations (dots in the equation below) but a fermion 
and an antifermion, we have to cope with the coupled equations 

M /,1 - M /,1 - 2i 

1 = ^iRP^^'>{Uf ,)+2iD<~f'>{Uf ,) '^I,l ~ 

M /,1 - M /,1 - 2i 

Thus, we can verify that the candidate string-like solutions m/,i = um,i ± i, uji = um,i T h and 
related isotopic rapidities of the singlet (8.9), with mm,i real, cannot satisfy (at large R) equations 

^^Before us, many authors shared this concern, as, for instance, [50, 20, 52, 36]. 


A(^^)(m/,i,M/-,i)... , 

S'(-^^)(m/-i,m/,i) ... . (9.1) 
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(9.1), since > 0 if liauf < 0 and viceversa^^. However, this remains true as long as g is 

finite, while this bound state can appear as a new ’particle’ when the value of g is strictly +cxd. In 
fact, the point g = +oo is rather peculiar and singular, as complex scaled rapidities {u = u/\/2g) 
all collapse into the real axis, thus making possible a solution of ABA equations with a stack with 
two (small) fermion-antifermion rapidities (besides the isotopic rapidities): the two constituents 
possess no binding energy and no ’breathing’, but they are just ’one on the other’. This is indeed 
a new (real) scalar, named here ’meson’, coming to life only in the classical string regime g = -|-cxo 

In summary, our following analysis of the ABA scattering on the GKP vacuum shows evidence 
for the existence of this bosonic particle with mass 2 as long as g ^ -|-cxo. And not only: the same 
mechanism at g ^ -|-cxd sustains the existence of a bound state of A: = 1, 2, 3,... mesons with mass 
rrik = 2k (zero static binding energy, as well as for gluon bound states). This is a new bosonic 
sector with respect to the classical (quadratic) string spectrum (and, a fortiori, to previous gauge 
theory analyses), but yet indispensable to be considered in the BSV series for 4D amplitudes, - as 
we shall see -, for making checks with and reproducing the string minimal area solution (in other 
words the Thermodynamic Bubble Ansatz (TBA)). On the other end if this is an important way 
to check the validity of the series, it also conhrms the pentagonal amplitude values and the 2D 
scattering factors entailing them. Eventually, the formation of mesons and bound states thereof 
shows a sort of con/fnemenf phenomenon at strong coupling as for the 4D amplitudes/Wilson loops, 
in that the contribution of the constituents, the fermions, to them is subtlety sub-dominant (as 
g —-|-oo, cf. [36] and below). In fact, this negligibility is not true for the 2D scattering amplitudes 
in themselves, but in their contribution to the 4D ones. 

A more mathematical understanding of the small fermion-antifermion state ought to arrive [59] 
from the collision of the poles into the integration (real) axis [60, 61, 62, 63, 36]: this will give us 
the opportunity to explain the meson bound states and hence the confinement under a different 
light. 


• Gluonic strings: 

A hrst example is provided by strings made up of gluons or, alternatively, barred gluons, as 
equations (E.6) and (E.7) suggest. In this case one remarks the emergence of complex of solutions 
characterized by length m and real centre 

u’fr = - I - 2k'), /c'= 0, ...,m - 1 ; (9.2) 

the very same structure may be built by assembling barred-gluon rapidities, too. We will study 
more extensively bound states of gluons in next subsection, where we will show that they can be 
also obtained starting from the BMN vacuum by considering stacks of roots of type 1, 2 and 3. 


^'^We take up the occasion to highlight that we have found in the previous part of the paper that this opposition 
of signs happens only for the small fermions 

^®Only at this value its rapidity, otherwise virtual [57], enters the physical domain [58, 57]. 

^®We are particularly grateful to I. Kostov and J.-E. Bourgine for explanatory discussions on this point. 
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• More bound states: 

Along with gluonic strings, whose structure is quite ordinary, the ABA equations also admit the 
existence of more peculiar kinds of complexes whose composition and length result completely 
determined by the SU(4) symmetry of the vacuum. In fact, the structure of the Bethe equations 
for the SU(4) spin chain (7.1), reflected in the equations for auxiliary roots (E.8, E.9, E.IO), 
prevents these strings from including more than two massive roots (or exceptionally three, as in 
(9.7) below), intertwined with isotopic roots which are spaced by a constant distance fully fixed by 
(7.1). The presence of isotopic roots is necessary for these strings to effectively represent solutions 
of the ABA equations and live in some definite scattering channel (see section 8), thus behaving 
in a broad sense like bound states which belong to some SU(4) ’isospin’ multiplet. Below, such 
peculiar strings are listed according to their composition and SU(4) behaviour. 

• Bound states of large fermions in the 6 channel: 

i i 

'^k T 2 '^a,k ^k ^b,k,± ^k T ; (9-3) 

the same structure occurs with antifermions too, upon substituting fermions with antifermions and 
the central a-root with a c-root. 

• Bound states of large fermions in the singlet channeT^: 

i 


UF,k 

— UM,k + i 

'^a,k '^Mjk T 2 

Ub,k — UM,k 

Up,k 

= UM,k — i 

i 

'^c,k '^M,k 2 



(the complex conjugate of (9.4) is a solution, too). 

• Bound states of scalars in the 15 channel: 

i 

^k i 2 ^b,k ^k ■ (^•^) 


• Bound states of scalars in the singlet channel: 

Uh,k,± = Uk±z ut,k,± = Uk±- Ua,k = = Uk ; (9.6) 

it is important to point out that these strings made of holes do not survive to the strong coupling 
limit in the non perturbative regime, as they are destroyed by poles of (4.1). Indeed, it is a well 

^^Anyway we remark that bound states of this sort do not play any role in the strong coupling perturbative regime 
and, noticeably, their centres need to lie on the real axis in the region \uM,k\ < V^g, hence inside a square root 
branch cut in the large fermionic u-rapidity plane, so that perhaps they should not even be considered physical. It 
is thus far from being obvious that any relation exist with what in the following we will refer to as ’meson’ bound 
states, which exclusively subsist at g = oo and are made of small fermions, instead. 
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known result that in this regime the scalar dynamics is regulated by the 0 ( 6 ) non-linear a-model, 
which lacks bound states. 

• Mixed bound states of large fermions and gluons in the 10 channel: 

'^F,k,± ~ Uk i 2 ~ ■ (^•'^) 

(the same also holds for barred-gluons and large antifermions). 


• (Purely) Magnonic strings: 

Also three distinct kinds of massless strings made of isotopic roots only, one for each type, can be 
found: 


u 


a,k,j 


U 


B 

b,k,j 


U 


C 

c,k,j 



j = o,. 

-,71-1 ; 



J =0,., 

..,5-1 ; 

(9.8) 

+ ^{C-l-2j), 

j = o,., 

..,^7-1. 



9.2 Bound states of gluons 

On the BMN vacuum with a sea of M 4 roots bound states of excitations with rapidity 
can be constructed [20] as stacks involving type 1, type 2 and type 3 roots: 



Q.m 1 ^ 

= % +5' 

(m — 2 

1 

to 

k' 

= 0 ,., 

.., m — 2 

U2,k 

q.m 1 ^ 

= “fc + 0 ' 

(m — 1 

1 

to 

k' 

= 0 ,., 

.., m — 1 


U 3 ,k = + ^(m - 2/c'), /c' = 0, ...,m. 


Analogously, bound states of gauge helds with rapidity are obtained from (9.9), with 
g ^ g and Ui, U 2 , u^, Uq, Uf. In presence of bound states of gluons Bethe equations should be 

modihed as follows. 

Bethe equations for bound states of F+± {Ng {N^g) is the number of bound states of F+± {F+±) 

^®In addition to the string configurations listed above, several further complexes of solutions could be found, 
although strictly speaking they should not be considered actual bound states, since they are not endowed with real 
valued momenta; an example is offered by strings made of one single scalar and one fermion (or antifermion) bound 
together, whose distance gets fixed by the SU(4) symmetry to and which could be probably related to a similar 
state described in [36]. 
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1=1 j=l 


1=1 j=l 


H 


h=l 


Nr 


with length /: their centers are indicated with {uf’')) are 

+00 +CXD 

1 ^ ,.KF“(.r)+rf(.r) ') n 11 sif («r. “f') n “/.) • 

(9.10) 

i=i j=i 

where momentum and defect are given by 

= 

2Dif(u) = 


Nf 

n ^F,) n W'™. “Fj) n “/f) n «/j). 

i=i i=i i=i j=i 


+°° dv 

— [x{v, u\m) + x(-n, M|m)] 


' —CX) 
/*+oo 


ClBEs{v) 


dv 


d 


— [x{v, u\m) + x{-v, u\m)] j-P{v) 


271 


dv 


(9.11) 

(9.12) 


and the various scattering factors are listed in Appendix C. 

Exchanging g with g we get Bethe equations for bound states of F+±. The other equations can 
be obtained from the equations written when simple gluons are present (and collected in Appendix 
E) by means of the replacements: 


Ng +0O Ng +0O 

j=l 1=1 j=l j=l 1=1 j=l 


(9.13) 


We now show that equations (9.10) and others, which constrain centers of the string (9.9), are not 
independent of equations describing excitations on the GKP vacuum, but actually can be obtained 
from these by considering strings involving gluons. 

In order to obtain (9.10), we shall hrst consider equations (2.92) for gluons. 


Na 




U 


k' 


g,9 I ,• ^F g 

' A3'K..<)n" 


j=l,^k' 


Q Q 

ul, — 




+ f o(9F)( a 

i^red \^k'^^F,j) 


N 


n 


f n,9 


+ 2 Q{9f)(,,g 


H 


=1 '^k' 2 


Sred (“fc' > “/,i ) n 


(9.14) 


h=l 


% 


n 4) n s'” ’(“i. «F,y) n «/,,), 

j=i j=i j=i 


where Sred (C.13, C.19) stand for the S factors deprived of the rational factors 
appearance of the rational factor 

Q Q \ • 

Ul, -Uj +l 

u{, -Uj - i 


_g±ixo_ xpe 

(9.15) 
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suggests that strings of gluons rapidities with the same real part and imaginary parts separated of 
i form. In specihc, 


“fc' = “fc”" + 2 ^"^ “ ^ ’ fc'= 0, ...,m - 1, 

= uf’™'+ ^(m - 1 - 2/c'), fc'= 0, ...,m - 1, 


(9.16) 

(9.17) 


where the real centers of the strings are in the region < \/2g, i«ri < n/ss- 

Performing the products over k' we arrive at the equation 


-^Vg a rn q . i / i i \ <7,m Q \ i f i \ 

1 ^ JRP^r^^ {ul'”')+2iD[^'> TT -Uj + :^{m + l) Uf. -Uj + ^{m-l) 

^ ^ II q rn q if i 1 \ 

_ —{im _L i oi 


y - t(m + 1) ur - - t(m -1) 


q(.99) ( 9,m 9^ _ 

P’red,m\^k i ^j) 


g,m I im 

“fc “ ^F,j + — ^(gF) 


N. 


f g,m 
U 


n ^k '^F,j -r 2 cigF) r g,m „ \ TT + 2 c(gf) / g,m . \ TT Q{gs)l g,m n 

g,m _ _ im^red,m\^k ’“^0/ J_J_ g,m _ _ im^red,m\^k ’ “/j 1 J_ J_ l“fc ) “/ij 


H 


j=i '' 


Up 


j=i 




/i=i 


Vg iV,, Nj 


j=l j=l j=l 


(9.18) 


where 


SlTd]m{'^,v)= n + , sl^*\u,v)= JJ S^^*\u + il,v). 


m — l 
2 


dff*) 


m—l 

2 


(9.19) 


1=-^ 


m—l 

2 


It is now immediate to recognize the scattering factors between a bound state of Ff± with center 
u and ’length ’m and a fermion or a scalar: 




u-v+^-f 


>5yL(M,'y), Si^^\u,v) = 


u — V — 
AgF) 


u 


... I im 

"-slf'L («,«), 


u — V — 


im ‘~^red,m\ 


, S^^f\u,v) = sll-2miu,v), S^^^\u,v) = u). (9.20) 

By means of (C.35) one shows that these factors equal the ones appearing in (9.10). Equations 
(9.18) are then completed by taking into account that rapidities appear into strings (9.16, 

9.17). Because of the properties 


1—1 g,m 


n 


Ut 


u 


g,i 


(/ — 1 — 2k') + |(m + 1) — u 


9,1 

j 


k'=0 '^k 


g,m 


U 


g,i 


{l-l-2k')-{{m + l) ul 


Ua 


■ds!„w’”. «?•' + -1 - 2k’)) = uf), 


1-1 


n 


fc'=0 


«r.«f + ^(' -1 - 2k')) = s!s\ur.uf‘) 


l(l-l-2k') + i(m-l) 
i(l-l-2k')-l(m-l) ' 

(9.21) 

(9.22) 


which follow from (C.35), one hnally hnds that equations (9.10) are reproduced. In analogous 
fashion, other expressions, derived from the BMN vacuum, can be reproduced, as well, starting 
from equations on the GKP vacuum. 
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The energy. 

We end this part by spending some words on the energy of the gluon and of its bound states 
(the same will apply to its barred companion). Immediately, we can put it in the general form at 
any g 

= 2m- f ^ [x{v,u\m) + $(u)] , 


dv 


(9.23) 


which we rewrite as 
E^\u) =2m + 

2tt 


dk 


ITT 


n=l 


1-00 02—6 2 


y/2ix{u + ^) 


{'yi{V2gk) - sgn{k)Y-{V2gk)) 


271 

ik 




(9.24) 


27r 


\k\ 


V2ix{u - ^) 


Jn{V2gk) - — Jo{V2gk)e 2 


By means of the relation (3.40) in [20] and of the Fourier transform 

In K n 


the expression above becomes 


^[cosfciie-‘t-l]- f 

Jo 


k 1 — e~ 
in agreement with [20]. 

Perturbative strong coupling regime. 

In this limit we can use the following results 


\/2ix{u =F ^) 
” dk 

k — 1 


(9.25) 


[cos kue ^2 — 1], (9.26) 


exp [—ix{u, v\m, /)] = exp 


\/2ml 


+ 0 ( 1 / 9 " 


and 


x{w,u\m) + ^{w) = 


m 


V^g{u - w) 

x{w,u\m) + ^{w) = 0{l/g), when |r(/| < 1. 


g{v - u) 

+ 0{l/g‘^), when |t(;| > 1 


(9.27) 

(9.28) 

(9.29) 


Repeating all the steps we did for one single gluon, we write the scattering factor between two 
bound states of m and I gluons, respectively, in the perturbative regime (up to terms 0 ( 1 / 5 ^^)) as 


Si»f>(«,9)=exp 


iml 


\/2g{u — v) 




2 V 1 -M 


1/4 


1 — V 
1 + h 


1/4 


+ - 


1 fl-u\ 


1/4 




2 \ 1 + M/ 


1 + u 
1 — h 


l/4^ 


(9.30) 

For what concerns energy (9.26) and momentum (9.11) of a bound state of m gluons, we have the 
simple relativistic result 

E^\6) = 2m cosh 6 , Pj^\6) = 2m sinh 6 , (9.31) 

in terms of the hyperbolic variable 9 dehned by m = \/2(/tanh26*. As already written, the results 
above apply to the barred gluon as well. 
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9.3 The meson and its bound states 


As hinted before, at infinite coupling in the perturbative regime a bound state of a small fermion 
and a small anti-fermion arises. As anticipated, we will call this state ’meson’. 

Now, we elaborate a proposal for the complete set of ABA equations satisfied by mesons. One 
meson is described by the stack (9.4) in which, of course, the index F is replaced by /, i.e. 


i 


Ufx — UM,k + i 

^a,k “1“ 2 

Ub,k — UM,k 

f,k = UM,k — i 

i 

"^c^k "^M^k 2 * 



(9.32) 


We already emphasised that this stack does not give rise to a genuine bound state for finite g, but 
describes a virtual state: in this perspective, we can formally write the ABA equations describing 
the scattering between virtual states (9.32) with themselves and with other excitations at generic 
g. Of course, we should start by deducing - in the usual manner - the scattering between two 
virtual mesons: 


S'(^^)(m, v) = ^ — V -\- {u + i^v + — i, n -f + i^v — — i,v — i) 


u — V — i 


u — V + i 


u — V — i 


dw 


d 


-expu / —Xm{w,u\1)—xm{w,v\1) 


27r 


dw' 




(9.33) 


with the dehnitions 


m— 1 
2 

XM{v,u\m) = E XMiv^u + ik\l) with m > 1 

k=-^ 

Xm(v,m|1) = -XH{v,u + i)-xniv.u-i) . (9.34) 


Then, we can introduce also the other excitations and derive, at least formally, the ABA equations 
involving Nm mesons 


H 


N„ 


h=l 3=1 

N-g Nf Nf ^ 

■ JJ S^^^\uM,k, U^j) JJ S^^^\uM,k,Ufj) JJ S^^^\uM,k, Ufj) JJ S^^^\uM,k, UM,j) , 




j=i 


i=i 


j=~id¥^k 


where P^^\u) = -|- i) -I- P^^\u — i), D^^\u) = -|- z) -|- D^^\u — i) and the scattering 

factors read 


! exp 


u — V — 


V, u\l) — i 


dw dQ 
2 tt dw 


v,w) xm{w,u\1) \ , 


(9.36) 
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S^^3\u,v) = 


3i 


dw 


d 




u — V + 
u — V + ^ 


(9.37) 


S^^~3\u,v) = 


Xm{w,u\1)^ [x(w,v| 1) + $(x)] - 


u -V - ^ f f dw 

-sT^^Pi* / , 

u — V — ^ y. J 271 dw 


(9.38) 


S^^-^\u,v) = 


u — V + i f f dw I \ d 

- exp <^-* / —xm{w,u\1)—xh{w,v) + 

u — V I / zvr dw 




(9.39) 


S^’^-^\u,v) = 


u — v (. I dw . d , , 

-: exp <^-* / —xm{w,u\1)—xh{w,v) + 

u — V — t I ,1 271 dw 

d^ 


+ i 


dw dz 
271 271 


CL "I 

Xm{w,u\1)-—-Q{w,z)xh{z,v) \ 
dwdz J 


(9.40) 


Eventually, these ABA equations (9.35) will acquire the aforementioned physical meaning in 
the infinite coupling limit where it is possible to parametrise um = V^g coth 29m in terms of the 
hyperbolic rapidity, and then P^^\9m) = 2sinh6'M, D^^\9m) = —Inf? and 


s^^^\eM,Uh) = 


exp 




_\/2g coth20 m - 1 


S^^-^\9M,9f) = exp 
S^^^'\9M,0f)=exp 


y/2g (coth 26m — coth 26j) 

i 1 

y/2g (coth 26m — coth 26j) 


O 


O 




(9.41) 

(9.42) 

(9.43) 

(9.44) 

(9.45) 


S^^^\6m,6'j^) = 


exp 


i^/2 


g (coth 20 m — coth 20)^) ^ \g'^ 


[S^"\6m,6'j^)]\ (9.46) 


^®When writing them, we have deliberately neglected the multiplication by {u,v) and {u,v) (which 

could be immediately obtained respectively from (9.39) and (9.40) via the usual substitution —xh{w,u) —> 
Xf{w,u) + 4>(ry)), since we know (c/. (4.50) and (4.47)) that this contribution from the large fermions, forced 
to live in the giant hole regime at strong coupling, would be exponentially small. 
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Explicit expressions for the scattering S^'^^\9M,Uh) is given by (4.71), with v = coth20M and for 
all the other S factors (entering the previous r.h.ss.) in Appendix C.3. 

Now, if we come back to (9.33), we can notice that, in analogy to equations (9.1) for small 
fermion/antifermion (and to (9.14) for gluons), a rational pre-factor in front of the scattering 
factor is present. This suggests that the same mechanism, which produces (virtual) mesons and 
bound states of gluons should also produce (virtual) bound states of mesons: in fact, these will be 
represented by strings with the customary form 

UM,k = u^ + ^{m-l-2k), k = 0,...,m-l, (9.47) 

where all roots, once scaled, collapse towards the real axis at inhnite coupling. At generic coupling, 
these virtual states enjoy scattering phases which, formally, do not distance themselves from the 
ones previously found for actual particles; indeed, when considering a process involving a bound 
state of I mesons with center u and a bound state of m mesons with center v, the resulting scattering 
matrix emerges from the fusion of (9.33): 


q(MM) 


{u,v) 


1—1 m—1 


n n ■5'“"’ ('“ + h' - 1 - 2j ).« + hm -l-2k)) = 


j=0 k=0 




/■+°° dw d 

exp\-ix{u,v\l,m) + i / —XM{w,u\l)—XM{w,v\m) 






(9.48) 


with the definitions (9.34). Similarly, from a generalisation of (9.36)-(9.38) we can deduce the 
scattering phases concerning bound states of m mesons against scalars or bound states of I gluons: 


c{Ms)i 


u-v + ^^ 
u-v-^f 


exp 


. . f dw dQ , ^ / I X 

v,u\m)-i / — — {v,w)XM{w,u\m] 


(9.49) 
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u — V -\r ~\~ 2k “1“ 2) f 

-=--- exp< i 

M-u-f(/ + 2fc-2) 


dw 


XM{w,u\m) 


d 

dw 


[x{w,v\l) + <h(x)] - 


/ dw dz 
2tt 2tt 


XM{w,u\m) 


d^ 

dwdz 


Q{w,z) [x{.z,v\l) + <h( 2 ;)]|. 


(9.50) 
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XM{w,u\m) 


d 

dw 


[x{w,v\l) + 4>(x)] - 


/ dw dz 
27r 27r 


XM{w,u\m) 


d^ 

dwdz 


e{w,z) [x{z,v\l) + <h( 2 ;)]|. 


(9.51) 
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It is crucial for the following (TBA) noticing the drastic simplihcation of the fusion into a multi¬ 
plication when g = oo (zero shifts) 

g{MM) ^ ^ ^g_52) 

along with the scattering phases between a bound state of I mesons and a bound state of m gluons 
which become 


qiMg} 

^Im 


(0, e') 



e') = e)]-^ = e')f^ 


(9.53) 

Explicit expressions for (9.52, 9.53) in terms of hyperbolic variables are given by formulae (C.45, 
C.46) in Appendix C. 

We need to conclude this section with the dispersion relations. At least in principle, we can 
write the energy and momentum of the virtual meson (of mass = 2 ) and of its m bound state at 
hnite g, albeit these are not stable particle: the energy as 



2m + 7^)(M) 


(9.54) 


with 



dk 

cos ku 

k 

sinh 1 

dk 

cos ku 

k 

sinh 1 

dk 

sin ku 

k 

sinh 1 

dk 

sin ku 

k 

sinh 1 


^coshfc e 2 ^ — cosh 7 ^(\/ 25 f/c) for m odd, 
(cosh fc — l) 7 ® (v^^ffc) for m even; 


^cosh ^ — cosh k e 2 Y-i.'^gk) 

(1 — coshfc 'yti'/^gk) for m 


for m odd, 
even, 


(9.55) 


(9.56) 


both in terms of the functions 7 ^ dehned in (B.19, B.20). The (rest) masses of these resonances 
are protected at all values of the coupling to be the purely additive values Eln\u = ±cxo) = 2m 
(zero binding energy), which are indeed stable as long as the momentum is zero. Yet, these acquire 
full (any momentum) stability in the strong coupling limit where boost invariance is recovered. 
In fact, if we enter the perturbative regime g -|-cxo, u = V2gu, with |h| > 1, we obtain the 
relativistic dispersion relations 


E^\6) = 2m cosh 6 ^, = 2msinh0, (9.57) 


in terms of the hyperbolic rapidity 0 , dehned via u = coth 26. 
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10 Pentagonal amplitudes at strong coupling (perturba¬ 
tive regime) and confinement 


An important application of the above scattering data, which implies a non-trivial check of them, 
is the construction of the so-called pentagonal amplitudes, P [34, 35, 36, 37, 38]. The latter, in 
their turn, are the building blocks of an inhnite expansion - the BSV series - of the gluonic (MHV) 
scattering amplitudes. In this section, we want to compute the pentagonal factors, P, relevant at 
large so to prepare the analysis of the BSV series (at strong coupling) in next section. 

The BSV series is a sum over the (intermediate) multi-particle states, where the particles may 
be, - at generic finite coupling -, scalars, fermions, gluons and bound states thereof, as analysed 
above. The simplest example is provided by the six-particle amplitude (or, in other terms, the 
equivalent hexagonal Wilson loop) 


H-CX) 


1 


N r 


Whex = 1H 

N=0 ai 


ajv 


k=l 


dUk 


paAuk)e ^ 


-TEai^{uk)+ia-paf.iuk)+imaf.(l} 


X 


X R 




{0\Ui...UN)Pai...aN{-UN--- “ Ml|0) , 


( 10 . 1 ) 


which is expressed by means of the measures ^a^Ui) (corresponding to quadrangular amplitudes) 
and the multi-particle pentagonal amplitudes Pai...aN{0\ui...UN), representing the transition from 
the vacuum to an intermediate state with N particles of the kinds listed above, each one associated 
to a label a*. When we go to the strong coupling limit, we have to disentangle the integrations 
over internal rapidities by performing the limit g —?■ -|-cxd in the integrand. This procedure means 
that we have to add different contributions. 

The hrst one comes from performing the limit g —?■ -|-cxd with integration variables fixed. This 
part depends on excitations in the non-perturbative regime and is dominated by scalars, and may 
reserve very interesting surprises as anticipated in [37]. In fact, this contribution would come from 
a (genuinely) quantised string in and would elude the minimal area argument of the AdS^ 
string. However, this regime misses contributions from regions in which rapidities are large: these 
are recovered by adding the integrals in which integration rapidities are scaled before taking the 
limit g —)■ -|-oo. 

More precisely, if we scale the integration rapidity u = \/2gu, with u hxed, we have the two 
following regimes. If |h| > 1 we are in the giant hole regime. In this regime all the excitations 
behave in the same way. In particular, as we showed in Section 4, scattering phases —i In S between 
any pair of excitations are all the same and are all proportional to the coupling g (4.3). The same 
happens to energies and momenta [20]. This property is crucial, since it implies that contributions 
to scattering amplitudes coming from integrations in these regions (all scaled rapidities Ihl > 1) 
are exponentially suppressed. 

Instead, things are different if |h| < 1, i.e. in the (string) perturbative regime (for all particles 
except scalars: the rapidity of the latter in this regime does not scale, but instead u = w + y lnm( 5 f), 
where m{g) ~ ^ as seen above). In this regime, energy, momentum and scattering factors 

are expanded in inverse powers of the coupling constant g. Additional structure is added when 
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expressing the pentagonal transition between a M particle state to a iV particle state in terms 
of the one-particle-to-one-particle transitions because of the matrix representation carried by the 
single particle (thus the singlets makes an exception to this). In this operation polynomials in the 
rapidity appear as denominators, taking into account the different representations to which the S 
matrices can belong. For instance, in the case of the hexagon (10.1), extensively discussed below, as 
we start from the GKP vacuum, we need consider only pentagonal amplitudes to the other possible 
singlet states. In particular, this polynomial is a monomial in the case of the transition (from the 
vacuum) to a two particle state of a fermion and an antifermion (which, though, belong to the 4 
and 4, respectively). This monomial ’squares’ in the integrand of the amplitude contribution to 
( 10 . 1 ) 


—mIO) 


1 1 
{u — vY + 4 PPP[u\v)PPP{v\u) 


( 10 . 2 ) 


Instead, the transition from the vacuum into the two scalar singlet is even more depressed, albeit 
the rapidity does not scale (for a scalar, but is added a (^-depending constant). In fact, the P factor 
contains at the denominator a polynomial of degree 2 multiplied by {cf. [35, 36] for details), 
and then the two scalar contribution to the hexagonal amplitude writes down: 


= 3 


f dudv 


g—T[i?® (n)+i?®(i;)]+i(T[p® (ul+p® (ti)] 

1 (27r)2 gY{u 

— vfi + 4][(m — vfi + 1] 

PW)(M|y)pM(p|M) 


(10.3) 


where we ought to consider that = 0{g) and P*'^*^(M|f) = 0{l/g). Hence, this integral turns 

out to be of order wjffj = 0{g^), then subdominant with respects to semi-classical approxima¬ 
tion (contributed by the gluons, for instance). Actually, while the scalar contributions are really 
subdominant in the perturbative regime, on the contrary fermion ones behave in a subtle manner: 
in fact, the lorentzian function in front of (10.2) would entail a contribution from the singularity 
u — V = P^plijg pinching the real axis when g -|-cx) [36]. But in our picture this is is the con¬ 
tribution given by their bound state, the meson indeed. Moreover, also the greater multi-fermion 
coalescence are taken into account by the multi-meson and meson-bound-state contributions, c/. 
below. In summary, we are in the presence of a phenomenon in which the fermions coalesces at 
least in a fermion and antifermion couple and disappear from the spectrum as free particles: as 
anticipated, this is a sort of confinement typical of MHV gluon scattering amplitudes/Wilson loops 
at strong coupling, not evident at hrst glance from the 2D scattering factors. To be fully precise, 
although the string theory minimal surface conhrms this disappearance, and the appearance of the 
meson [50, 6] and its bound states [39, 41, 40], nevertheless a detailed multi-fermion description is 
missing so far [59]. 

Concluding this preamble, these polynomials in the denominator produce in the general case 
negative powers of the coupling constant after scaling^^ the rapidities and thus ’depress’ the ampli¬ 
tude. Of course, these polynomials are absent if the excitations belong to a S'f/(4) singlet (see also 
[35, 36] for a detailed analysis of the two particle case). Therefore, we can argue that the leading 
contributions in the perturbative regime are due to particles behaving as singlets under SU{4:), 
indeed. They are gluons and their bound states, as already proven by the detailed two particle 

^^Rapidity of scalars do not need to scale. 
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analysis of [36]. But, at strong coupling, we have shown necessary to add mesons and their bound 
states to the spectrum, as well. 

Now, the pentagonal amplitudes P enjoy at general coupling a series of axioms depending 
on the S-matrix entries. Therefore for the latter we need to use our previous (strong coupling) 
perturbative expansions at leading order and ’solve’ the axioms. For exposition’s sake, we give 
in the following the complete list of the P factors (gluon-gluon, gluon-meson, meson-meson and 
bound states, contributing at leading order), leaving the details of their derivation in the Appendix 

We start from the gluon and then the bound states of ^ of them. In this gluonic sector the 
rapidity enjoys (at perturbative strong coupling) the parametrisation u = \/2g iaii\i29. Thanks 
to this, the three axioms (6-8) in [34] for the gluon {g) and its barred companion {g, the other 
component of the massless spin 1 held) simplify their arguments: 

pi99)[-e\ - e') = p^33\e'\6), p^33\-e\ - o') = p^^~^\e'\e) , 

p(99)(0|0') = A(®^^)(0,0')p(99)(0'|0), p(^'5)(0|p) = S^33)^Q^Q'^p{99)^Q>\Q') ^ 

p(99)(0 - m/2\e') = P^33\d'\6 ), 

and we can solve them with input the leading order expansion of the gluon-gluon scattering matrix 
(C.40, C.41). We obtain 

= 1 + 2^20 ‘sinhpg - 290 “ ^')1 + . (lO ^) 

and the symmetric channels P^^^\9,9') = P^^^\9,9'), P^^^'>{9,9') = P^^^'>{9,9'). The constant a 
may equal ±1: its precise value is not hxed by the axioms, but by the comparison with data derived 
from the Thermodynamic Bubble Ansatz (TBA) in [41]. As we wrote above, formulae (10.7, 10.8) 
with a = 1 have been already reported in [34]. 

For what concerns P factors of gluon bound states, we may conjecture, along the lines of 
the previous equations (10.6) for the single gluons, the following functional relations as axioms: 
= PS\D,ff), P‘S\e,ff) = PiSHe,ff) and moreover 

-«') = p‘if(e',e ), - irr/2,«') = Ptf («',«), 

Pt\o,ff) = s'S\e,ff)PiS\e',e ), p/f>(»,»') = s[f{,e,ff)p‘S\ff, 0 ). ( 10 . 9 ) 

Moreover, we recall that the F-matrix factors are simply multiplicative at perturbative strong 
coupling: S^l'\9,9') = [S^^^^(9,9')]"^^ and S^P(9,9') = [S^^^^(9,9')]"^^. Therefore, solutions to 
(10.9) should enjoy the same property, which entails upon expansion for large g 

amiP!pP(9,9') = 1 + Xh^2^^-2P) ^ ~ -ismh{9 - 9')] + 0(1//), (10.10) 

have to say that the expansion of the gluon-gluon P factors - formulae (10.7,10.8) - previously appeared in 

[ 34 ]. 


(10.4) 

(10.5) 

( 10 . 6 ) 
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(^miPl^i\0,e') = 1 - 2^0 + *sinh(0 - 6')] + 0(1//), (10.11) 

or in barred rapidities u = tanh 26 


(^mlPl^l\u,u') = 1 + 


iml 


2V2g 
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( 10 . 12 ) 


amiPl^f{u,u') = 1 
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2y/2gu - u' 

l + uV/^ 

1 — u) VI + m' 


;(1 + ^) 


l + u' 
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1/4 


1 — u 
l + u 


1 / 4 . 


(10.13) 


Overall constants ami = cizm can be equal to ±1 and are constrained by the comparison with the 
TBA of [41]. 

Let us now consider the meson and its bound states, and in particular recall that for all of 
them the rapidity enjoys the perturbative parametrisation u = \/2(/coth26'. As we discussed 
before, these are self-conjugate particles and this property allows us to postulate the following 
set of functional relations (which now will be meaningful only in the perturbative strong coupling 
regime, where the particle does exist) for the single meson P factor: 


6 ') = ^ 

p^^^\e,e') = s^^^\e,e')p^^^\e',o) , 

p{MM)i^Q _ ^^^ 2 , 6') = 6) , (10.14) 


where 5'A^) is given by (C.45). We write the solution of (10.14) as 

1 i sinh 26 sinh 26' 


/3P(^^)(0,0') = 1 


^/2g sinh(26' — 26' 


y\/2cosh y-P+0(l/^2), (10.15) 


where (3 = +1. For mesons bound states, we have anew the multiplicativity of the scattering 
factors in the perturbative regime, namely Sl^^\6,6') = 6)')]^^ Which, in its turn, 

imply the same property on P factors, i.e. upon expanding at large g 


AW («.«') = !- 


ml i sinh 26^ sinh 2P 
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or in barred variables u = coth 26 


/2cosh +0(1//), (10.16) 
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where Pmi = Am can equal ±1. 

Eventually^^, we consider the scattering between (bound states of) mesons and (bound states 
of) gluons. We are now looking for functions Pm^\ which may conjecturally 

satisfy the functional properties (meaningful only at perturbative strong coupling) 
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(10.19) 


We write solutions to these equations in the form P^P = 1 + + O (y), where 

, , . (Mg) (gM) (Mg) (gM) , 

can equal ±1, fmi = Tim = Tmi = Tim and 


A'W9)(0,0') = K^^~9){9,9') 
K(9M)^0>^0^ = K^~9^\9',9) 


sinh 29 cosh 29' 
V2 cosh(20 - 29') 
sinh 29 cosh 29' 
y/2 cosh(20 - 29') 


[sinh(6' — 9') +i cosh(6' — 9 ')], (10.20) 
[sinh(6* — 9') — i cosh(6' — 9 ')], (10.21) 


or alternatively in the barred variables 
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( 10 . 22 ) 

\ 1/4. 


w 


(10.23) 


^^Even if contributions of small fermions to amplitudes is suppressed at strong coupling with respect to gluons 
and mesons, we give also the strong coupling limit of their P factors. We refer to formulae (38) of [36] and use 
formulae (C.42, C.49, C.50) for the (strong coupling) perturbative regime of the fermion-(anti)fermion scattering 
factor and its mirror, respectively. We eventually obtain 


[p^P\e,e')f 
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11 Hexagon at strong coupling 

11.1 Aim and assumptions 

Now, we want to compute the hexagonal Wilson loop as the series proposed by [34, 35, 36] 
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( 11 . 1 ) 


at strong coupling. As argued above, in this regime intermediate states which contribute are gluons 
and their bound states, together with mesons and their bound states. All of them are singlets and 
then for their pentagonal amplitudes a simple product and inversions hold when changing a rapidity 
from in to out: 
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( 11 , 2 ) 


This formula entails an easy product to appear inside the hexagonal amplitude: 
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(11.3) 

where the indices label the species of different particles (including bound states): this is the 
formula we want hrst to match with initially, and then to sum up. 

For the gluon and the bound states of i of them, rapidity may be parametrised as n = 
\/25ftanh26*. Then, we can recall their relativistic energy and momentum (9.31) and notice how 
they are purely additive 


Ef{u) = \/2£cosh0 = iEf{u ), Pi{u) = V^f'sinh^ = ip{{u) 
In the following, the expression will be written in terms of 

yx 


(11.4) 


27r 


= V2g 


in order to facilitate the comparison with the existing literature on this subject. Gluonic measures 
appearing in (11.3) is given by (£ = 1 is the gluon) 


= —/Arm =_ - _— 

277 ’ 


(11.5) 


with given by (10.10). In order to have agreement with TEA it is enough to choose = 
(—1)^“^. Then, in terms of the single gluon measure (at strong coupling) 


p^{u) = -1 + ... 


( 11 . 6 ) 
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(11.7) 


we can express the i gluon bound-state ones 

du Q. . \/A , do 

— uAu) = —(—1) -IT— 

^ 27r^ ’TiPcoAA2e 
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27r P ^ ' 271 £2 

On the other hand, let us remind that the meson and its bound-states enjoy the rapidity parametri- 
sation u = \/25fcoth26' with relativistic energy and momentum (9.57) which are purely additive 
(zero binding energy) as well 

E^{u) = 2m cosh 0 = mE^{u ), Pmi^) = 2msmh9 = mp^{u). 


The measure for bound states of m mesons is 
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( 11 , 8 ) 


(11.9) 


e'^g 


with given by (10.16). We choose (3mm = (—1)™ ^ in order to have agreement with TBA. 

Then, similarly to gluons, we obtain the meson bound-state measures 
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( 11 . 10 ) 


where the single meson measure (at strong coupling) 


p^{u) = -1 + 


( 11 . 11 ) 


We shall point out that, for compactness sake, the paths of integration are not explicitly written 
in what follows, so we need to describe now the choices we adopt throughout the rest of the section. 
The integrations over gluon rapidity u can be computed along the interval which 

turns to the whole real axis when the map to hyperbolic rapidities 9 is performed u = ^ tanh(26*). 
Otherwise, the integrations for meson u-rapidity is intended on a section of a straight line below 
the real axis, lying in the small fermion sheet [41] [36] (we remind that large fermions in this 
regime behave as giant holes, hence their contributions are exponentially suppressed because of 
their dispersion laws [20] and scattering factors (4.50) ), namely (-|-oo — ie, +^) U (—^, —oo — ie), 
where 0 < e < 2. When switching to meson 9 rapidity according to the map u = ^ coth(26'), the 
integration may be set to run along the straight line (—oo -|- ip, -1-cxd -|- ip), where p (according 
to [41]) can be hxed in terms of conformal ratios, as shown below (11.32). Eventually, the same 
straight line is chosen for the integrations on the gluonic 9 rapidity. 


11.2 One particle 


Let us start from one particle contribution. With ’one particle contribution’ we mean that in 
(11.3) we consider only one insertion, which can be a gluon, a meson or bound states of thereof. 
Therefore, we can write down 
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which at strong coupling reads 
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where Li 2 (z) = X] + is the dilogarithm function. 


m=l 


11.3 Two particles 

Let us pass to two particle terms = wjfj^ + in which we distinguish three 

contributions: gluon-gluon, meson-meson and gluon-meson. 

Gluon-gluon 

Let us start from gluon-gluon: 
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At strong coupling the symmetric product of P factors (10.7, 10.8), entering (11.14), enjoys the 
property 
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We can now perform the sums over £i, £2 and get 
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A'S(9i.92)i‘'(9.)i'’(«2), (11.19) 


where we used the short notation 


L®(0) = In 


_j_ \/ 2 t cosh 6+i\/2cr sinh 6 j | _j_ g—* 0 —\/ 2 t cosh (9+i\/2a'sinh 6 


( 11 . 20 ) 


Meson-meson 

The meson-meson contribution is written in a completely analogous way: 
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Expression (10.16) for the strong coupling limit of mesonic P factor implies the property 
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We can now perform the sums over mi, m 2 and get 
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where we used the short notation 
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, (11.24) 
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Next step is to consider the meson-gluon contribution 
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Now, at strong coupling, with the redefinitions ui = \/ 25 fcoth 26 'i, U 2 = \/ 25 'tanh 202 expressions 
( 10 . 20 , 10 . 21 ) imply the property 


where 


97r 

= l-^(-m/fi"»>(9i.92). (11.27) 


cosh 26*2 sinh 26*i sinh( 6 *i — 6 * 2 ) 


= a'<"»>(»i.«2) + a(»"'(»2,»i) = Ai;")(»2,«i) = V2 

(11.28) 

Remembering the measures (11.7, 11.10) and the forms of energies and momenta (11.4, 11.8) and 
performing the sums, we arrive at the expression 
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11.4 Comparison and checks with TBA 

We now compare our previous predictions at strong coupling with TBA outcome. We use (F.42- 
F.46) of [41]: these expressions depend on the functions e( 6 * — i(p), e{9 — i(p), which satisfy integral 
equations 

.(9-2^) ^ ^l(9)-/-^A<»)(9.9W) + /^^AgS>(9^9)L(^ 

(11.30) 

7(9-2^) ^ 72^(9)-/-^A<*S>(9.9W) + /-^A.-^ 

where the integrations are performed along the straight line R -|- iip, and the functions L{6), L{9) 
are defined as 


L{e) = In (1 + , L{e) = In [(l + (l + 


while [41] 
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The relation between and the cross-ratios can now be clarihed by means of the functions e, e, 
since 


In 


M 2 M 3 


= e(f7r/4 — i(p) , 


In 


U2 


1 = e{-i^) . 


We expand (11.30) at large E{6) and plug such expansions in (F.45, F.46) of [41], by using the 
Taylor expansion (|/| 1) 


Li 2 (—e ^ — Li 2 (—e / ln(l-f-e ^) + 0{p). 
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where we used the notations 
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Summing this result with (F.42), (F.43), (F.44) we get 
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AiJS'(9i,92)rE(9.)iE(92), (11.34) 


(11.35) 


(11.36) 


Now, in order to compare with our results, we parametrise the cross-ratios ui, ua, W 3 as in [35] 
(formula (157)), i.e. 


— = 1 + , — = 1 + (e-" + (e"" -F , 
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Consequently, if we define fi = we obtain the relation 


ji -\- jj, 


1 — Ml — M2 — Ms 


(11.38) 


In addition, the function E{9) may be written in a more transparent form as 

£’(0) = \/2r cosh6* — iV^cTsinh^ . (11.39) 

Plugging (11.38, 11.39) into (11.36), the following relation 

in.. = = exp ^ 

\ fe=42 / fc=42 V / fc=42 

(11.40) 

is in agreement with expressions for ^hex computed in last subsection. 


11.5 Re-summation of the BSV series 


The agreement displayed above between the series written in [34] for hexagonal Wilson loops and 
the TEA for scattering amplitudes [39, 41, 40] can be made even tighter, since it is not restricted 
to one and two particle contributions, but instead it does also extend to any number of particles. 
Even better, the BSV series for the hexagon (11.3) can be fully re-summed by exploiting some 
standard techniques: eventually we will reproduce (as for the strong coupling regime) the TEA 
(in the form elaborated in) [41]. In the following we will produce the main steps, but leave some 
further details and generalisations for an incoming publication [59]. 

The expression to sum up is the simple manipulation of the initial formula, (11.3), strictly 
speaking valid for singlets, which we re-call here for practical reasons 


+00 


N r 


A^=0 ai fliv k=l 


duk 




-I N 

n 

Kj 


■ Pai,ajiUi\Uj) Paj^aiiUj\Ui) 

(11.41) 


where the indices a^, label the species of different particles (including bound states). 


In general, not only at strong coupling, we may use a path integral trick of the type as in [61, 62], 
but then we should integrate eventually the extra p-£eld(s) [59]. Thus, we better perform on the 
above series (11.41) a similar trick without the p £eld(s), the Hubbard-Stratonovich transform [64]. 
The latter makes use of the well know identity of functional gaussian integration in the presence 
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of a linear source term 


N 

JJ" g{^(ai)K)>^(aj.)K')) ^ ^g^(ai)(«l) . . . g^(a^) («iV) ^ ^ 

i<j 


(11.42) 


(^g{^K)K)^K)K)> ~ fQr allowing the summation to act on the single exponential of the r.h.s.^^. 
This means that we need also to relate the pentagonal amplitudes Pa,b{u\v) to correlators (and 
then to the kinetic part) of the gaussian held X(^a) in this way 


_^_ = g(^(a)W^(6)h)> _ 

Pa,b{u\v) Pb,a{v\u) 


(11.43) 


For instance, we associate gluons to the (huctuating) held = X^{u), whereas we denote their 

bound states as X^^^{u), where i stands for the number of components. Thus, the ’linearisation’ 
of the exponent is complete, namely we can recast the hexagonal Wilson loop (11.41) into a shape 
aiming at re-summing the diherent contributions: 


Whex = {exp 


du 


Q--rEa{u)+iapa{u)+ima<t> I ^ ^ 


(11.44) 


With this series (on any particle kind, a) inside, (11.44) can be interpreted as a Kac-Feynman 
path integral (partition function) for any value of the coupling, nevertheless it is at strong coupling 
A —)■ oo that a tangible simplihcation occurs since it can be summed up. For simplicity’s sake, 
we will initially include only gluons and their bound states, then adapt our derivation easily to 
mesons and their bound states. Eventually, we will consider the general system (at strong coupling 
only), composed of gluons, meson and bound states. For, in this regime, the bound states enjoy a 
series of simple, useful properties: their energies and momenta (11.4) are simply additive, as so is 
the relation X^^<^ = aX^^^ implied by the peculiar limit form iii 2 Kiym{u,v) of the (bound state) 
gluonic kernels in (11.15) via (11.43) on the bound state helds X^^y, and hnally the measures 
(11.7) exhibit a peculiar square at denominator. Altogether these properties turn out to be crucial 
to re-sum the gluonic part of the hexagonal Wilson loop (11.44), wjf2;, in a handy shape, and they 
bring up the dilogarithm function Li 2 (a:) (tuned by the third property, cf. [59] for more details): 


= (exp <1 - / — ). 


27r 


(11.45) 


23 


The following formula is the infinite dimensional d ^ oo version of 


^gXlJlgX2J2 . . . gXdJd 




dX, X X.TyXj Y:XiJi k X 


= 1 — p, i,j = l 


with propagator G = where we choose, for the continuum limit of the external field Ji —>• J{u), the configuration 
of point-like sources J{u) = X]i=i — Ui), cf. [59] for details. 

^^We mention the talk held by B. Basso at IGST 2013 in Utrecht concerning one gluon (without bound states) 
and [63] for useful suggestions. We wish also to notice the possibility of interpreting the free boson c = 1 2D CFT 
(Coulomb gas) correlation function formulm by means of this one. 
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Now, we can make explicit the gaussian measure in (11.45) as a kinetic term so to read W^ll. as a 
quantum mechanics partition function for the held X^{u) 


= [ VX^ , 


(11.46) 


where the action directly expressed in terms of the hyperbolic rapidity 6, has the form 


5 ( 9 ) [X5] = - j dode'x\e)T^{e,e')x3{e') + 

+ [^ /i®(0') [Li2(-e-^(^')+*'^e^'’(^')) + Li2(-e-®(^')-'^e^^(^')) 


(11.47) 


with E{6) coinciding with the derived (11.39). Of course the kinetic kernel T^{6', 9") is the inverse 


j do' G^{e, e')T3{e', 9") = s(9 - 9") 


(11.48) 


of the Green function 

G^(9,9') = (X^(9)X^9')) 25. (11.49) 

Remarkably, the action (11.47) is proportional to g, which is going to +oo, so making possibile 
the applicability of the saddle point with classical equation of motion: 


X^9)- f —G^(9,9')/u^(9') log (l + e^®(^')e-®(^')+*'^)(l + e 
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XS(e’)^-E(9')-i 


) =0, (11.50) 


where the Green function, at strong coupling, can be easily related to (the symmetric part of) the 
gluonic pentagonal amplitude 

Oyr 

G^e ,«') = -^ 9') + 0(1/A) . (11.51) 

The introduction of the ’pseudo-energy’ e{9) via the relation e{9 — iip) = E{9) — X^{9), leads 
us to the special version of the TEA equations for gluons only, the hrst of (11.30) in which we 
fully neglect e{9 — i(p), i.e the meson contribution. In other words we have found an action (a 
Yang-Yang functional) whose differentiation give rise to equations in TEA form [29, 28], without 
thermodynamics. In this respect the generation of the Li 2 (x) function via summation on bound 
states is of fundamental importance. 

As for the mesonic sector, the reasonings outlined above can be easily repeated by substituting 
the gluon and bound states thereof with the meson and bound states thereof, respectively 2®. In 

could realise that directly G^{9,9') (instead of T^{9',9")) appears in the action by means of an Hubbard- 
Stratonovich transform which introduces some gaussian field p® coupled to X^\ in this way we will end up with 
the usual form of the Yang-Yang potential of Nekrasov-Shatashvili [60] as it would be following ah initio the path 
integral trick contained in [61, 62], c/. [59]. 

^®As anticipated about the bound state analysis in section 9, a more mathematical understanding of the contri¬ 
butions of the mesons, as small fermion-antifermion state, and their bound states should be given in future [59] 
with a mechanism where the poles pinch the integration axis [60, 61, 62, 63, 36]. 
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first place, we associate the fields {6) to bound states of mesons, each represented by the single 
meson X^s^{6) = X^{6) by means of the relation X^-^{6) = (.X^s^{6). From the identihcation 


pM, 




= e 




v\u) 


(11.52) 


it follows that the meson-only hexagonal Wilson loop assumes a shape analogous to (11.44) 

and can be re-summed at all coupling, even though a remarkable simplihcation occur at strong 
coupling, owing to the properties of the mesonic kernel: 


<“> = 


exp 




(11.53) 


Again, the meson hexagonal Wilson loop can be associated to a partition function, dehned via the 
action [X^] 

= l j dO'd 9 ” X^{ 0 ')T^{ 0 ', 0 ")X^{e") + j ^p^(0")Li2(-e-^^(®")e^""(®")) 

(11.54) 


which, under extremisation, gives the equation of motion: 


x^{e) 


dO' 


G^{9,e')p^{0') log 


1 + ^x’^{e')^-V2E(B') 


= 0 


(11.55) 


where the mesonic Green function has been introduced 


j d 9 ' G^{ 9 , 9 ')T^{ 9 ', 9 ") = 6{9 - 9 ”) . 


(11.56) 


If we dehne the function e{9 — iip) = \/2E{9) — X^{9), we obtain uniquely the mesonic TEA 
equation, i.e. the second of (11.30) where we discard the gluonic contribution (considered above). 


Complete system 

After the considerations outlined above for incomplete systems, made of a single type of particle 
(and relative bound states) at one time, we can now cope with the complete system, including 
gluons and mesons together, by arranging the gluonic and mesonic helds into a vector, and the 
measures as well: 

^ . ( 

(the label a takes the values a = 1,2; the sum convention on repeated 
complete hexagonal amplitude can thus be expressed as 

Wne. = j VX, VX2 , (11.58) 


^ (11 57) 

p^{9) J ^ 

indices is assumed). The 
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where the action reads 


5[X] = - / d9d9'X%9)Ta,{9,9')X\9’) + 


(11.59) 




/ii(6>) Li2(-e 




+ /^i(^) Li2(—e 




fi2{9) Li2(-e 


-V2E{e) gX2(0)^ 


and the Green functions, now arranged in a matrix, can anew be associated to the pentagonal 
amplitudes 


G‘^\9,9') 


K%t{9,9') Ki^ri\9\9)\ 

^/\\k'^^^\ 9,9') Ki^r^\9,9') ) 


(11.60) 


and allow us to dehne the kinetic 2x2 matrix as its inverse: 


j d9' G‘^\9,9') Th,(9', 9") = 5“ 6(9 - 9"). 


(11.61) 


Upon extremising the action iS[X], we obtain 

X‘^(9) - G“'(0, 9') fi,{9') L\9') = 0 , (11.62) 

where we have assumed the definitions 


L^(9) = log [(1 + ^ ^xHe)^-E{e)-i 4 >-A ^ ^ j^g 


1 + e 


X 


■2(0)g-v/2i?(0) 


(11.63) 

These equations of motions (11.62) match the TB(ubble)A equations (11.30) of [41], provided we 
identify the pseudo-energies as e(9 — i(p) = E(9) — X^{9) and e{9 — i(p) = ^/2E{9) — X^{9), which 
simply translates into 


L(9) = L\9), L{9) = L\9). 


(11.64) 


Since the action (11.59) possesses the diverging prefactor iS[X] oc \/A, the hexagonal Wilson loop 
Whex (11.58) is dominated by the classical conhguration, achieved by imposing the equations of 
motion on the fields, and therefore with the aid of (11.62) we can rewrite the kinetic term in the 
action (11.59) as 


^ j d9d9'X\9)Ta,(9,9')X\9') = 

y/X f d9d9' L\9)L\9') ^^X f d9d9' L\9)L‘^{9') 

2e j (27r)2 cosh26' cosh20' cosh(0 — 0') 2e J 2e‘^ sinh 20 sinh 20' cosh(0 — 0') 
^ f d9d9' L\9) L\9') sinh(0 - 0') 

27r J TT^ sinh 20 cosh 20'cosh(20 — 20') 


Eventually, the sum of the kinetic term (11.65) and the potential part, given by the second line of 
(11.59) computed on the solution on the saddle point equations, amounts to the critical Yang-Yang 
functional ^ YY^r , which has been computed in [41] by adding together the right hand sides of 
the formulae from (F.42) to (F.46). 
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12 Conclusions in perspective 


We have derived the complete set of Asymptotic Bethe Ansatz (ABA) equations referring to the 
GKP vacuum instead that to the half-BPS state (Beisert-Staudacher equations [2]). These describe 
the dynamics of all the elementary excitations over the GKP vacuum (gluons, fermions and scalars), 
but they also admit solutions in the form of complexes of Bethe and/or auxiliary roots, the so- 
called strings or stacks. The latter are the bound states, among whose the most important are the 
bound states of the elementary particles (the other are bound states of the auxiliary or isotopic 
roots, yet important for the spectrum TBA and so on). In this way, we have performed the ’fusion’ 
of the fundamental (elementary and isotopic) excitations, which is in its whole an alternative way 
to perform the bootstrap of A-matrices (c/. for instance [65] for a review). 

Moreover, we outlined this system of algebraic equations at all coupling values, also including 
weak and strong (in different dynamical regimes) coupling. Above all, we have mainly focused on 
the scattering phases between all kind of particles at any coupling, but also the new feature of two 
defects has arisen in the form of new scattering phases for any flavour. Then, we have devoted a 
meticulous care to the behaviour of the scattering factors in the three possible, - non-perturbative, 
perturbative and giant hole -, regimes which allow different large g expansions: for all these three, 
we obtained explicit expressions of all the scattering factors. 

If the momentum of any particle enters the ABA equations, the energy/anomalous dimension 
is the hnal object expressed via a solution of these equations. And we could conhrm for these first 
two conserved charges the achievements by [20], but also have been led to consider all the higher 
integrals of motion (which do play a so important role in the construction of the dressing factor in 
the usual ABA on the BMN vacuum). 

A deeper look at the form of these new ABA equations brought to our attention an interesting 
property or identification for them: the sm(4) residual R-symmetry constraint the elementary 
particles to have as rapidities the inhomogeneities of a sm(4) symmetric spin chain of S'-matrices 
which belong at any lattice site to the characteristic representation of the particle, i.e. 1, 4, 4, 6 
(for gluons, fermions, antifermions, scalars, respectively). Thus, as anticipated in [23], the matrix 
structure of the ABA equations could be inferred from the SU (4) symmetry, but the specific form 
of the scalar factors and its gf-dependence must be computed explicitly. For instance, in this 
perspective, the two defects are simply two purely transmitting impurities which still respect the 
SU{4:) symmetry. Moreover, the particular ^f-dependence shows explicitly the decoupling of the 
six scalars in the non-perturbative regime and their approach to the 0(6) non-linear sigma model 
A-matrix in [53], being, besides, the defects of no importance in this limit. More importantly, we 
have seen from the fusion of a fermion and an antifermion the formation of a new particle in the 
g —)■ -|-oo perturbative regime: a meson. Then we also identified bound states thereof. 

At last, but not least we have been looking for confirmation and deep comparison of our careful 
strong coupling outcomes with the scattering amplitude/WL TBA [39]-[40] via the OPE or flux 
tube (BSV) series [41], [34]-[38] (see, also, e.g. [67] for the multi-Regge limit). In fact, the basic 
object of the latter, the so-called pentagon amplitude, can be expressed via the aforementioned 
scattering factors as proposed for the gluons in [34]. The bound states of the latter, the meson 
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and its bound states appear to be the only other relevant particles at leading order (the minimal 
area of classical string). Therefore, we have checked explicitly those features by re-summing the 
BSV series [34] in case of a null hexagonal Wilson loop: we have used the saddle point method at 
large g to obtain the critical equations coinciding with the TBA equations of [39, 41, 40]. Then, 
we have computed the action on them and obtained the same (critical) Yang-Yang functional (or 
free energy) as in [39, 41, 40]. Interestingly, the same set-up should be easily applicable to the 
computation of the heptagon WL. Nevertheless, it would desirable to have a more direct under¬ 
standing of the phase we dubbed confinement of the fermions, which disappear as free particles, 
inside the mesons and their bound states. 
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A Functions 


This appendix is devoted to the introduction of the functions we used throughout the paper. 
In the study of scalars we found convenient to use the following shorthand notations: 


$(«) = <I>o(m) ^h{u) , (p{u, v) = (po{u -v) + (Ph{u, v) 


with 


*ho(M) = , ^h{u) =-iln 

t — 2u 


1 + ^ 


2x-{u)'^ 


(j)o{u-v) = iln ^.^^ ^ , fniu.v) =-2i 

l — u + v 


1 + 

9^ 

2x+{u) 


(1 

In 

V- 


2xr{u)x {v) 


i6{u, v) 


(A.l) 


(A.2) 


, (A.3) 


6{u,v) being the dressing phase [10] and x{u) = | 
also 


1 + \/l ~ 




2x 

, x^{u) = x{u ± |). We used 

(A.4) 
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For what concerns gluon and their bound states, we used the function 

Xiv, u\l) = xoiv - u\l + 1) + xh (v,u - + xh fn, m + ^ 


where 


/in il 2u 2u , , 1 I ^ 2x-(v)x(u) 

Xo(ulO = * In _ 2 ^ = 2 arctan — , xb(v, m) = * In ^ — 


(A.5) 


2 x'l'(v)x(u) 


which enjoys the expression 




2x-(v)x(^u—^^ 

1 - 


(A,6) 


2x+(v)x(u+^^ 

Scattering factors involving gluons and their bound states are expressed in terms of the function 


. , X .X I a;+(n) — x (u — ^) \ / 2x(u-^) 

x{v,u\l) + ^{v) = ill! I —^-77T-^^-^1+iln ^ 


X iv 


X (m + I) — a;“(n) 
Finally, for large fermions we introduced the function 

Xf(u,v) = Xo(u - n|l) + Xh{u,v) = iln 


X^[V) — 

x~^{u) — x{v) 
^x{v) — x~{u) 

Scattering factors involving large fermions depend on the function 

x'^{u)—x{v) f x~(u 


(A.7) 


2x(^u+ 2 j 


Xf{u, v) + <F(m) = i In ■ 


i In 


x(v} — X [u] 


x^[u, 


(A.8) 


(A.9) 


Scattering factors for small fermions are obtained from scattering factors for large fermions after 
the substitution 


Xf{u, v) + $(«) -Xh{u, v) =i\n 


1 - 


_ 1 _ ^/h) 

2x'^(u)x(v) \ 1 I x^(u) 

* = z In * 


2x~ (u)x(v) 


1 - 


x~ (u) 


where 


Xf(v) = 


2x(v) 


1- A 1 


V 

y2 


(A.IO) 


(A.ll) 


B Useful formulae 


B.l Fourier transforms 


We collect here some of the Fourier transforms 


b+cx) 


m = I due-^^-f{u) 


(B.l) 
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of functions f{u) we use in the main text. 
For scalar we used 


$o(w) = 


i -\-2u 


^h{u) = —fin 


i — 2u 

2x-(u)'-‘ 


1 + 


2x+{u) 


i n\ j -iku^ / ^ 27 r _ 1 M 

<Fo(/c) = / due = — ttC 2 , 

J-00 


^H{k) = ' 2 ' [1 - Jo{^/2gk)] 


(B.2) 

(B.3) 


and also 


(j)o{u — v) = flnT 


i + u — V 


0o(fc) = 


27re-l^l 

ik 


i — u + V 

ifo{u-v) = ^^(t)o{u-v) = ^ -- 

2ti dv TT 1 + (n — vY 

I- 


0o{k) = -e 


- _p-FI 


(B.4) 

(B.5) 


(j)H{u,v) = -2i 

^H{k,t) = 


In 


2x'^{u)x {v) 


1 - 


+ i6{u, v) 


2x (ii)x+ (hi) 
+00 /‘+00 


—00 J —00 

_ Ul + |fe| 00 

2 " 


SfTT 


>6 2 


+ 


due * w) = 

I I 
00 00 

sgn(f) EE Cr,r+H-2i/ {g){-lY+'^(^Jr-i{V2gk)Jr+2u{V2gt) 

r=2 1^=0 

(/y* _ 1 (\/2 ^ t) _j_ 2( a/ 2 5^/c) 


(B,6) 


We remark that in previous literature integral equations concerning the scalar sector are often 
written by using the ’magic kernel’ K [10], related to <pH by 


4’H{k, t) + ^nik, —t) = Sin^g^e 2 K(^^/2gk, \/2gt ), f, /c > 0 . 
For what concerns gluon bound states, we introduced 
Xo{u\l) = 2arctan^ ^ Xo{k\l) = ! due~'‘^'^Xo{u\l) = 

and for higher loops the function 

Xiv, u\l) = xoiv - u\l- + ^) + Xh (v,u - + Xh (v, “ + ^ > 


(B.7) 


(B.8) 


where 


Xh(v,u) = iln 


1 - 

® \ 

2x~{v)x{u) \ 

1 - 

9^ I 


(B,9) 


2a:+ {v)x{u) 
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whose Fourier transform reads 


^+oo 

' —OO 
+ CX) 


r'+cx) 


du 


dve-^'‘\-^^^x{v,u\l) = 27r6{t + k)—e-'^'^ + 


271 _|;.|1+1 


ik 


Jn{V2gk)Jn{V2gt). 

n=l ' ' 


(B.IO) 


In getting (B.IO) we used the Fourier transforms 


^+oo 


due 


—iku 


'V2 


277 


X \u 


±h) 


= ±n ( -— I 9{±k)—e^2’^Jn{V2gk). 


W 


(B.ll) 


It is useful to Fourier transform x(u,m|/) and xh{v,u) with respect only to the variable v: 

^ + (B.12) 


due-*^"x(u,u|/) = 


H-co 


+ ‘E 


9 


+ 00 


' 2 'Jn{V2gk) + i'^ 


9 


^ \V2ix{u-fj 

r^°° u 271 ifei f 

/ due "^^xh{v,u) =i—e 

n=l V 


—f xiu 
'1=1 ^ ^ ^ 

Finally, for large fermions we introduced the function 

Xf{v,u) = xoiv -u\l) + xh{v,u) = iln 


^ y\/ 2 i x(m+ f)y k 

^ ' Jn{V2gk) 


' 2 'J„(\/25(fc), 


(B.13) 


x~^{v) — x{u) 


^x{u) — X (u) / ’ 

whose Fourier transform with respect to v is easily extracted from (B.12, B.13): 


(B.14) 


+CX) 


9 


s 27r i»=i 27r |fc| 

dve-^^'’XF{v, u) = e-V + i—e" V 

-00 ik k ^^\^/2ix{u) 


JniyV^gk). (B.15) 


B.2 BES and BES-like integral equations 

The BES integral equation for the density dsEsik) in Fourier space reads as 

2i/c 6Hik,0) ik f dt ^ , 

asEsik) = - - - - + _g_|fc|^ J t)aBEs{t) • (B.16) 

Owing to the parity properties &BEs{k) = ^BEsi—k), we can restrict this equation in the region 
k > 0. Introducing the kernel K 

^H{k,t) + ^H{k^—t) = Sm‘^g‘^e ^ K{y/2gk, V^gt), f, A: > 0 , (B.17) 

we have 

dBEs{k) = . , K{V2gk,0) -r-rr / dte~^K{V2gk, V2gt)aBEs{t) ■ (B.18) 

Sinn Sinn ^ Jo 
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We can connect to qnantities used in [20] by means of 


^BEs{k) — — 


TT 


sinh I - 


7+(\/25'^) 


A: > 0 


(B.19) 


and 


7+(\/25'fc) = 2^2n7|„J2n(\/2fif/c), -ft{V2gk) = 2^(2n - l)'^'^^_^J2n-i{V2gk). (B.20) 


n>l 


n>l 


The total density at order Ins is (Jins(A;) = —87r6{k) + asEsik) which satisfy the equation 


^ins{k) = -8E6{k) 


ik 


^)(Tlns(^) - 


(B.21) 


4(1 — e 1^1) J TT^ 

The Fourier transform of the density associated to the hrst generalised scaling function [12] satisfies 
the equation 

ik /■+°° dt 


a^^\k) = 


TT 


[e ' 2 ' - Jo{y/2gk)] + 


sinh 


FI 


1 — e“l^l 


47J-2 


^H{k,t) 271 + 


(B.22) 


Eventually, the density ’all internal holes’, which satishes equation (3.8) of [21] with L = 3 is 
solution of 


27re“l^l 

a{k; x) = (cos kx — 1) + 


ik 


<•+00 


1 — e“l^l 


dt 


^H{k,t) 27r(cosfa; — 1) + (j(f; x) . (B.23) 


B.3 Integrals 

In the one loop case we make use of the following integrals 
dv ^ ib + V d, r(a + iv — iu) 

I i„ 


In 


and 


2e ib — V dv T{a — iv + iu) 


dv ^ ib — V + w d, ic — V + u 

— In —-— In-= * In 


r(a + 6 + m) 

^In . , a,6>0 

i{a + b — lu) 

(B.24) 

iic + b) -u + w , 

Y 1 M 1 ’ ^ c > 0 . 

i{c + b) + u — w 

(B.25) 


27r ib + V — w dv ic + v — u 
In order to show that (bound states of) gluons do not couple to (type b) isotopic roots, we used 
the following results 

f'^°° dv ^ x+(n) - x (u - f) d, ^ + v-u' x (u + "4) - x{u'- i) 

2e a;(M+|)— x“(n) dv — v + u' x{u'+ i) — x (u — 

f'^°° dv ^ X (u + 4) - x~^(v) d. I + n - u' x{u'- i) - x (u - 
“ ' — In \ ^ — In 4 -: = t In - 4-s -^^ + 


i In 


27r X-{v) - X [u - 4'j dv ^ - v + u' 

I > 2 


X [u+ — x{u' + i) 


u — u' + ^ — i u — u' — ^ 


u — u' + ^ u — u' — + i' 


(B.27) 


dv , x^{u)—x~^{v) d . I + n — u' xiu'— i) — x (u) 

— In —44 -r4 ^ In 2-= t In ^ ^ ^ + 

2e x~[v)—x~[u) dv ^ — v + u' x+[u) — x[u'+ 1 ) 


i In 


u + u' 


+ u — u' 


(B.28) 

































In calculations for the strong coupling limit of scattering factors, we used the following integrals 


'-1 

rl 


(^) * = 


u + 1 


1 7rsgn(i;) 


u — 1 

1 + M ^ 4 


1-1 z — V — 1 ’ 

/■^ 11 
I dzPV- 


1 — u 
v\ > 1, 


\u\ > 1 


\u\ < 1 


z-v v'l - 2 ;" 


= 0 , |w| < 1, 


'\w\>l 


2tt w — u 


tc + 1\ 


1 fl+lt 4_L /'l±“i 

2 \Z—1/ y/2 \l—u) 


dz 


w — z \w — 1/ z — u 

1 fz-l\^ 1 


'p-|>i 2vr-_^^ _ j.Xf{v) - x{z) V^ + l 


u — z 


(B.29) 

(B.30) 

(B.31) 

(B.32) 

|h| < 1, 1^1 > 1, (B.33) 

(B.34) 


feifly + (%(^) - I (i^)' + (rrl)'I + ^ (%(^) +1) v^irl I (i+l)' “ (i^l) 


in 
1 4 


2xf{u){u — v) 

C Collection of scattering factors 

C.l One loop: explicit expressions 

We list here the scattering factors at one loop: 

• Scalar - Scalar 


S'^^^\uh,Uh') = 


r (I - iuh) r (I + iuh') r(l + iuh - iuh') 


(C.l) 


r (I + iUh) r (I - iUh') r(l - iuh + iUh') 

Formula (C.l) does agree with result (3.8) of Basso-Belitsky [55], but seems to be the inverse of 
(2.13) of Dorey-Zhao [54]. 

• Gluon - Gluon 

r (1 + TO - to)) r (I - to) f (I + to) 


Sq^^\u,v)= /q \ /q \ 

° ^ ^ F(1-zm + to))F(|+to)F(|-to) 

In addition, we have 

So^\u,v) = S^o^^\u,v) = S^o^^\u,v) 


( 99 )/,, _ 0 ( 99 )/,, _ ci99)(,, v + l oigg)/^, ^ \ _ \c{gg)(^. „,m-i 


u — V — l 


-, Fr(«w) = [^r(^,«)]- 


(C.2) 


(C.3) 


(Large) Fermion - (Large) Fermion 


S'f^\u,v) = 


F(1 + TO — to) F(1 — iu) F(1 + to) 
F(1 — TO + to) F(1 + iu) F(1 — to) 


(C.4) 
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and when antifermions get involved 


q(FF) 

‘^0 


{u,v) = Sq 


(FF) 



{u,v). 


(C.5) 


• Gluon - Scalar 


S^^"\u,Uh) = = S^^"\u,Uh) = [S^"^\uh,u)]~^ = 

r (1 + m - iuh) r (^ + iUh) r (I - iu) 
r (1 - m + iuh)) r (i - iuh) r (I + iu) 


(C.6) 


(Large) Fermion - Scalar 


S^^‘'\u,Uh) = S^^‘'\u,Uh) = [S^''^\uh,u)] 

r(i + iu- iuh) r(l - iu) r(i + iuh) 
r(| - iu + iUh) r(l + iu) r(i - iuh) 


-1 




(C.7) 


Gluon - (Large) Fermion 


F (i + m — iv) F (I — iu) F (1 + iv) 
F (I — w + iv) F (I + iu) F (1 — iv) 


and 


So^\u,v) = [Fj^®^(n,M)] ^ = S^^^\u,v) = [Fj^^^(n,M)] ^ = S^^^\u,v) 


u — V — i/2 
u — V + i/2 


(G.8) 


(G.9) 








C.2 All loops: expressions in terms of solutions of integral equations 

We list here the factors found in [23]. We start from the ’direct’ S factors: 


v) 


v) 



S^^u,v) 


exp[—-?©(«, n)], 

dw d 

— [Xf{w,u) + ^{w)] — [xFiw,v) + $(u;)] 


exp< i 


2f' 


dw 


f . . , dw 

-expf-*Xo(M-w|2)+« / 


(C.IO) 

(C.ll) 

(C.12) 


d 


2 ^ [x(w,n|l) + <^>{w)] — [x{w,v\l) + $(«;)] - 


r (jin (jz (j"^ 


(C.13) 


u 


SlZ\u,v) = -exp{-ixo(M-'y|2)} Sl.fJ{u,v), 


u — V — l 

dw 


exp I* 




d 


— OQ 

dw dz 


2 ^ + ^{w)] — [x{w,v\l) + $(n;)] - 

d^ 


- j 27 r ^ <^{w)]j^Q{w,z)[x{z,v\l) + $( 2 ;)] \ , 

[S~^3{v,u)]-^ = SlfJ{u,v). 


(C.14) 

(C.15) 






The ’mixed’ S factors are: 




v) 


S^^^\u,v) 


v) 

S^3^\u,v) 


S^^^\u,v) 


= [S^^''\v,u)] ^ = exp|--i[xF(M,'y) + + 

+ i J ^^{u,w)[xf{w,v)+ ^{w)]^ , (C.16) 

= S^^^\u,v), S^^^\u,v) = S^^^\u,v), (C.17) 

= Gxp^^i[x{v,u\l) + <^{v)]-i j ^^{v,w)[x{w,u\l) + <^{w)]^ , (C.18) 

= = -exp|-ixo(M - w|l) + 

r+°° dw d 

+ ij -^[xiw,u\l) + ^{w)] — [xF{w,v) + ^{w)]- 

r dw dz 

- i J ^^[x{w,u\l) + ^{w)]-^^e{w,z)[xF{z,v) + ^z)]j =, (C.19) 

= 1 ^lld\u,v) = -exp{-ixo(M-w|l)} Sll^\u,v) 

U V ^ 

= S^'^^\u,v), (C.20) 

= [S^^^\v,u)]~^ = exp^^i j ^[x{w,u\1) + <^{w)]-^[xf{w,v) + ^{w)]- 

- i J ^^[x{w,u\1) + ^{w)]^^^Q{w,z)[xf{z,v) + ^{z)]^ , (C.21) 

= (C.22) 


The S matrices involving small fermions are obtained from the corresponding ones for large fermions 
by means of the replacement 


Xf{,v,u) + ^{v) —)■ -xh{v,u). (C.23) 

All the scalar factors are expressed in terms of known functions listed in Appendix A and the 
’dynamical’ function Q{u,v) [23], which equals 

Q{u,v) = Q'{u,v) + P{v), (C.24) 

where Q'{u,v) and P{v) are found as solutions of the linear integral equations 


“+00 


Q'{u,v) = 0(m, n) + $(«) — / dw(p{u,w)Q'{w,v), 


P{v) = —<h(n) 




^[(p{v, w) - ip{v, -w)]P{w). 


(C.25) 

(C.26) 


Scattering factors involving bound states of gluons 
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If bound states of gluons are present, the factors involving the gauge held should be generalised 
as follows. 

The scattering factor between a bound state of gluons with length m and center u and a scalar 

is 


— exp|i[x(^',«|m) +‘l>(t;)] — i J ■^-^(v,w)[x{w,u\m) + . (C.27) 

The scattering factor between a bound state of gluons with length m and center u and large 
fermions is 

=-expJ^-ixo{u-v\m) + 




+ i I ^[x{w,u\m) + (!>{w)]-^^[xf{w,v) + 


271 271 

{u,v) = (n, m)] "^ = exp |i 


dw 

i [ ^^[x{w,u\m) + ^{w)]-^^^e{w,z)[xF{z,v) + ^{z)]^ , 


(C.28) 


(C.29) 


f+OO 


^^[x{w,u\m) + ^{w)]^[xf{w,v) + $(«;)] 


2f 


- i J ^^[x{w,u\m) + ^{w)]-^^-^Q{w,z)[xF{z,v) + ^{z)]^ , 


(C.30) 
(C.31) 

The replacement xf{w,v) + <h(tn) —)> —xh{w,v) gives the corresponding quantities for small 
fermions. 

The scattering factors between a bound state of gluons with length m and center u and a bound 
state of gluons with length I and center v are 


= exp^^-ix{u,v\mj)+i I ^[x{w,u\m) + ^{w)]^[x{w,v\l) + ^{w)]- 

(C.32) 


f —oo ^ ‘ dw 

i I ^^[xiw,u\m) + ^{w)]^^^e{w,z)[xiz,v\l) + ^z)]^ , 


r» _|_QQ 7 7 

= expj^y ■^lx(w,ujm)+ ^(w)] — lx(w,vjl)+ ^(w)]- 

- I I ^^[x{w,u\m)+ <^{w)]-^e{w,z)[x{z,v\l)+ ^{z)]\ , 


2f 2f 


dwdz 


where 


m—1 


x{u, v\m, 1) = Xo(w — n|m + /) — Xo{u — v\m — 1) + 2 xo(w — n|m + / — 27 ) 

7=1 


(C.33) 


(C.34) 
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In the particular case m = / = 1, since = —1, one recovers from (C.32) the gluon-gluon 

scattering factor (C.13). 

Remark The following relations hold, for u real and < 2g‘^, 

m—1 

2 

x{v,u\m) + ^{v) = ^ [x('y,M + */|l) + *h(r’)], m>l, (C.35) 

7 m — 1 

2 


where + *^|1) has to be understood as analytical continuation of x(n,M|l). 

To prove this statement, we refer to (A.7) and remember the following properties 

lim xiu — ie) = lim —-- , 

e^o+ £^■ 0 + 2x[u + le) 


(C.36) 


which are valid for u real and v? < 2g‘^. Therefore, when the complex variable u crosses the real 
axis in the region —\/2g < Reu < \/2g, the function x{u) is analytically continued in g^j2x{u). 
With the help of this property, relation (C.35) is easily shown. 


C.3 Strong coupling and mirror in hyperbolic rapidities 


In the strong coupling perturbative regime, the scattering matrices for gluons, fermions and mesons 
can be suitably recast in term of hyperbolic rapidities, according to the following identities (written 
up to OiXjg^') corrections): 

Gluons: 

vF = \f2gu^ , = tanh(20), (C.37) 

Fermions: 


Uf = \/2guf , Uf = coth.{26) or else 2xf = ta.nh.6, 


Mesons: 


Explicitly, we obtain: 


um = y/2guM , um = coth(26'). 


^( 99 ) ( 0 , 0 ') = 


exp 


V2g 

^( 99 ) ( 0 , 0 ') = (i 
^(//)(0,0') =exp 


+ 


cosh 29 cosh 29' 


O 


tanh26* — tanh20' 2sinh(0 —0') J \g^ 


S^>^^\9,9f) = exp 


2^/2g sinh(20 — 20 

i 2 cosh(0j — 0) -|- y/2 


+ 0 


Ag tanh20 — coth20/ \g 


(C.38) 

(C.39) 

(C.40) 

(C.41) 

(C.42) 

(C.43) 
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s^sf)^e,ef) = exp 


i 2 cosh( 6 *j — 9) — \/2 (9 ^ ^ I 

Ag tanh 20 — coth 20 / \g‘^) J 


S(^^){e,e') =exp 


i sinh 26 sinh 26' 
y/2g sinh( 6 ' - 6') 


(C.44) 

(C.45) 


^(^9)(0,0') = =exp 


i cosh( 6 * — 6') 
g tanh 26' — coth 26 


(C.46) 


Mirror transformations: 

The mirror rotation should be implemented in different ways on the scattering phases, depending 
on the kind of particle the transformation is acting on. For instance, in the scalar case it is achieved 
by means of a shift u —?■ u'^ = u + i- For gluons, the mirror transform is performed via a closed 
path across the complex rapidity plane = u), passing through a cut, so that actually the initial 
and hnal points do not he on the same sheet. Dehning a procedure for the mirror rotation on 
fermions is more involved, and for this purpose we refer to [36]. A complete all coupling study of 
all the mirror S matrices can be found in [ 66 ]. Nevertheless, as long as the perturbative strong 
coupling regime is concerned, the mirror rotation gets simplihed: for gluons, scalars and mesons it 
amounts to an imaginary shift in the hyperbolic rapidities 6'^ = 6 + i^- For instance, we have: 


S<l9»,«')=S<»'(« + i^.«') = exp 


i 

Wg 


1 

tanh 26 — tanh 26' 


cosh 29 cosh 29' 
2i cosh( 6 * — 9') 


+ 0 





(C.47) 




2i 


exp 


y/2g tanh 26* — tanh 20 ' 
cosh 26 cosh 26' 


oil 


^/2g tanh 26 — tanh 26' 2i cosh (0 — 6') 


+ 0^T2 


(C.48) 


Si^l\6,6')=exp 


i sinh 26 sinh 26' ,. , 


r- w /i- ; 7 T (* sinh(0 — 6') — 2) + O [ — 

2^/2g sinh (20 - 20 ') ^ ^ \g^ 

j 1 sinh 20 sinh 20 ' 


s'"j {»,»') = exp 


i 4:y/2g cosh (0 - 0 ') 


oil 


n{MM) 

^mir 


( 0 , 6') = exp 


1 sinh 20 sinh 20 ' 
^/2g cosh (0 - 0 ') 




o[Mg) 

^mir 


( 0 , 0 ') 


c{Mg) 

^mir 


( 0 , 6') = exp 


1 sinh (0 — 0 ') 
g tanh 20 ' — coth 20 


(C.49) 

(C.50) 

(C.51) 

(C.52) 
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C.4 On the factor 0 


Since the scalar-scalar factor 0 is the building block for all the scattering factors, we give here 
some alternative constructions for it. Following what was done in [23], we dehne the function 


M{u,v) = 


0(m, v) + 0(m, —v) 


(C.53) 


which stores all the information on 0, since (see (2.23) of [23]) 

Q{u,v) = M{u,v) — M{v,u). 


(C.54) 


Then, we use (2.21) of [23] to express M{u,v) in terms of densities and a{u;v). Formulae 

(3.7, 3.8) and Neumann expansions (3.9) allow to arrive at 


1 T{1 + lu + iv)T{l + lu - iv) rU-tu) 

M(u,v) = -ln-:0;-^^-r4 -7 dn —^ {- + 

2 T{1 — tu + tv)T{l — tu — tv) T + iu) 


(C.55) 


OO OO 


n=0 p=l \ V / 


X (u + |(1 2n)) 


x{u- |(1 2n)) 


The hrst line of (C.55) is the one loop contribution; the second line is the higher than one loop 
correction. We can manipulate (C.55) in order to get alternative expressions for M{u,v): for 
instance. 


t T{1 + lu + iv)T{l + tu - tv) rU-tu) 

M(u,v) = -ln-0-(-0-(- -Fdn—0f -i- 

2 r{l — tu + tv)r{l — tu — tv) F (I-I-to) 


(C.56) 


El'S!" + 

p=l 


T -iV2g sine-iu) 
-n 27ri F (I — z\/25fsin6 *-|-to) 


and 


* T{1 + tu + tv)T{l + tu — tv) T(^ — tu) 

M(u,v) = -ln-0-00-(- + ^ln—0-f + 

2 r{l — tu + tv)r{l — tu — tv) F -|-in) 


(C.57) 


EEi 


- 1 ) 


n+p 


OT!(n + p) 


/ n r 

! (^) r 


{2n+p-l) 


- tu 


- 0(2n+p-l) 


D Some calculations on how to get scattering factors 

D.l Scalars vs others 

The ’fermionic’ contribution to Z 4 ^{v) — P reads 

Z^{v) - P\f = F^{v,UF,j) + j ^-^^{w)F^{w,ufj) . (D.l) 
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Inserting the integral equation (2.63) for we arrive at 


/ 7 7 

'^'dw ^ • (D-2) 

Remembering the integral equation (2.62) for 0' and then, after an integration by parts, the 
equation (2.63) for F^, we get 

/ dll) d 

■^J^Q'{w,v)[xF{w,UF,j) +<h(w)] (D.3) 

= XF{v,UF,j) + <h(u) - j ^-^Q{v,w)[xF{w,UF,j) + $(w)] = f In (u, Mpj) . 

The same calculation can be done for large antifermions. For gluons we can repeat the same 
reasonings after the substitution xf{v,u) —)■ x(n,M|l). For small (anti)fermions we perform the 
substitution xf{v,u) + $(n) —?■ —xh{v,u). 

For scalar-scalar factor we have to consider 

Z4^{v) - P\s = Q'{v,Uh) - ^{Uh) + j ^-^^{w)Q'{w,Uh). (D.4) 

Now, we remember the integral equation (2.62) for 0', which allows to write the identity 

/ dw d ^, f dw (p(uh,w) — (p(uh,—w) X 

__4(x„)e'(u.. u,) = -j T (^). (D.5) 

Therefore, using the integral equation (2.67) for P we get 

Z^{v) - P\s = 0'(n, Uh) + P{uh) = 0(n, Uh) = i ln[-S'^'*®)(n, Uh )]. (D.6) 


D.2 Non scalars vs others 


Taking as a prototypical example the case of fermions, after multiplication by e®'^ we have to cope 
with 

/d'7) fj 

— [XF{v,UF,k) + ^{v)]j^Z4{v) , (D.7) 

which is the master relation from which the various scattering factors originate. 

For instance, the fermion-fermion factor is 


exp 


/ dij d 

— [XF{v,UF,k) + ^v)]—F^{v,UF,j) 


Mixing equations (2.62) and (2.63), we get 
F^{v,u) = xf{v,u) + <l>(n) 


dw d 
271 dw 


Q'{v,w)[xf{w,u) + ^{w )], 


(D.8) 
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(D.9) 







which, inserted in (D. 8 ) gives the fermion-fermion factor as reported in Appendix C: 


exp< i 


dw d 

■[XFiu!,UF,k) + ^{w)] — [xF{u!,UF,j) + $(w)] - 


J —oo 

dw dx 


2n' 


d^ 


For what concerns the fermion-scalar factor, it receives contribution from 
exp 

Expression (D.ll) has to be multiplied to the factor 


(D.IO) 


i J ^[XF{v,UF,k) + ^{v)]-^Q'{v,Uh) 

= exp 

i J ^[XF{v,UF,k) + ^{v)]-^e{v,Uh) 



(D 


XF,k _ ^iXF{uh,up^k) 

xt - XF,k 


(D.12) 


present in equation (2.80) and to the factor present in to get the full fermion-scalar 

factor. 


E Bethe equations 


We list the complete set of Bethe Ansatz equations we found in this paper: 
• Scalars 




Kb _ \ i ^ 


i=i 


'^h ^b,j 2 


[Uh,Ujj 


h'=l 

h'^h 

Nf 


i=i 
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i=i 


■ n UF,) n ^F,) n «/,.) n 

j=i j=i j=i j=i 

Large fermions 


(E.l) 
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:KF(-’)(.,,.)+2iD(-)(„,,.) TT -fF.t - TT 2 , 

^,UF,k-Ua,j-t2^^ 

7 = 1 ’ ' h=l 


Nf 


JJ S^^^\uF,k, Ufj) JJ S^^^\uF,k, Upj) U S^^^\uF,k,Uf,j) JJ S^^^\uF,k, Uf,j) 

j=l j=l j=l j=l 


N„ 


n s^^<^\uF,k, Mf) n s^^~^\uF,k, wf) 

i=i i=i 


(E.2) 
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Large antifermions 


1 _ TT n,.i + i/2 yj . 

2=1 - Uc,j - 112 


Nf 


Nf 


Nf _ _ _ 

YlS^^^\uF^k,UF,j)'^S^^^\uF^k,UFj)Y[S^^^\uF^k,Uf,j)Y[S^^^~\uF^k,Ufj] 


1=1 

N„ 


1=1 


1=1 


1=1 


Y[S^^^\up^k,u’^j)YlS^^^\up^k,u^j) 

1=1 1=1 


(E^3) 
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Small fermions 

1=1 


^a,j ~F i l2 
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u 


a,3 


i/2 J-J- 

' h=i 




Ni, 


Nf 


Nf 


1=1 

Ng N-g 

1=1 1=1 

Small antifermions 


n wfj) n Ufj) jj Ufj) 

1=1 1=1 1=1 
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1 = JJ 


^f,k ^c,j i/2 
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IVf 


IVs 


Aff 


Aff 
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Gluons 


Af„ 


Afs 


H 


I ^ giFp(«)K)+ 2 iF(«)K) jj s^^^\ul,u^j)l[S^^s\ulu^)l[S^^^\uluh)- 


i=i,i^fc 


1=1 


A^f 


Affi 


Nf 


h=l 

Nf 


■ n “F,i) n WF,i) n “/.i) n w/g; 

1=1 1=1 1=1 1=1 


(E.6) 


Barred gluons 


N„ 


N-a 
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1 ^ ^iRP^^Hul)+2iD(oHul)Y[S^~^^\ulu^.) JJ S^~^~^\ulu^)l[S^~^^\uluh)- 
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Nf Np Nf 
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Isotopic roots 

Np 


n Ua,k — UF,j + TT Ua,k ~ Ufj ^ ij^l 
1! . - 1! T-^ ■ - ? /9 A A 


n 


Ua,k ^CL,j “1“ i 1 r ^a,fc ^b,j i/2 


. TT ~ ~ (E.8) 

a,j i ■_^ 'Ua,k ^b,j 


j—\ i/2 ^a,fc ^f,j i/2 ^a,k ^ 

TT Ub^k — Uh + i/2 ^ -i-i- Ub^k ~ Uaj — i/2 -i-i- Ub^k ~ Ucj — i/2 -i-i- Ub^k ~ Ubj + i (^E 9) 

Ub,k '^h i/2 Ub,k '^a,,j “1“ i/2 Ub,k “1“ i/2 Ub,k ^b,j i 

jvfc 






TT j + i/2 -p-|- Mc,fc Ujj + i/2 ^ -p-|- Uc,k ~ Uc,j + i tt Uc,k ~ Ub,j — i/2 

Z^Cjfc ^F,j i/2 ._^ Uc^k ^f,j i/2 ^ Uc^k ^c,j i '^c,k ^b,j i/2 


. iZ 6 
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